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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 
typed using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 

H. Jeffreys, of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
1.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I.A.U. practice. 


4 ams.—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destro 


wan Tables.—These should be arranged so that they can be printed upright on 
page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—When papers are submitted in revised form it is especiall 
requested that they be accompanied by the original MS. r 


Reading of Papers at Meetings 
8. When submitting papers authors are requested to indicate whether they will be 


willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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THE PHILOSOPHICAL MAGAZINE 


A JOURNAL OF THEORETICAL, EXPERIMENTAL AND APPLIED PHYSICS 
The Philosophical Magazine is an independently owned journal for the publication of 
papers in Natural Philosophy; under this term the Editors include experimental and 
theoretical physics, theoretical mechanics and certain branches of applied mathematics. 
This important monthly publication has a world-wide circulation, and has been 
published continuously since 1798. 


Annual Subscription £8, post free Price per part 15s., plus postage 


ADVANCES IN PHYSICS 


A QUARTERLY SUPPLEMENT OF THE ‘ PHILOSOPHICAL MAGAZINE * 
On Ist January, 1952, the Philosophical Magazine published the first number of a new 
quarterly supplement. The aim of this supplement will be to give those interested in 
physics comprehensive and authoritative accounts of recent important developments. 
It is felt by the Editor that, in view of the rapid advances in many branches of physics, 
scientists will welcome a journal devoted to articles of this type. Contents for Vol. 1, 
1952, and Vol. 2, 1953, are now available. 


Annual Subscription £2 15s., post free Price per part 15s., plus postage 


ATOMIC SCIENTISTS’ JOURNAL 
(Former title * Atomic Scientists News ’) 
THE OFFICIAL PUBLICATION OF THE ATOMIC SCIENTISTS’ ASSOCIATION 

The Journal contains articles by acknowledged authorities on various aspects of atomic 
energy, and provides a medium in which atomic scientists can discuss their work and 
its wider implications. It also provides a channel of communication between the 
specialists and the general public. 

Annual Subscription £1 12s. 6d., post free Price per part 6s., plus postage 


Payment from overseas should be made, for preference, by Foreign Money Order 
obtainable from Post Offices. 

Cheques drawn on London Banks or dollar cheques drawn on American Banks 
are also acceptable. 


' 


Taylor & Francis Publications 


ANNALS OF SCIENCE 


A QUARTERLY REVIEW OF THE HISTORY OF SCIENCE 


Founded in 1936, Annals of Science has come to occupy a secure place in the esteem of 
the historian of science as well as in the mind of everyone interested in the history of 
science in modern culture. It is recognised as an indispensable part of modern culture. 
Each volume consists of approximately 350 pages in royal octavo, copiously illustrated 
with half-tone plates and text-figures. 


Annual Subscription £3 3s., post free Price per part 17s. 6d., plus postage 


THE SCIENTIFIC WORK OF RENE DESCARTES (1596-1650) 


By J. F. SCOTT, B.A., M.SC., PH.D. 

This book puts the chief mathematical and physical discoveries of Descartes in an 
accessible form, and fills an outstanding gap upon the shelf devoted to the history of 
philosophy and science. The careful treatment that Dr. Scott has accorded to this 
work of Descartes is welcome, and will be an asset to all libraries. Publication is 
recommended and approved by the Publication Fund Committee of the University of 
London. 

212 pages, 7 in. x 10 in. Price £1 net 


NOTES ON PRACTICAL PHYSICS 


By PROFESSOR G. STEAD, M.A., D.SC., F.INST.P. 
Professor of Physics, Guy’s Hospital Medical School, University of London 
The book provides enough information to enable the average student to start an 
experiment without waiting for the assistance of a demonstrator. 
Demy 8vo. 109 pages, 98 text-figures Price 7s. 6d., plus postage 


GLAISHER’S HYGROMETRICAL TABLES 


SEVENTEENTH EDITION 
Adapted to the use of the Wet and Dry Bulb Thermometer. 
Price 12s. 6d., plus postage 


DIARY OF ROBERT HOOKE, M.A., M.D., F.R.S., 1672-1680 


Epirep By H. W. ROBINSON AND W. ADAMS 
** This vivid record of the scientific, artistic and social activities of a remarkable man 
during remarkable years has too long remained in obscurity.’”-—Extract from foreword 
by Sir Frederick Gowland Hopkins, O.M., President of the Royal Society. 
Recommended for publication by the Royal Society, London 
First Published 1935 Price £1 5s. net 
Limited stock still available 
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Taylor & Francis Publications 


MATHEMATICAL WORK OF JOHN WALLIS, D.D., F.R.S. 


By J. F. SCoTT, PH.D., B.A. 

~’ His work will be indispensable to those interested in the early history of The Royal 
Society. 1 commend to all students of the Seventeenth Century, whether scientific or 
humane, this learned and lucid book.”-—Extract from foreword by Professor E. N. da C. 

Andrade, D.Sc., Ph.D.. F.R.S. 

Recommended for publication by the University of London 

Price 12s. 6d. net 
Limited stock still available 


CORRESPONDENCE AND PAPERS OF EDMOND HALLEY 


ARRANGED AND EDITED BY EUGENE FAIRFIELD MACPIKE 
First published on behalf of The History of Science Society by Oxford University Press. 
Now re-issued by Taylor & Francis, Ltd. 
Price £1 Is. net 
Limited stock still available 


MEMOIRS OF SIR ISAAC NEWTON’S $ LIFE 


By WILLIAM STUKELEY, M.D., F.R.S. (1752) 
This work is taken from an original manuscript now in the possession of the Royal 
Society, London. 
First Published 1936 Price 5s. net 
Limited stock still available 


HEVELIUS, FLAMSTEED AND HALLEY 


By EUGENE FAIRFIELD MACPIKE 
This work contains an account of three contemporary astronomers and their mutual 
relations. 


Published by arrangement with The History of Science Society 
First Published 1937 Price 12s. 6d. net 
Limited stock still available 


THE JOURNALS AND MONOGRAPHS OF THE 
QUEKETT MICROSCOPICAL CLUB 


Now printed and published by Taylor & Francis, I td. Full details will be sent on 
request. 
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incorporating ‘ The Journal of Botany ’ 


Probably few independent scientific journals have such a long history as the * Annals ’. 
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Natural History and Journal of Zoology, Botany, Mineralogy, Geology and Meteorology) 
and one from 1838 (The Annals of Natural History, or Magazine of Zoology, Botany and 
Geology), it has appeared continuously since 1840 under its present title: The Annals 
& Magazine of Natural History, including Zoology, Botany and Geology. 


Annual Subscription £6, post free Price per part 10s., plus postage 


THE JOURNAL OF BOTANY 
British and Foreign 


Epitep By A. B. RENDLE, D.SC., F.R.S., F.L.S. 
Late Keeper, Department of Botany, British Museum (Nat. Hist.) 
All parts of this important monthly botanical journal are available up to the year 1942. 
Demy 8vo. Price per part 2s. 
First Published 1877 


THE NIDIFICATION OF BIRDS OF THE INDIAN EMPIRE 


By E. C. STUART BAKER, C.LE., O.B.E., F.Z.S., ETC. 


In the four volumes of this work the author has collated the information recorded by 
previous observers, but his personal experience, gathered during thirty years’ residence 
and field-work in the Indian Empire, has enabled him to add a wealth of entirely new 

matter as well as to verify many points which had been matters of surmise or conjecture. 
Publication of this work brought the knowledge of the nidification of Indian birds 
up to date. Each volume contains approximately 500 pages and several plates. 
Medium 8vo., cloth boards, gold lettered. 


Price complete £5, post free Price per volume £1 10s., plus postage 
First Published 1932 Limited stock still available 


A HANDBOOK OF THE BIRDS OF EASTERN CHINA 
(Chihli, Shantung, Kiangsu, Anhwei, Kiangsi, Chekiang, Fohkien & Kwangtung Provinces) 


By J. D. D. LA Toucue, C.M.Z.S., C.F.A.0.U., M.B.O.U., ETC. 


This work, published in two volumes, contains a full description, with synonyms and 
distribution, of over 750 varieties of birds. Illustrated with plates. Volume I in five 
parts; volume II in six parts. 


Medium 8vo. Price per part 7s. 6d. net 
First Published 1925 Limited stock still available 


BRITISH AND IRISH BOTANISTS 


SECOND EDITION 
COMPILED BY JAMES BRITTEN, F.L.S., AND GEORGE R. BOULENGER, F.L.S., F.G.S. 
REVISED AND COMPLETED BY A. B. RENDLE, M.A., D.SC., F.R.S. 
Formerly Keeper of Botany, British Museum 

A comprehensive list, including all who have contributed to the literature of botany. 
The growth of the work may be estimated by comparison of the original edition 
comprising 246 pages, with the present volume of 338 pages. 

Demy 8vo.. cloth boards, gilt lettered Price 16s., post free 
First Published 1931 Limited stock still available 
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THE IBIS 
THE QUARTERLY JOURNAL OF THE BRITISH ORNITHOLOGISTS’ UNION 

The * Ibis’ deals with the whole field of ornithology. In addition to work on the 
ecology, behaviour and taxonomy of individual species and groups, a special effort is 
made to secure full authoritative contributions on general subjects. At intervals numbers 
are devoted entirely to particular subjects of widespread interest, such as migration. 
At present the annual volume of the * Ibis’ contains 500-600 pages of origina! matter. 
Each issue contains a large number of notices of recent ornithological literature. 
The no ices are classified by subjects. are critical. and are designed to keep members 
informed of recent work in all branches of ornithology throughout the world. They 
are so drafted as to present the important features of each publication. 


Annual Subscription £3 10s. 0d., post free Price per part 17s. 6d., plus postage 


A HANDBOOK OF THE BIRDS OF ICELAND 


By Masa U. HACHISUKA, F.Z.S., ETC. 
Member of the British Ornithologists’ Union and of the Ornithological Society of Japan 
This most valuable work on the ornithology of Iceland contains descriptions of 162 
varieties of birds. Illustrated with seven plates. 
Super-Roy. 8vo., cloth boards, vi 128 pages Price 12s. 6d. net 
First Published 1918 Limited stock still available 


A VERTEBRATE FAUNA OF THE MALAY PENINSULA 
from the Isthmus of Kra to Singapore (including the adjacent islands) 
Epitep By H. C. ROBINSON, C.M.Z.S. 

Director of Museums, Federated Malay States 


Reptilia and Batrachia 


By GEORGE A. BOULENGER, D.SC., PH.D., F.R.S. 
Published under the authority of the Government of the Federated Malay States. 
This work comprises xiii-+ 294 pages with 79 text-figures and a map. 
Medium 8vo., cloth boards, gold lettered Price 15s. net 
First Published 1912 Limited stock still available 


THE FLORA OF THE PRESIDENCY OF BOMBAY 


By THEODORE COOKE, C.I.E., M.A., M.A.L., LL.D., F.L.S., F.G.S., M.INST.C.E.L 
Formerly Principal of the College of Science at Poona and Director of the Botanical Survey 
of Western India 
Published, in two volumes, under the authority of the Secretary of State for India in 

Council. Only volume II, which is printed in five parts, is now available. 
Roy. 8vo. Price: Parts I and II 9s. each, Parts III and IV 8s. each, Part V 12s. 
Limited stock still available 


Tayler & Francis Publications 


THE FAUNA OF BRITISH INDIA 
including Ceylon and Burma 
11ST OF VOLUMES STILL AVAILABLE 
VERTEBRATA 
Birds 
SEcoND Epition. By E. C. STUART BAKER 
Vol. I. [Passeres, Fam. I. Corvidae—VIII. Troglodytidae]. Pp. i-xxiii, —_ 8 col. pls., 


text-figs. 30/-. Aug. 24, 1922. 
Vol. II. [Passeres, Fam. IX. Cinclidae—XVII. Regulidae|, Pp. i-xxiii, 1-561, 8 col. pls.. 

text-figs. 30/-. April 30, 1924. 
Vol. Ill. [Passeres, Fam. XVIII. Irenidae—XXXIII. Eurylaimidae]. Pp. i-xx, 1-489, 

7 col. pls., map, text-figs. 30/-. March 20, 1926. 


Vol. IV. [Coraciiformes]. Pp. i-xxiv, 1-471, 7 col. pls., text-figs. 30/-. July 28, 1927. 
Vol. V. Accipitres, Columbae, Pterocletes, Gallinae, Hemipodii|. Pp. i-xviii, 1-469, 
col. pls., text-figs. 30/-. March 21, 1928. 
VI. (Grallae, Charadriiformes, Steganopodes, Tubinares, Herodiones, Phoeni- 
copteri, Anseres, Pygopodes]. Pp. i-xxv, 1-499, 3 pls., text-figs. 30/- 
March 26, 1929. 
VII. [Synonymical Catalogue, Passeres—Grallae|. Pp. i-vili, 1-484. 30/-. 
March 30, 1930. 
Vol. VIII. [Synonymical Catalogue, Grallae—Pygopodes: Corrigenda and Addenda ; 
Index]. Pp. i-iv, 485-801. 15/-. Sept. 25, 1930 


Reptilia and Batrachia 
SECOND EDITION. By MALCOLM A. SMITH 


J. {Loricata, Testudines]. Pp. i-xxviii, 1-185, 2 pls., map, text-figs. 15/-. 
March 27, 1931. 


Vol. 


Vol. 


Vol. 


Vol. II. [Sauria]. Pp. i-ix, 1-440, 1 pl., 2 maps, text-figs. 30/-. Feb. 7, 1935. 
ARTHROPODA 
Lepidoptera 
MOTHS 
Vol. Ill. |Fam. 26, Noctuidae (Subfam. Focillinae, Deltoidinae); 27, Epicopiidae: 


28. Uraniidae; 29, Epiplemidae; 30, Geometridae]. By G. F. HAMPSON. 


> i-xxviii, 1-546, text-figs. 28/-. Feb. 21, 1895. 

Vol. V [Sphingidae]. By R. D. BeLt and F. B. Scott. Pp. i-xviii. 1-537, 15 pls.. 

text-figs. 32/6. June 15, 1937. 
BUTTERFLIES 


SECOND EprTIon. By G. TALBOT 


Vol. I. {Papilionidae, Pieridae]. Pp. i-xxix, 1-600, 3 pls., map, text-figs. 35/-. 
March 8, 1939. 
Vol. II. {[Danaidae, Satyridae, Amathusiidae, Acraeidae]. Pp. i-xv, 1-506, 2 col. pls., 
map, text-figs. 55/-. Dec. 31, 1947 

Coleoptera 
ADEPHAGA 
CARABIDAE. By H. E. ANDREWES 
Vol. I. (Carabinae}. Pp. i-xviii, 1-431, 9 pls., text-figs. 22/6. May 15, 1929. 
Vol. I. (Harpalinae—I}. Pp. i-xvi, 1-323, 5 pls., map, text-figs. 22/6. Oct. 23, 1935. 
STAPHYLINOIDEA 
STAPHYLINIDAE. By MALCOLM CAMERON 

Vol. 1. [(Subfam. Micropeplinae, Oxytelinae, Oxyporinae, Megalopinae, Steninae, 


Enaesthetinael. Pp. i-xvii, 1-471, 3 pls., map, text-figs. 30/-. March 31, 1930. 


Vol. IT. [Subfam. Paederinae]. Pp. i-viii, 1-257, 2 col. a text-figs. 15/-. 
Feb. 28, 1931. 
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Vol. {Subfam. Staphylininae, Trichophyinae, Termitodiscinae, Pygosteninae, 

Wokerotenet Pp. i-xiii, 1-443, 4 col. pls., text-figs. 30/-. March 30, 1932. 

Vol. IV. Part I. [Subfam. Pseudoperinthinae and Aleocharinae (part)}]. Pp. i—xviii, 

1-410. map, text-figs. 25/-. Aug. 11, 1939. 
Part II. [Aleocharinae}]. Pp. 411-691, 3 col. pls., map, text-figs. 15/-. 

Aug. 11, 1939. 


CLAVICORNIA 


ee Languriidae and Endomychidae. By G. J. Arrow. Pp. i-xvi, 1-416, 
1 col. pl., map, text-figs. 30/-. March 21, 1925. 


PHYTOPHAGA 


CHRYSOMELIDAE. By S. MAULIK 
Vol. If. [Hispinae and Cassidinae}. Pp. i-xi, 1-439, text-figs. 21/-. Aug. 9, 1919, 
Vol. III. [(Chrysomelinae and Halticinae}. Pp. i-xiv, 1-442, map, 25/-. 
May 20, 1926. 
Vol. IV. [Galerucinae|. Pp. i-xvi, 1-648, 1 col. pl., map, text-figs. 35/-. 
Jan. 30, 1936, 


RHYNCHOPHORA 


Curculionidae |Part |. Brachyderinae, Otiorrhynchinae|. By Guy A. K. MARSHALL. 
Pp. i-xv, 1-367, text-figs. 21/-. Nov. 28, 1916. 


LAMELLICORNIA 


SCARABAEIDAE. By G. J. ARROW 
Part II. [Rutelinae, Desmonycinae, Euchirinae]. Pp. i-xiii, 1-387, 5 pls., text-figs. 
21/-. May 6, 1917. 
Part III. (Coprinae|. Pp. i-xii, 1-428, 13 pls., map, text-figs. 30/-. Dec. 15, 1931. 
Part IV. [Lucanidae and Passalidae|. Ppp. i-xi, 1-274, 23 pls. (1 col.). 50/-. 
April 28, 1950. 


Diptera 


Vol. LL. [Brachycera}, Vol. I. |Stratiomyiidae, Leptidae, Nemestrinidae, Cyrtidae, 
Bombyliidae, Therevidae, Scenopinidae, Mydaidae, Empidae, Lonchopteridae. 
Platypezidae|. By E. BRUNETTI. Pp. i-ix, 1-401, 4 pls., text-figs. 35/-. 

May 28, 1920. 

Vol. IIL. [Pipunculidae, Syrphidae, Conopidae, Oestridae|. By E. BRUNETTI. Pp. i-xii, 
1-424, 6 pls., text-figs. 30/-. March 1, 1923. 

Vol. V. [Culicidae (Tribes Megarhinini and Culicini)|. By P. J. BARRAUD. 
Pp. i-xxvii, 1-463, 8 pls., text-figs. 30/-. March 14, 1934. 

Vol. VI. (Calliphoridae]. By R. SENior-WxHiTe, DAPHNE AUBERTIN and J. SMART. 
Pp. i-xiii, 1-288, map, text-figs. 18/-. March 28, 1940. 


Rhynchota 
By W. L. Distant 
Vol. VI. |[Homoptera: Appendix (Cicadidae, Fulgoridae, Membracidae, Cercopidae, 
Jassidae (pt.))]. Pp. i-viii, 1-248, text-figs. 14/-. March 31, 1916. 
Vol. VII. [Homoptera: Appendix (Jassidae (pt.)); ge Addenda (Penta- 
tomidae, Coreidae, Berytidae, Lygaeidae)}. Pp. i—viii, 1-210, text-figs. 
May 9, 1918. 
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Orthoptera 
Acridiidae. By W. F. Kirsy. Pp. i-ix, 1-276, text-figs. 14/-. June 9, 1914. 


Odonata 


By F. C. FRASER 
Vol. 1. (Coenagriidae|. Pp. i-xiii, 1-423, map, text-figs. 25/-. March 1, 1933. 
Vol. Il. [Agriidae and Gomphidae]. Pp. i-xxiii, 1-398, 4 col. pls., text-figs. 25/-. 
Oct. 29, 1934. 
Vol. Ill. [Cordulegasteridae, Aeshnidae, Libellulidae]. Pp. i-xi, 1-461, map, 2 pls., 
text-figs. 30/-. Dec. 21, 1936. 
MOLLUSCA 


Vol. Ill. [Land Operculates (Cyclophoridae, Truncatellidae, Assimineidae, Helici- 
nidae)}|. By G. K. Gupe. Pp. i-xiv, 1-386, 2 pls., text-figs. 35/-. April 5, 1921. 
Vol. IV. [Freshwater Gastropoda and Pelecypoda]. By H. B. PRESTON. Pp. i-xi, 
1-244, text-figs. 14/-. March 31, 1915. 


WORMS 
Oligochaeta 
{In | Vol.] By J. STEPHENSON. Pp. i-xxiv, 1-518, text-figs. 30/-. June 30, 1923. 


Hirudinea 


{In 1 Vol.) By W. A. HarpinG [Rhynchobdellae| and . Percy Moore [Arhyncho- 
bdellae]. With an Historical Preface by the Editor, A. E. SHIPLEY. Pp. i-xxxii, 
1-302, 9 col. pls., map, text-figs. 25/-. March 23, 1927. 

Cestoda 
By T. SOUTHWELL 


Vol. I. [Cestodaria, Eucestoda (excl. Taenioidea)|. Pp. i-xxxi, 1-391, map, text-figs. 
May 29, 1930. 


22/6. 
Vol. II. [Taenioidea}. Pp. i-ix, 1-262, text-figs. 15/-. Dec. 29, 1930. 
Nematoda 


By H. A. BAYLIs 
Vol. 1. [Ascaroidea and Strongyloidea]. Pp. i-xxxvi, 1-408, map, text-figs. 25/-. 
March 23, 1936. 
Vol. II. [Filarioidea, Dioctophymoidea and Trichinelloidea]. Pp. i-xxviii, 1-274, map, 
text-figs. 17/6. Aug. 18, 1939. 
COELENTERATA, ETC. 
Freshwater Sponges, Hydroids and Polyzoa. By N. ANNANDALE. Pp. i-viii, 1-251, 
5 pls., text-figs. 14/-. Sept. 21, 1911. 
PROTOZOA 
By B. L. BHATIA 


Protozoa: Ciliophora. Pp. i-xxii, 1-493, 11 pls., map, text-figs. 30/-. Aug. 7, 1936. 
Protozoa: Sporozoa. Pp. i-xx, 1-497, 2 pls., map, text-figs. 30/-. Nov. 29, 1938. 
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MEETING OF 1953 MARCH 13 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Noel Elliott, B.Sc., 20 Eastholm, London, N.W.11 (proposed by H. Wildey); 
Joan Faunch, B.A., 6 Highmore Road, Blackheath, London, S.E.3 (proposed 
by P. J. D. Gething) ; 
*Josephine Mary Gilloch, B.Sc., Haulfryn, Park Street, Denbigh, North 
Wales (proposed by W. H. McCrea); 
Hans Q. Rasmusen, Voerslevgaard, Voerslev, Denmark (proposed by 
W. H. Steavenson); and 
Richard Quintin T'wiss, M.A., Sc.D., 53a Pall Mall, London, S.W.1 
(proposed by A. C. B. Lovell). 


The election by the Council of the following Junior Members was duly 
confirmed :— 


Paul Laurance Byard, 54 ‘Temple Fortune Lane, London, N.W.11 
(proposed by C. W. Allen); 

Peter Chadwick, Pembroke College, Cambridge (proposed by J. Hospers); 

Kenneth Hugh Fea, 4 Manor Drive, Hinchley Wood, Esher, Surrey 
(proposed by C. W. Allen); and 

Michael Godfrey Fisher, West Wood, Carleton, Pontefract, Yorkshire 
(proposed by C. W. Allen). 


Seventy-eight presents were announced as having been received since the 
last Meeting, including :— 

G. Abetti, Keplero (presented by the author) ; 

S. Chapman and T. G. Cowling, The mathematical theory of non-uniform 
gases, 2nd edition (presented by the authors) ; 

A. T. Doodson and H. D. Warburg, Admiralty manual of tides (presented 
by A. 'T. Doodson); and 

A. J. Thompson, Logarithmetica Britannica, being a standard table of 
logarithms to twenty decimal places. Part II (presented by the author). 

* Transferred from Junior Membership. 
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Meeting of 1953 April 10 
MEETING OF 1953 APRIL 10 


Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed : — 


Alan Frederick Collins, Xanadu, 85 Highridge Road, Bishopsworth, Bristol, 3 
(proposed by W. H. Steavenson) ; 

M. Zia-ud Din, M.A., Ph.D., Punjab University Observatory, Lahore, 
Pakistan (proposed by R. A. Newing); 

William Geddes, M.A., Stratford, Caldy, Cheshire (proposed by A. C. 
Sanderson) ; 

Frank Gwynne-Davies, M.A., 6 Bryn Road, Brynmill, Swansea, Glamorgan- 
shire (proposed by F. J. Sellers); 

James Nisbet, B.Sc., 72 Bridge of Weir Road, Linwood, near Paisley, 
Renfrewshire (proposed by W. M. Smart); 

*Brien George ‘Tunmore, B.Sc., Emmanuel College, Cambridge (proposed 

by D. E. Blackwell) ; 

Geoffrey Herbert Woodward, 19 The Steyne, Worthing, Sussex (proposed 
by T. A. D. Lawton); and 

Marshal Henry Wrubel, Ph.D., Indiana University, Bloomington, Indiana, 
U.S.A. (proposed by S. Chandrasekhar). 


The election by the Council of the following Junior Members was duly 
confirmed :— 


Alan Arthur Day, Trinity College, Cambridge (proposed by S. K. 
Runcorn); and 

William Douglas Hodgson, 28 Shaftesbury Street, Burnley, Lancashire 
(proposed by R. Stoneley). 


The President announced that the Council had elected the following 
Associates of the Society :— 


Viktor Amazaspovich Ambartsumian, Byurakan Observatory, Armenia, 
U.S.S.R.; and 
Albrecht Unséld, Institut fiir Theoretische Physik, Kiel, Germany. 


Eighty-two presents were announced as having been received since the last 
Meeting, including :— 

R. v. d. R. Woolley and D. W. N. Stibbs, The outer layers of a star (presented 
by R. v. d. R. Woolley) ; 

R. A. Lyttleton, The comets and their origin (presented by A.S.L.I.B.); 

H. Slouka (editor), Astronomy in Czechoslovakia from its early beginning to 
present times (presented by the Czechoslovak Astronomical Society) ; 
G. A. Shajn and B. F. Gaze, Atlas of diffuse gaseous nebulosities (presented by 

the Academy of Sciences of the U.S.S.R); 

L. J. Carter (editor), The artificial satellite: Proceedings of the Second 
International Congress on Astronautics (presented by the British Inter- 
planetary Society); and 

P. Jordan, Schwerkraft und Weltall (presented by the publishers). 


* Transferred from Junior Membership. 


THE VELOCITY DISTRIBUTION OF SPORADIC METEORS 


IV. EXTENSION TO MAGNITUDE +8, AND FINAL CONCLUSIONS 


Mary Almond, 7. G. Davies and A. C. B. Lovell 


(Communicated by A. C. B. Lovell) 


(Received 1953 March 3) 


Summary 


The final paper in the series on the velocity distribution of sporadic meteors 
determined by the radio echo technique presents the results for meteors of 
magnitude down to +8. The results are similar to those described previously, 
and give no evidence for a significant hyperbolic velocity component. ‘The 
scale of magnitude for the entire series has now been established, and the 
velocity distribution shows no change with meteor magnitude. Evidence 
is given that nearly 90 per cent of sporadic meteors must be moving in 
elliptical orbits with periods of about two years, and that the remainder 
have velocities close to the parabolic limit. It is impossible to reconcile 
these results with the theory that sporadic meteors are moving with hyperbolic 
velocities from an interstellar origin. 


1. Introduction.—In 1948 a series of radio echo measurements of the velocity 
of sporadic meteors was commenced in anticipation that a decisive answer could be 
given to the controversial question of their origin. In Part I* of this series of 
papers the results of the measurements during the autumn of 1948 and 1949, and 
the spring of 1950, were described. ‘The investigation was then extended to fainter 
meteors and the results for the autumn of 1950 and spring of 1951 have been des- 
cribed in PartII.+ The calculation of the theoretical velocity distribution for 
comparison with the experimental results was given in Part III.{ The present 
paper is the final one in the series. It describes the extension of the work to even 
fainter meteors during experiments made in the autumn of 1951, and gives a 
somewhat more detailed interpretation of the results than was possible in the earlier 
papers. In particular, the scale of meteor magnitudes, which was uncertain when 
Parts I and II were published, has now been established for the entire series. 

2. Apparatus.—The main purpose of the final experiment was to extend the 
autumn apex measurements of 1948, 1949 and 1950 to fainter meteors. It was not 
easily possible to increase further the transmitter power, or the size of the aerials. 
‘The necessary increase in sensitivity was therefore achieved by rebuilding the 
aerials described in Part II at an angle of 45° to the horizontal, the main axis being 
in an easterly direction as previously. Although the maximum gain of the aerials is 
decreased slightly, due to the loss of ground reflections, the mean range of the 
observed meteors is reduced from 400km to 150km. ‘The minimum electron 

*M. Almond, J. G. Davies and A. C. B. Lovell, M.N., 111, 585, 1951. 


+ M. Almond, J. G. Davies and A. C. B. Lovell, M.N., 112, 21, 1952. 
tJ. A. Clegg, M.N., 112, 399, 1952. 


29* 


412 Mary Almond, }. G. Davies and A. C. B. Lovell Vol. 113 


density detectable varies as R®*, and the net result is that the limiting magnitude 
for the new equipment was 1-5 magnitudes fainter than that used previously (see 
Section 4). 

The parameters of the apparatus were as follows :— 


Wave-length A =8-20 metres 
Peak transmitter power P=240kW 
Aerial beam | 
width to 
half power 
Aerial power gain over isotropic 
source G=28 
Receiver noise level €) =1°3 x 10 8 watts 


azimuth = +12° 
| elevation = +12° 


3. The 1951 apex experiment.—The basis of the experiment which was carried 
out between 1951 November 18 and December 6 was the same as in the 
previous apex experiments. The details of the observations, in which 
362 velocities were measured in 21 hours, are given in Table I. 


TABLE | 


Details of apex experiment; autumn 1951 
(Wave-length 8:20 m) 


Limits Total No. of Hourly No. of 
of watch time echoes rate velocities 
h 
22-08 15 
25-08 16 
08-08 12 
25-08 15 
19-06 28 
00-08 15 
12-08 I1 
20-08 
16-08 
14-08 
20-08 
19-08 
24-06 
12-06 
20-07 
27-08 


h 
° 
I 
2 
I 


NH 


Totals 21 03 2181 104 


‘The data for the individual velocity measurements are given in Table II, and 
the number distribution of velocities in Fig. 1. As in the previous experiment, the 
aerial had a subsidiary lobe radiating energy at a high elevation in an easterly 
direction. Meteors recorded in this lobe are distinguished by their short range, 
and 64 in this category are included in the measurements of Table II.* These 
are shown unshaded in Fig. 1. The smooth curve shows the theoretical distri- 
bution calculated as described in Part III, on the assumption that the meteors are 


* In this aerial system there was also a weak lobe at low elevation. Meteors occurring in this 
have ranges greater than 300 km, and one meteor of the 64 was excluded on this account. 


Date 
1951 
h 
Nov. 18 07 
19 06 
20 06 
21 06 
13 
08 
23 06 
25 06 
26 06 
27 06 
28 06 
29 06 
30 06 
06 
Dec. 
6 06 
362 
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moving at the parabolic velocity limit and in random directions. It is normalized 
to the 298 meteors recorded in the main aerial beam (that is, to the shaded areas of 
the histogram in Fig. 1). The errors of measurement have been combined with 
the theoretical distribution as described in Parts I and II. 


160 


Fic. 1.—The measured velocity distribution (histogram) for the 1951 apex experiment. The 
non-hatched areas show the influence of the antihelion component. The theoretical distribution, 
calculated on the assumption of uniformly distributed parabolic orbits, is shown in the smooth curve. 
The theoretical curve is normalized to the 298 meteors recorded in the main aerial beam (that is, 
to the hatched area). 

Ordinates: Number of meteors. 
Abscissae: Velocity in km/sec. 


4. The magnitude range.—Since the publication of Parts I and II, theoretical 


and experimental work on the scattering of radio waves from the ionized meteor 
trails has elucidated many of the problems which previously caused uncertainty in 
the interpretation of the received echo.* It is now clear that the equation origin- 
ally derived by Lovell and Clegg+, which was used to obtain the electron density 
from the echo amplitude in Parts I and II, is applicable to trails in which the 
electron line density is less than about 2 x 10!electrons/em. Except for the 
extreme tail of the histograms shown in Figs. 6, 7 and 8 of Part 1, the densities were 
below this value in the present series of experiments. Hence, renewed 
confidence can be placed in the data regarding the electron densities of Parts I 
and II. 

A more serious cause of uncertainty expressed in Parts I and II lay in the 
conversion of these line densities to visual magnitudes. Fortunately, this doubt 
has now been removed by Greenhow and Hawkins{, who have determined 
experimentally the relation between visual magnitude and electron density. From 
their results the following relation between zenithal magnitude M and electron 
line density « can be derived for values of x <10!"electrons/cm : 

M =35 — 2°5 logy). (1) 

‘The use of (1) to convert the electron densities of Parts I and II to visual 
magnitudes actually gives the same results as listed in Table VIIa of Part II for the 
1949 and 1950 experiments deduced from the estimates made before the work of 
Greenhow and Hawkins was complete. 


* 'T. Kaiser and R. L. Closs, Phil. Mag., 43, 1, 1952. J. 5S. Greenhow, Proc. Phys. Soc. B, 65, 


169, 1952. 
+ A. C. B. Lovell and J. A. Clegg, Proc. Phys. Soc., 60, 491, 1948. 
t J. S. Greenhow and G. S. Hawkins, Nature, 170, 355, 1952. 
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TaBLe II 
Data for velocity measurements during apex experiment; autumn 1951 


Electron line 
y=AR*/G density (a) 


(Signal to (R in units x 20"° 
cognate of 95 km) electrons/cm 
noise ratio) 

path 


Velocity Range 
VtdV R 
km/sec km 


Time 


3 
45°3+2°7 
53°4+311 
4°3 
355231 
62:°8+ 3:1 
47'6+23 
53°2+5°1 
56-7+3°3 
61:0+ 3:0 
48-7+2°7 
449+ 2°5 
59°4+2°0 
62:9+2°5 
54°7+3°2 
68-4+5°4 
57°6+43 
60°1+4°1 
40°3+2°9 
57°1+3°5 
62°8+4°9 
54°3437 
32°3+1°9 
67°34 3°4 
53°6+2°5 
57°9+43 
347443 
48°44+2°7 
43°212°2 
37°7423 
58:3+2-9 
96°27 2°7 
22°3+08 
34°742°3 
37'5+18 
45°8+2:-2 
65°7+4°0 
30°8+1°9 
3°3 
47°6+2°9 
42°94 2:2 
60-9 + 3°5 
60:0+2°8 
65:1+5°0 


414 
Date 
1951 
a I 2 4 5 6 7 
h m 
Nov. 18 07 25 100 II 
07 28 280 10 5°71 20°2 
07 37 170 12 1:78 6-3 
08 10 150 15 2°86 10°! 
19 606 28 100 14 
06 29 110 13 
baal 06 36 200 II 2°04 7:2 
06 41 120 13 6°84 24:2 
ae 06 43 180 9 1°40 49 
07 00 180 2°02 a1 
07 10 110 5 
07 14 190 13 2-19 77 
07 15 200 14 2°59 9:2 
07 19 150 12 2:29 8-1 
07 20 240 6 1°87 6-6 
07 26 190 13 2°17 
°7 39 200 14 2°59 9:2 
©7 44 180 9 1°40 4°9 
07 48 180 II 6-0 
©7 55 200 12 2°22 79 
07 56 240 13 4°06 14°4 
08 06 220 II 2°47 8-7 
08 07 160 14 2°19 
08 13 150 10 1-90 6:7 
20 606 30 200 12 79 
06 34 200 14 2°59 9:2 
06 45 140 13 3°47 12°3 
06 58 200 6 39 
07 03 200 13 2°41 8-5 
07 03 180 13 2:02 7°1 
07 14 180 10 1°55 5'5 
07 27 195 12 2°12 75 
07 29 180 12 1°86 6:6 
07 31 180 13 2°02 71 
07 34 170 13 1°93 6°8 
°7 35 180 14 217 77 
40 200 II 2°04 7°2 
43 160 13 2°03 
07 45 110 9 
21 06 34 160 13 2°03 72 
06 43 155 14 2°39 8-5 
06 51 120 II 5°79 20°5 
06 58 170 13 1°93 6:8 
07 05 110 II 
07 07 460 II 
19 120 13 6:84 24:2 
°7 27 180 13 2°02 71 
°7 34 245 10 3°23 11-4 
©7 37 240 13 406 14°4 
°7 49 200 13 2°41 8°5 
07 52 160 9 1°41 
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TABLE II (cont.) 


Electron line 
Velocity Range Amplitude density (a) 
Date Time A . : 10 
UT R (Signal (R in units x10 
95 aa km/sec km gn of 95 km) electrons/cm 


noise ratio) path 


3 


§2°342% 200 II 2°04 7°2 
08 03 64°6+2°7 170 14 2°07 73 
08 04 180 13 2°02 
08 07 160 13 2°03 72 
22 08 07 50°44 5°4 240 7 2°19 77 
08 10 34°8+2:1 160 13 2°03 
08 11 56-7+2°7 180 12 1°86 6-6 
23 12 68:7+ 4:2 210 13 2°60 
06 12 54°1+2°0 200 12 2:38 79 
06 12 62:5+2°9 245 2°26 8-0 
06 12 47°24 160 14 2°19 | 
06 13 32°4+1°4 170 II 1°63 5°8 
06 15 60:9+ 4°0 190 6 1°00 
06 19 65°6+3°7 200 4 0°74 26 
06 21 66-6443 160 15 2°34 8-3 
06 22 100 12 
06 23 100 13 
06 28 qo-7+5§'1 280 8 4°57 16:2 
06 33 20'°5+1°6 100 5 
06 39 3°5 180 9 1°40 4°9 
06 39 38:5+2°0 140 14 3°73 13°2 
06 40 68-3+5°2 240 9 2°81 10°0 
06 41 50'2+1°9 200 12 2:22 79 
06 50 69:6+4°4 190 12 2°00 71 
06 51 62°9+5°1 110 14 
06 51 70'2+46 180 14 2°17 77 
06 51 53°5+2°0 180 13 2°02 71 
06 56 100 12 
07 02 38-7424 190 9 1°50 5°3 
07 02 36-5+1°3 110 12 
07 05 64:6+2°0 180 15 2°33 8-2 
07 05 41:°6+3:°0 110 8 
07 06 39°'2+1°9 180 12 1°86 6-6 
07 07 53°6+3°2 160 7 1°09 
07 67°342°7 230 15 3°75 13°3 
07 12 63°9+1°8 190 14 2°33 8-3 
07 14 64:6+4'8 240 9 2°81 
07 24 59°244°7 140 6 1-60 5°7 
07 24 62°7+ 4°9 170 6 0°89 3 
07 25 56:2+1°8 160 15 2°34 8-3 
07 28 130 10 3°64 12°9 
07 29 58:5+43 190 6 1-00 
07 31 66-0+ 200 9 1°67 5°9 
07 32 66-:1+ 4°2 160 9 1°41 5°0 
07 36 28-7+1°7 110 12 
07 38 445-4416 180 8 1°24 4°4 
07 40 100 13 
07 43 190 6 1°00 3°5 
07 45 64:2+3°8 245 6 1°94 6-9 
47 58:14+2°3 180 II 1-71 6:0 
07 50 66:5 130 15 5°45 19°3 


46°3+1°4 


I 2 | 4 5 6 | 
h m 
07 50 170 15 2°22 79 
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TABLE II (cont.) 


Electron line 
Velocity Range Amplitude _ density (a) 
ate ‘Time V+aV A (Ri 
km/sec km : : of 95 km) electrons/cm 
noise ratio) 
path 
I 2 3 4 5 6 7 
h m 
Nov. 23. 07 52 §4°8+4:2 150 12 2°29 8-1 
07 54 42°6+1°9 160 13 2°03 7-2 
07 55 58:6+2°8 200 14 2°59 9:2 
07 55 180 10 1°55 5°5 
07 56 65:0+2°6 200 II 2°04 
08 03 68-9+ 4:2 210 13 2:60 9:2 
08 09 64°6+2°6 190 15 2°50 8-9 
25 06 27 54°5+2°7 150 6 40 
06 30 39°2+2°1 200 5 33 
06 30 65°5+4:°2 180 10 1°55 5°5 
06 34 63:°8+ 3:2 200 10 1°85 6-6 
06 36 s7°si27 180 13 2°02 
06 38 56-9+2°3 160 13 2°03 72 
06 45 56-4+471 190 8 1°33 4°7 
06 47 54°64+3:8 240 10 3°12 II‘! 
06 48 30°7+2°7 170 5 2°6 
06 48 33°6+1°9 170 6 089 3°1 
06 50 33°9+1°8 110 8 
06 53 22°7+1°9 160 9 1°41 5:0 
06 54 64:6+ 3°3 160 II 5-92 6-1 
07 00 56-5+2°7 180 10 1°55 
07 OI 36:84 2-0 7 
07 O1 36 150 6 1°14 
07 05 49°4+28 170 7 1°04 3°7 
07 08 39°'2+1°9 200 9 1°67 5°9 
07 11 74°0+4°9 150 14 2°67 9°4 
07 12 ro 8 
07 13 68:-7+ 5°4 200 14 2°59 9:2 
07 46°7+3°5 110 II 
07 29 50°2+5°5 190 9 1°50 5°3 
07 29 54°343°9 200 13 2°41 8°5 
07 33 38:34+2°7 110 12 
07 37 52°5+2°9 210 12 2°40 8°5 
07 40 64°5+4°9 190 10 1°67 
07 63°7+3°9 180 12 1°86 6-6 
O7 44 31'9+2°1 110 II 
45 24°5+1°0 100 5 
07 46 58:6+2-6 250 9 3°10 
07 50 68:4+4°4 160 II 
07 55 56-1+3°1 190 II 1°83 6°5 
08 00 70°9+3°4 250 13 4°48 15°9 
08 03 55°8+3'5 240 9 2°81 10°0 
08 04 56-4+2°9 110 12 
08 12 44°6+1°9 200 12 2:22 79 
26 06 20 3494271 150 7 1°33 4°7 
06 22 140 12 3°20 ¥i°3 
06 24 61°7+3°1 180 II 1-71 6-0 
06 33 34°84 1-4 130 12 4°36 15°4 
06 38 190 6 1°00 
06 39 120 10 5°26 18-6 
06 42 62:2+1°8 240 12 5-95 


06 46 46:8+1°6 200 13 8:5 
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TABLE IT (cont.) 


Electron line 
Date Time Velocity Range 4 y= density (a) 
V+dV R (R in units 
1951 U.T (Signal to 
km/sec km : : of 95 km) electrons/cm 
noise ratio) 
path 
1 3 5 6 7 
Nov. 26 06 56 62°34+2°2 170 14 2°07 73 
06 58 59'2+2°8 200 14 2°59 9:2 
07 00 75:0+2°8 180 15 2°33 8-2 
07 00 66°8+2°7 205 15 10°3 
07 02 34°8+1:-2 100 6 
07 04 110 6 
07 05 37°9+1°'8 170 II 1°63 5°8 
07 09 67°5+3°3 250 14 4°83 
07 43°2L1°4 100 14 
07 10 180 9 1-40 
07 ‘II 65:3+4:°2 160 15 2°34 8-3 
07 13 60°3+4°5 180 14 2°17 
07 21 62°5+3°9 160 10 1°56 5°5 
07 23 36-5+1°6 200 6 E52 3°9 
23 29°6+1°9 100 8 
07 23 s81+44 180 5 2°7 
07 24 16 110 9 
07 28 49'2+3°6 160 10 1°56 55 
07 29 60°5+4°6 190 8 1°33 47 
07 32 69:'1+3°8 140 13 3°47 12°3 
07 34 67:2+34 190 1°83 6°5 
©7 34 240 10 4:12 
07 35 56-2+3°3 190 8 1°33 4°7 
07 41 190 8 4°7 
°7 43 67°5+3°4 180 15 2°33 8-2 
07 46 56-4+ 4-2 200 10 1°85 6-6 
07 51 54°5+2:2 150 14 2°66 94 
07 52 57-324 180 15 8-2 
07 53 190 7 1°17 471 
07 54 54°64+23 110 12 
©7 55 190 5 0°83 2-9 
57 55°7+3'1 220 7 1°57 
07 58 63:2+4:°0 160 13 
07 58 59°4+44 200 8 1°48 
08 31°7+1°6 100 13 
08 03 54°3+6-2 200 15 2:78 9°8 
08 04 58-8+ 3°5 190 8 1°33 4°7 
08 10 46-2+2°1 150 9 1°71 
08 11 24°‘8+1°1 100 6 
08 13 58-7+2°8 200 II 2°04 7°2 
27. +06 18 33'8+2°0 150 5 0°95 3°4 
06 19 65°7+3°4 150 13 2°48 8-8 
06 25 32°9+1'8 100 14 
06 26 39'9+1°8 110 9 
06 33 54°1+3°2 160 9 1°41 5°0 
06 40 61:6+3°6 240 12 3°75 13°3 
06 45 5:0 140 10 2°67 
06 45 26-7+2°0 200 5 0°93 3°3 
06 50 54°8+471 170 7 1°04 
06 51 68:0+ 4°4 150 14 2°67 9°4 


06 57 58:1+2°8 180 13 2°02 
07 02 66-3+3°7 II 
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Velocity 
VtdV 
km/sec 


3 
az6+23 
67:1+3°3 
374416 
68-6+ 4°4 
47°5+3°2 
51°4+26 
57°5+28 
60°6+ 3:0 
66:5+2°8 
61-1+3°0 
67°1+5°3 
60-9+7°0 
34°74 2°3 
3°3 
40°6+2°9 
60°8+ 2°6 
20°44 1°5 
61:0+3°0 
60:2 + 3°7 
62:942°7 
33°4+2°2 
64°6+ 3:2 
2°71 
29°94 1°5 
39°4+2°0 
60°4+ 3°8 
50°4+4°4 
53:°7+4'8 
58:8+34 
39°5+1°9 
70'8+2°7 
58:2+3°3 
43°5+2°0 
61:7+4°7 
43°5+3°0 
45°7t2°4 
62:0+ 
60°5+ 2°5 
48-4+1°5 
29°21 + 
69:0+ 3°6 
28-7+1°4 
62°44 3:2 
54°8+3°1 


Range 
R 
km 


Amplitude 
A 


(Signal to 
noise ratio) 


y= AR*?/G 
(R in units 
of 95 km) 


Vol. 113 


Electron line 
density (a) 
electrons/cm 
path 


1951 U.T. 
2 6 7 
Nov. 27. 07 03 210 II 2°20 78 
07 12 240 9 2°81 
07 13 110 10 eee eee 
07 18 190 II 1°83 6°5 
07 24 200 13 2°41 
07 26 100 9 
07 28 160 13 2°03 72 
= 07 42 180 12 1°86 6-6 
©7 45 190 15 2°50 8-9 
°7 53 170 14 2°07 Es 
07 55 180 10 1°55 
07 56 160 II 1°72 61 
07 58 240 10 3°12 
08 o1 110 13 eee 
: 08 13 120 ¥3 6°84 24°2 
06 20 120 10 5:26 18-6 
06 20 150 13 2°48 8°8 
06 21 100 10 
06 21 160 10 1°56 5°5 
06 25 110 
06 25 170 14 2:07 73 
06 35 160 13 2°03 72 
06 36 170 10 1°48 5°2 
06 43 100 13 bee eee 
06 43 210 II 2°20 78 
06 45 110 
06 49 100 12 
06 51 150 13 2°48 8-8 
. 06 52 190 8 1°33 47 
07 00 160 10 1°56 5°5 
07 00 190 12 2°00 71 
07 00 160 II 6-1 
: 07 00 160 ‘12 1°87 66 
07 13 200 6 3°9 
07 22 170 12 1°78 6-3 
07 24 180 6 0°93 3°3 
07 24 200 15 2-78 9°8 
07 30 245 10 3°23 114 
07 38 110 14 
07 41 170 10 1°48 5'2 
07 44 160 12 1°87 6-6 
07 52 200 12 2°22 79 
°7 53 190 12 2°00 
°7 54 170 12 1°78 6°3 
07 56 180 13 2°02 71 
07 58 100 Il 
07 59 180 12 1°86 6-6 
08 02 190 7 41 
: 08 06 170 14 2°07 72 
08 10 160 12 1°87 66 
08 13 190 14 2°33 8-3 
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TABLE II (cont.) 


Electron line 
D Velocity Range y= AR*"*/G density (a) 
ate ‘Time R A (Ri 
km/sec km cule ented of 95 km) electrons/cm 
path 
I 2 3 4 5 6 7 
h m 
Nov. 29 06 21 36-44+1°4 160 10 1°56 5°5 
06 28 60°:1+2°4 200 10 1°85 6-6 
06 35 160 II 1°72 6-1 
06 37 60°5+2°5 180 12 1°86 6-6 
06 39 65:7+2°8 160 15 2°34 8-3 
06 41 110 14 
06 42 60:9+ 3:0 180 13 2°02 
06 43 60'1+ 46 170 13 1°93 6°8 
06 47 37°0+2°4 150 II 2°10 7°4 
06 50 12°32 110 
06 53 56-4+4°0 210 7 1°40 
06 54 696+34 245 9 2-90 10°3 
06 55 160 2°03 
06 55 38-8+1°7 110 II 
06 56 34°341°7 130 7 2°55 
06 59 64°6+ 4:1 150 10 1-90 6-7 
7 02 62:8+ 4:0 150 6 1°14 4°0 
07 04 67:°0+ 44 120 9 4°73 16°8 
07 140 9 2°40 8°5 
07 10 61°64+7°4 150 8 1°52 5°4 
07 10 55°9+3°2 200 13 2°41 8-5 
07 16 43°$+2°2 110 8 
07 23 33°6+1°0 160 13 2°03 
07 26 32°2t+1'°6 100 13 
07 30 58-3+6-7 210 9 1:80 6-4 
07 31 40°5+1°9 100 5 
©7 35 29'2+1°6 140 II 2°93 10°4 
07 36 65°34 3:2 170 9 1°33 4°7 
07 36 30°6+2°1 170 II 1°63 5°8 
37 36:1+2°0 210 5 1°00 3°5 
©7 37 66:1+ 8-0 160 12 1°87 6-6 
07 39 67°5+2°8 200 13 2‘41 8-5 
97 39 60°5+3°0 160 10 1°56 
44 60'9+ 3°1 160 13 2°03 
07 44 4°6 170 13 1°93 
©7 45 52°7+2°0 170 13 1°93 6°8 
©7 45 75°3+49 180 8 1°24 4°4 
07 48 65°4+4°0 210 12 2°40 
07 48 24°9+1°0 100 12 
©7 49 65°1+ 3°4 160 9 1-41 5° 
07 50° 34°5+16 160 6 0°94 3°3 
07 50 659+ 2°8 190 15 2°50 8-9 
©7 50 51°7+2°6 245 9 2°90 10°3 
07 52 69°9+5°8 120 15 7°89 27°9 
°7 53 160 10 1°56 5°5 
07:54 +31 245 15 4°84 
07 55 70°3+4°5 170 7 1-04 3°7 
07 57 44°6+2°4 190 
07 58 72°2+3°7 200 14 2°59 9:2 
07 58 140 15 4°00 14°2 


07 58 64°1+4°6 280 9 5°14 18-2 
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Dec. 
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TABLE II (cont.) 


Velocity 
V+dV 
km/sec 


3 


1°7 
62:1+7°5 
39°64 2°5 
73°7+46 
60°3+3°8 
64°8+3°3 
63°74 3°3 
53°6+41 
37°38+2°8 
56-:2+2°6 
34°34 
64°5+5°0 
33°52 
a7 
64°3+4°0 
22:2+0'8 
59°84 3°5 
24°7+1°3 
1°6 
32°74 
69°5+5°3 
67-1t4°2 
43 
51°8+1-6 
yorts8 
61°4+3°7 
43°44 1°5 
45°5+3°1 
51°5+2°9 
43°4424 
64°3+4°0 
45°7+4°5 
3-3 
37°5+23 
§6:0+2°7 
49°1t1°9 
66°3+3°2 
50°7+4°0 
41'64+1°4 
21:8+0°7 
54-7441 
65:6+2°8 
60°4+ 3°7 
33°22 
62-7+4°9 
56°7+2°3 
58-9+2°9 
32°3+1°9 


Range 


km 


Amplitude 
A 
(Signal to 
noise ratio) 


5 


y=AR®?/G 
(R in units 
of 95 km) 


Vol. 113 


Electron line 
density (a) 
x 20*° 
electrons/cm 
path 


7 
3°5 
6°7 
2°4 
15°5 
8-8 


6-6 


1951 Lime = 
I 2 4 6 
h m 
Nov. 29 07 58 190 6 1°00 
08 02 110 6 
08 04 150 10 1°90 
08 05 190 4 0°67 
08 07 240 14 4°37 
08 08 150 9 I°7! 
08 09 180 13 2°02 
08 10 140 13 3°47 
m 06 26 150 13 2°48 
06 27 190 13 2°17 
06 27 110 6 
06 28 200 10 1°85 || 
06 29 110 14 
; 06 41 110 14 
: 06 44 200 Il 2°04 7°2 
06 45 160 13 2°03 72 
06 46 100 14 on Gs 
06 48 180 ? >2°69 >9'5 
- 07 20 210 8 1°60 5°7 
07 22 100 6 was 
6 06 42 170 13 1°93 6°5 
06 45 120 13 6°84 24°2 
06 49 245 10 3°23 11-4 
06 50 180 12 1°86 6-6 
06 51 190 9 1°50 
07 05 170 15 2°22 7°9 
07 06 140 10 2°67 9°4 
07 08 190 12 2°00 a 
07 12 170 14 2°07 a3 
15 190 12 2°00 
07 180 9 1°40 4°9 
07 15 150 10 6-7 
07 20 180 8 1°24 4°4 
07 26 240 II 3°44 12°2 
07 26 100 12 ose oe 
07 26 160 13 2°03 72 
: 07 29 180 14 2°17 74 
07 37 160 13 2°03 72 
07 38 245 13 4°19 14°58 
07 38 120 8 4°21 14°9 
07 40 100 16 
©7 45 170 5 0-74 26 
©7 45 170 13 1°93 6-8 
07 57 160 15 2°34 8-3 
07 58 170 14 2°07 7°3 
07 58 100 10 
08 o1 160 10 1°56 5°5 
08 04 180 15 2°33 8-2 
o8 14 170 II 1°63 
08 15 150 8 1°52 5°4 
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(a) The magnitudes of the 1951 apex meteors.—The process of reduction of the 
received echo amplitude to electron line density «, for the measurements described 
in this paper, was the same as in Part II. A quantity y = AR*2/G is calculated for 
each meteor. A is the amplitude of the radio echo from the first Fresnel zone 
compared with the receiver noise level, R the range of the meteor and G the 
effective gain of the aerial system. G depends on the range R and on the azimuth 
of the meteor with respect to the main beam of the aerial system, and the method 
of computing y has been described in the Appendix to Part I. The appropriate 


100 200 300 


Fic. 2.—Relation between most probable sensitivity § and range R, for the aerial system 
used in the 1951 apex experiment. From these curves the values of y are calculated for each meteor. 


Ordinates: Most probable sensitivity s. 
Abscissae: Range in km. 


20 


0 4 8 12 16 20 24 28 


Fic. 3.— Number distribution of electron line density (x) for the 1951 apex velocity measurements, 
The distribution is for the 298 meteors recorded in the main aerial beam (corresponding to the 
hatched distribution in Fig. 1), but one meteor for which the only information is x >-9°5 x 101° is 
not included. 


Ordinates : Number of meteors. 
Abscissae : Line density « electrons/cm path ( x 10"). 


sensitivity-range curve for the aerial used in the present experiment is given in 
Fig.2. The line density « for each meteor is then calculated from the Lovell—Clegg 
formula 
me? € 
a= Pe’ (2) 
where the values of «,, P and A are those given in Section 2. For each individual 
meteor the values of y and « are given in Table II, columns 6 and 7, respectively. 
The corresponding number distribution against electron density is shown in 
Fig. 3 for the 298 meteors occurring in the main beam of the aerial. ‘These values 
of « can be converted to zenithal magnitude M by using (1) and it is found that the 
peak of the curve in Fig. 3 corresponds to magnitude + 8-0. 
In Part II an alternative method of estimating the equivalent visual magnitude 
was described. This involves a comparison of the hourly rates for a visual 
observer and for the radio echo apparatus. In the case of the present experiments, 
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the procedure gives a value for the zenithal magnitude at the peak of the distri- 
bution curve of +8-5. 

(6) Comparison of electron density and magnitude distribution for all the experi- 
ments.—The difference in the estimates of meteor magnitude by the electron 
density method and the visual rate method amounts to only 0-5 magnitude. Thisis 
comparable with the uncertainty in the magnitudes derived by either method, and 
since the visual method cannot be used for the antapex experiments it is convenient 
in comparing the distributions for the entire series to use the magnitudes calculated 
from the line density by equation (1). The relevant information from Parts I and 
II is combined with the results of the present paper in Table III.* 


TaBLe III 


Distribution of electron densities and zenithal magnitudes 
in the apex and antapex experiments 


Number of meteors 


Line density 
a electrons/cm 
path ro!” 


398-631 

251-398 

159-251 

100-159 

63°1-100 

39°8-63°1 

25°1-39'8 

15°9-25'1 

10°O-15'9 
6°31-10°0 wae 
3°98-6°31 
2°51-3°98 eee eee 
1°59-2'51 


Zenithal 
magnitude 
M 


1948 1949 1950 1950 I95I 1951 Totals 
apex apex antapex apex antapex apex antapex apex 


5. Discussion. (a) Variation of velocity with magnitude.—In Parts I and II the 
various velocity distributions were examined for any change of velocity with meteor 
magnitude, with negative results. During the four apex experiments the velocities 
of 951 meteors were measured. 84 of these occurred in side lobes of the aerial 
system, and the electron density of a further 12 was uncertain. ‘The remaining 
855 for which electron densities are known have now been divided into magnitude 
groups as in Table III, and the velocity distribution plotted for each group. The 
results are shown in Figs. 4 (a to f) from which it is clear that no significant change 
occurs in the velocity distribution over a range of five magnitudes. 

(b) The velocity distribution observed in the apex experiments.—The velocity 
distribution of Fig. 1 is similar to those of Parts I and II and shows that in the case 
of these fainter meteors there is still no indication that any appreciable proportion 
of the meteors are moving with velocities in excess of the parabolic limit. There is, 
however, a marked discrepancy between the observed velocity distribution and 
that calculated on the assumption that the meteors are moving in random directions 

* Meteors occurring at short range in the side lobes of the aerial system (non-hatched areas in 


the velocity distributions) are excluded from Table III, as also are those meteors whose electron 
densities are not accurately known (prefixed with >, in the appropriate tables). 


ie. 38s +4 

+ 

: 7 15 15 238 +6 
13 31 
157 
2 67 
a5 4 93 +8 
I I +9 
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24 


24 


op 
30 


10 30 50 70 90 10 30 50 70 


(a) Distribution for all magnitudes (855 meteors). (b) Distribution for magnitudes +3 and 
+4 [35 meteors (magnitude +4) with « 
between 100 and 25110! electrons/cm ; 
1 meteor (magnitude + 3) with « between 251 
and 631 electrons/cm). 


24 


0 0 t+ +++ 
10 30 50 70 90 10 30 50 70 90 
(c) Distribution for magnitudes +5 


(170 meteors with « between 39°8 and 
100 10!” electrons/cm). 


(d) Distribution for magnitudes +6 
(238 meteors with x between 15-9 and 
39°8 10! electrons/cm). 


24 


24 


10 30 50 70 90 10 30 50 


(e) Distribution for magnitudes +7 
(317 meteors with x between 6°31 and 
15°9 10!° electrons/cm). 


(f) Distribution for magnitudes +8 
(93 meteors with « between 2°51 and 
6-31 x 10!" electronsicm and 1 meteor 
with a between 1°59 and 2°51 x10! 
electrons/cm). 


Fic. 4.—The velocity distribution for different magnitude groups observed in the apex 
experiments of 1948, 1949, 1950 and 1951. A _ total of 951 velocities was measured in these 
experiments; but 96 are excluded from these distributions either because they occurred in a side 
lobe of the aerial system or because the value of electron density was uncertain for other reasons. 


Ordinates : Percentage of total in group. 
Abscissae : Velocity in km/sec. 
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at the parabolic limiting velocity. ‘The most probable velocity in the theoretical 
distribution is near the parabolic limit, that is 72 km/sec, whereas the experimental 
curve peaks at about 60 km/sec. A similar discrepancy was also found in Part Il 
for the 1950 apex measurements, and in Part I for the 1948 and 1949 apex 
measurements. * 

This discrepancy could, of course, arise from instrumental causes, or from 
deceleration of the meteors. The question of instrumental limitation was dis- 
cussed in Part I, and with the 8 m apparatus it was shown that no instrumental cut-off 
existed, at least for velocities less than 140km/sec. Any appreciable deceleration 
would influence the measurements on the shower meteors, and the absence of this 
effect is clearly demonstrated by the results of measurements with the same 
apparatus on the major meteor showers. In particular Hawkins and Almond+ 
have summarized the results for the Perseid shower, where this apparatus gave 
the velocity as 60-5 km/sec, the computed orbit being in close accord with that of 
Comet 1862 III, with which the meteor shower is known to be associated. It can 
be concluded, therefore, that the discrepancy between the theoretical and 
experimental curves is a real effect, arising because of some error in the two basic 
assumptions of the theoretical treatment. ‘These will be examined in turn. 

(i) The calculation is based on the assumption that the meteor paths are distri- 
buted at random. If this is not the case then there is clearly a possibility that the 
velocity distribution found experimentally might be explicable in terms of a 
particular asymmetrical distribution with the heliocentric velocities remaining at 
the parabolic limit. Consideration of this case indicates that the observed velocity 
distribution could be produced if there were a very large concentration of radiants 
from an area near the zenith, about 45° from the apex and in ecliptic latitude N. 45°. 
However, such a concentration has never been observed, either in the visual} or 
radio echo § work on the sporadic meteor radiants. 

(ii) The calculation is also based on the assumption that the heliocentric 
velocities of the meteors are at the parabolic limit. In order to explain the dis- 
crepancy there seems to be no alternative to the conclusion that the velocities must 
be considerably less than parabolic. Simple considerations indicate that under 
the conditions of the apex experiments, a random distribution of meteors with 
heliocentric velocities of about 34 km/sec would produce the observed distributions. 
The bulk of the sporadic meteors would therefore appear to be moving in elliptical 
orbits with periods of about two years—a conclusion which could hardly be in 
greater conflict with the hyperbolic orbits required to explain the results 
of Hoffmeister|! and Opik.{ 

(c) The velocity distribution observed in the antapex experiments.—The antapex 
experiments were designed primarily to demonstrate the fact that no instrumental 
cut-off was influencing the results. However, the distributions show a similar 
discrepancy between theory and experiment as in the case of the apex experiments. 
The broader shape of the theoretical distribution makes it more difficult to analyse 

* As mentioned in Part II, the theoretical curves in Part I are incorrect due to an error in the 
summation for the theoretical distribution. The correct theoretical curves are similar to those of 
Part II and the present paper, with the most probable velocity at nearly 72 km/sec. 

+ G. S. Hawkins and M. Almond, M.N., 112, 219, 1952. 

tJ. P. M. Prentice (to be published). 

§ Jodrell Bank (to be published). 


|| C. Hoffmeister, Die Meteore, Leipzig, 1937. 
4 E. J. Opik, Tartu Obs. Publ., 30, No. 5, 1940; ibid., 30, No. 6, 1941. 


No. 4, 1953 The velocity distribution of sporadic meteors 425 


the cause of the discrepancy, but the following considerations indicate that the 
results are in general agreement with the conclusions reached above. The visual* 
and radio echot surveys of the sporadic radiants have shown that there is a marked 
concentration of direct orbits near the ecliptic. If the radiants produced by these 
orbits were included in the strip of sky to which the apparatus is sensitive, a pre- 
ponderance of radiants near the antapex would be observed, and these would 
produce a peak in the velocity distribution at 17 km/sec for parabolic orbits. As 
few meteors are observed with velocities less than 20 km/sec it is evident that these 
meteors do not enter the distribution; their radiants lie at points on the ecliptic at 
some distance from the antapex, outside the strip observed in these experiments. 
As each antapex experiment was performed with the same apparatus as the previous 
apex experiment it is possible to calculate the total number of meteors that would be 
expected in the antapex experiments from the rate in the corresponding apex 
experiment on the assumption of random distribution, taking into account the 
concentration of radiants towards the apex and the relative times of watch. In the 
case of the 1950 antapex experiment the number expected is 38; 36 meteors were 
actually observed in the main lobe of the aerial. In the case of the 1951 antapex 
experiment the expected number is 21, compared with 53 actually observed. This 
would indicate that the distribution of radiants along the strip is not far from random 
and again leaves low heliocentric velocities as the only plausible explanation of the 
discrepancy between the theoretical and experimental curves. t 

(d) Meteors seen in side lobes of the aerial system.—Each of the aerial systems 
used had a subsidiary lobe radiating energy at high elevation in an easterly 
direction. Whilst the power radiated in this lobe was relatively small, meteors 
were seen in it at short range, and hence the sensitivity was relatively high. ‘These 
meteors can be distinguished by their short ranges, and have been separated in 
Fig. 1 of this paper and in Figs. 1 and 2 of Part II by leaving them unshaded.§_ The 
most probable radiant point for these meteors is low in the west, that is from near 
the antihelion point in the apex experiment, and from near the helion point in the 
antapex experiment. ‘These are the regions where the radiant survey indicates 
marked concentrations of activity. Since the exact elongation from the apex of 
the radiant points is uncertain it is not possible to compute the heliocentric 
velocities, but estimates of the probable elongation show that the meteors are also 
non-hyperbolic. All theoretical curves have been normalized to distributions not 
containing these meteors, since they do not come from the strip of sky considered in 
deriving the theoretical curves. Whilst the proportion of these meteors is much 
higher in the antapex than in the apex experiments, it should be noted that the 
number per hour of watch is lower in the antapex experiments. 

(e) The high-velocity tail.—Of the 867 meteors measured in the main lobes of 
the aerials in the apex experiments, 7 yielded velocities in excess of 80 km/sec. 
While this is less than 1 per cent of the total, it is important to examine whether 
they can be attributed to errors of measurements, or do, in fact, represent a hyper- 
bolic velocity component. Full details of the measurements on these meteors are 


* J. P. M. Prentice, ibid. 
+t Jodrell Bank, sbid. 
t As the strips lie far from the ecliptic except at the apex and antapex, the result is not in 
conflict with the ecliptical concentrations found in the radiant survey. 
§ Although not indicated in the diagrams a similar distinction applies in the cases of Figs. 1, 2 
and 3 of Part I. 
30 
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given in Table IV. Theerrors quoted are calculated from the number of pulses in 
each zone measurement as described in PartI. Alternative methods of computing 
the errors have since shown that the one adopted gives a rather high value for the 
error. However, it allows a direct comparison to be made with the errors obtained 
by the same method in other velocity ranges. 


TaBLe IV 
Details of velocities in excess of 80 km/sec 


Electron 
density 
10!/em 
path 


29 
30 
136 
119 
42 
26°8 
4°9 


Amplitude 
x receiver 
noise level 


Velocity Range 
Date 
1948 Nov. 80-5+ 2°9 275 
1949 Nov. 83:3+ 6-9 380 
1949 Nov. 82-4+ 5:0 700 
1949 Nov. 96-9+ 12°6 615 
1949 Dec. 87-0+ 96 460 
1950 Nov. 05 81 700 
1951 Nov. 07 83°5+11°1 180 


Meteors with velocities between 65 and 75 km/sec have been selected and the root- 
mean-square value of the errors quoted in Table II, and in the appropriate Tables 
of Parts I and II, have been calculated. The results are given in Table V. 


TABLE V 


Mean errors at 70 km/sec 


Experiment 1948 1949 1950 1951 Total 
No. of meteors II 62 152 
R.m.s. error km/sec 4°3 4°6 


The r.m.s. error of the 7 meteors of Table IV is 8-6 km/sec, nearly twice the mean 
value for meteors in the range 65-75km/sec. It is therefore clear that these 
meteors have unusually large errors, and their true velocities may well lie within 
the parabolic limit. 

The number is, however, considerably more than would be expected, even if 
the errors of measurement used in deriving the theoretical curves were doubled, 
and the observed velocities did conform to the theoretical curves. As only some 
10 per cent of the total appear to have velocities close to the parabolic limit, the 
amplitude of the theoretical curve should be reduced to one-tenth its present size 
and a further curve superimposed with a peak at a lower velocity. The effect of 
this would be to lower still further the amplitude of the tail of the distribution 
above 80km/sec. The number of velocities above 80km/sec that might be 
expected on statistical grounds from the errors, and assuming parabolic velocities, 
is 0-15. If the errors are doubled, as indicated from the above discussion, the 
number becomes 1-5. However, the evident lack of parabolic velocities in the 
distribution will reduce this probable number to abouto-2. Itis therefore improb- 
able that the high velocities of all the 7 meteors are due to errors, and while 
the total number is too small to reach any definite conclusion it is likely that 
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somewhere between 1 and } per cent of the 867 meteors observed in the apex 
experiments had velocities in excess of the parabolic limit. 

6. Conclusion.—In the work described in the present series of four papers, 
the velocities of 1 095 sporadic meteors have been measured. ‘The measurements 
cover a range of visual magnitudes from +3 to +8. ‘There is no indication that 
any appreciable number of meteors have velocities in excess of the parabolic limit ; 
on the contrary the results can be explained only by assuming that nearly 90 per 
cent of sporadic meteors are moving in elliptical orbits with periods of about 2 years. 
About Io per cent have velocities close to the parabolic limit. The results do not 
exclude the possibility that less than I per cent may have hyperbolic velocities—a 
result which may imply an interstellar origin or may equally well be explained on 
the basis of perturbations within the solar system. ‘The results of these measure- 
ments are in good agreement with the radio echo measurements made by 
McKinley.* 

The work was initiated in an attempt to settle the controversy over the origin of 
the sporadic meteor component and the result is entirely in favour of the belief that 
the sporadic meteors are essentially permanent members of the solar system. The 
widespread belief in the hyperbolic theory has arisen mainly from the work of 
Hoffmeister and Opik. Hoffmeister’s conclusions were based on studies of the 
diurnal and annual variation of meteor numbers, and in the light of recent 
surveys of the distribution of sporadic radiants by Prentice, and by radio echo 
measurements at Jodrell Bank, he has shown how his observations can be 
reconciled with these recent results.+ Opik’s results were based on velocity 
measurements with the rocking mirror during the Arizona expedition and later in 
Tartu. It is cleariy impossible to reconcile these results with the material 
presented in this series of papers. 
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PHOTOGRAPHIC PLACES BY A POSITION-LINE METHOD 
G. Merton 
(Communicated by the Director of the University Observatory, Oxford) 
(Received 1952 December 31) 


Summary 


A method is given for deriving, from an approximate position of the image 
of an object on a photographic plate and the accurately measured distances 
between its image and those of the comparison stars, the corrections necessary 
to give the accurate equatorial (or other) coordinates of the object. A simple 
graphical method, akin to the navigator’s “ position-line’’ method, is 
described for obtaining these corrections. Approximations and tables, to 
facilitate the computations, are added in appendices. 

The solution is applicable to the spherically curved plates of such types 
of reflecting telescopes as Schmidt’s, as well as to the flat ones of refractors 
(even those of large field); and the method obviates the need for special means 
to measure the plates in rectangular coordinates. 


Introduction.—The methods at present in use for deducing the equatorial 
coordinates of a celestial object recorded on a photographic plate, from nearby 
stars of known positions, require the measures to be made in rectangular 
coordinates. This necessitates the use of some form of micrometer specially 


arranged for making the measures in two coordinates at right angles; recent 
experience has shown that an instrument of this kind is not always readily 
available, especially where an individual worker is concerned or where very 
large plates are in use. 

There is thus need for a method which does not employ rectangular 
coordinates and for which an ordinary short screw micrometer, or some form 
of fine scale, will suffice to make the necessary measures. Such a method is 
described in the present paper. Its essential feature is that accurate measurement 
is confined to that of the linear distances of the object from each of the images 
of three or more adjacent comparison stars. Then, just as in the position-line 
method of navigation, these accurate measures, together with some approximate 
values of other quantities required, enable the observer to locate the object 
within a ‘“‘cocked hat” with an accuracy which depends essentially on the scale 
of the photograph and the accuracy with which the equatorial (or other) 
coordinates of the comparison stars are known. 

Because of radial symmetry, the method is applicable to photographs taken 
at focal surfaces of any radius of curvature, and in particular may be usefully 
applied to those obtained with Schmidt and Maksutov telescopes. Moreover 
it can be used at large off-axis distances where the dependences method, and 
most other two- and three-star reduction methods, fail in accuracy.* Further, 
should images of the object fall at the edges of two overlapping plates, measured 
distances of it from stars on one plate can be used in conjunction with those from 
other stars on the other plate, since in this method each measured distance is 
independent of the others. 

* See final remarks in Appendix C (p. 449). 
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It is possible that the procedure proposed may show some small gain in 
accuracy over conventional methods for three reasons: single direct measures 
of distance are made involving, that is, one coordinate instead of two; the 
measures include the object, whose position is to be found, each time and not 
only once in a set of measures; and lastly, these measures are made in various 
position angles and this should help to decrease the effects of setting errors which 
arise from irregularities in the image of the object concerned. 


Part 1. Theory of the method 


1.1. Angular distances from measures.—We proceed first to derive the 
equations needed to obtain the great circle angle 6 on the celestial sphere between 
two objects, the distance s between the images of which on a spherically curved 
photographic film or plate, at S, and S,, has been measured. The equations 
for a flat plate will then result as a special case of the general solution. 


Fic. 1. 


Let O be the centre of the object glass or mirror of the camera, OC the 
optical axis with C the plate centre and f=OC the focal length (Figs. 1 and 2). 
The photographic recording surface is assumed to be an area on the surface of 
a sphere of radius 7, with its centre R on the optical axis. The measured 
distance, s, between S, and S, is assumed for the present to be the true measure 
of the great circle arc between these points. 

Then 


s= (1) 
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where Aj, is the angle subtended, at the centre of the sphere, between 5, and Sg. 
Denoting the chord 5,5, between them by s,, 


=arsin } (2) 


In the case of a flat plate r-+co and the measured distance s replaces s, in 
the equations that follow. 

Applying the cosine formula to the plane triangle S,OS,, where c,=S,O 
and ¢,=$,0, 

= + — cos 8, (3) 

from which 
5. 

— (4) 

The distance c in each case is found from the plane triangle SRO, where A 
is the angle subtended at the centre of the sphere between S and C, and RO is 
a straight line. If we treat r and f as always positive then we must put 


RO=r¢+f, (5) 
the positive sign being for cases where the photographic surface is convex to 
the incident light (Fig. 1), and the negative where it is concave (Fig. 2). Thus 

(6) 

Since the distance g=CS is measured in each case from the plate centre 
(optical axis) along an arc on the sphere, we have 


A=g/r. (8) 


Introducing an instrumental constant «, defined by 


e=1tfi/r, (9) 
and an auxiliary angle ¢, defined by * 


sin? }6 = 


tan = 2k = sin =, (10) , 


f 


c=fsec 9, (11) 
and (4) then gives the equation required: 


2 
sin? =} cos COS E — (sec — sec (12) 


This equation may be modified, to simplify the computational work, by 
introducing approximations and tables based on them. The discussion of this 
aspect however is postponed to Appendix A since the exact equations (12), 
together with (2), (9) and (10), may be used as they stand to find 6. 

A case not so far considered is where the centre of curvature R lies between 
the objective and the plate. x? is then negative, cosA= —g/r; and the solution 
has to be modified by defining the auxiliary angle y by 
f , 
sin p=2 =sin=. (10a) 


* See Appendix A (A.5, p. 444), and Table I in Appendix B. 


we can rewrite (7) simply 


S,=2rsin 
dy 
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Then 
c=fcos 9, (11a) 
and 


2 
sin? = } sec sec E — (cos — cos | ; (x? negative). (12a) 


In the case of the flat plate of a photographic refractor r—-co, x=1 and 
(10) becomes 


tan p=g/f, (x =1), (13) 


where g is the base of a right-angled triangle, and is the angle SOC in the 
figures; and equation (12) is immediately deducible from the plane triangle 
S,OS8,, bearing equation (11) in mind. 

1.2. Accuracy of measurement of g.—It is necessary at this point to consider 
the accuracy to which the values of ¢ are needed, and to show that the method, 
at any rate when applied to the cases normally arising at present, is not 
complicated by any need to make additional accurate measurements, namely of 
the distances g from the plate centre. 

The maximum value of ¢ is given by an image at the edge of the plate and, 
in the case of the flat plate, is half the angular field of the plate. In the case of 
the ordinary type Schmidt telescope « = \/2 and g given by (10) is larger than 
for a flat field by about this factor; but, on the other hand, the half field in 
present-day instruments of this type does not exceed 34°. In practice therefore p 
is unlikely to exceed 6° and will generally be much less, since objects are not 
usually measured at the edges of plates. 

Taking (12) we may rewrite it 

2 


4sin? 36 = COS COS Po + 2— (SEC COS COS |, SEC Ho), (14) 


and also 
tan? + 3( — (15) 
1 —tan* + pa) tan? — 


2 
Differentiating (14) 
As, 
4A8@ . sin cos COS COS Po 


COS? — COS? Po 
— Ag, tan (53 COS COS Po 608 ) 


33 cos? — cos? 
—Ag,tan (55 COS COS + (16) 


Neglecting the second-order terms of (¢,— 2) in (15) we use the equation, in 
the simplified form omitting the last term, to divide the left-hand side of (16) 
by 4sin? 36 and the right-hand side by (cos g, cos )s*/f?. Also, since s, is small, 
we may replace it by the measured arc s, so that (16) can then be written 


> 2 —Agy,tan 9, E + (tan? — tan? 


—Ag,tan E — (tan? — tan? | (17) 
Next, to the first order in g/r, (10) may be written 
tang=«gif, P2)- (18) 


| 
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Differentiating this and multiplying by tan y cos? » 


Ag. tang= sin* (19) 


which can be substituted outside the square brackets of (17), while (18) is 
substituted inside them. Thus, finally, 


Aé = tan 30 int 1 {1420 (8 g 


where 7 = }(g,;+g2); and the equation for the mean error of @ is 


(\6)? tan? 46 [4 (=) + sin! {; + 2K? (8 


In assessing the permissible error in Ag we may insert the maximum 
numerical value of (g,—ge)/s, namely unity, which occurs when the two images 
lie on the line joining them to the plate centre. Then consider the Ag term 


rewritten in the form 
(2 +2xt2) (21) 
8g g 


When g, and therefore y=xg/f, is small the error in g can be very large, 
compared to the error As in the measurement of s, because the factor sin? p 
is very small. As g increases s/g->o and g/g—>1; consequently the error in g 
can then be of the order of (cosec? g)/x? that in s to give the same error in 0. 
Thus when g=6°, for an image at the edge of a Schmidt camera field of 84° 
diameter, g can have an error 46 times that in s to give the same error in 0. 
As will be shown later (Appendix A.6, p. 444) the measurement of g with a 
millimetre ruler will usually suffice. 

1.3. Approximate position and its correction.—For the next step we shall 
require an approximate position (%», 59) for the object, in order to compute its 
angular distances (p;, p2, p3) from the comparison stars, and compare them with 
those (4,, 92, 8,) obtained by (12) from the measures (s,, 59, 53). In the subsequent 
correction of this approximation a diagram will be needed, for the graphical 
method of solution, having lines drawn from the position of the object, as origin, 
in the directions of the comparison stars. 

In order to ensure accuracy in the diagram, and sufficiently correct orientation 
of each direction line, the best procedure will probably be to plot the tabular 
positions of the stars on as large a scale as may be convenient.* If, then, from 
each comparison star an arc be drawn with radius 0, so that all three arcs 
intersect, the intersections will indicate the position of the object, and its 
coordinates taken from the diagram will then yield its approximate right 
ascension and declination. Lastly, lines are drawn joining the star positions to 
the point representing the object, and they will be sufficiently correct in orientation 
for the diagram to be used for dealing with subsequent corrections. 


* In the case of a large-scale plate it may be feasible to use the plate itself, laid on graph paper, 
to give the direction lines in fairly correct orientation, and to lay off the measured distances @ along 
them, from the object as origin, and mark the star positions. 


ae 

1 A 

— sin? (20) 
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‘The angular distance (p), of each comparison star from the approximate 
position, may now be computed as follows. In the spherical triangle formed 
by the comparison star, object, and pole, we have 


COS p; = Sin sin 8; + Cos cos 5, COS — (22) 
and thus 
sin® }p, = sin® — 8,) + cos cos 8, sin? — %,), (23) 


and similarly for pg and ps. An approximation and a table for simplifying the 
calculation of these distances are given in Appendices A and B at the end. 
1.4. Next, the differences 


D,=0,+A0,;—p,, and similarly D,, Ds, (24) 


are found which represent the errors in the assumed position of the object. 
The quantities A@,, etc. have been introduced here to take into account the 
small corrections to be made for the effects of factors which have still to be 
discussed—generally only the effects of differential refraction need to be 
included. 

The correct coordinates («, 5) for the object could now be obtained from 
the residuals by differential correction. It is however proposed to use a 
graphical method for obtaining the corrections to the assumed position, since, 
with the diagram already available, it will be more convenient and more 
informative. 

1.5. Position-line method of correcting assumed place.—The method of 
correction now proposed is closely akin to that used by the navigator to fix his 
position from altitude observations of stars and a dead-reckoning (D.R.) position 
of his ship; the equivalent items will be indicated in brackets. Relative to each 
of the comparison stars the position of the object (ship) will lie on the circum- 
ference of a circle with radius equal to 6 (angular distance from sub-stellar 
point = zenith distance). The difference (intercept) between @ and p, computed 
from the assumed (D.R.) position, being in each case small, only a short piece of 
the arc of each position circle need be drawn on the diagram (Mercator chart) 
in the neighbourhood of the intersections with the others. These arcs can 
usually be represented by straight lines at right angles to the radii, and the 
object (ship) is assumed to lie on these lines (position lines) and to be at the 
centre of the triangle of intersections (the “‘ cocked hat’’). 

Each difference, or intercept, given by (24) is therefore marked off from 
the object, as origin in the diagram, along its appropriate direction line, positive 
away from the comparison star, but on a scale much larger than previously used, 
say I unit=o’-o1; and a position line is drawn at right angles through the 
point so indicated, so that the three position lines intersect to form a triangle. 
Fig. 3, which refers to the example in Appendix C, will illustrate this procedure. 

The position triangle may be expected to be small in size, probably less 
than o’-1, and the point inside it which is equidistant from the position lines 
(the incentre) may be taken as indicating the true position of the object relative 
to that assumed. From it the corrections (Ax’, AS’) may be read and applied 
to the assumed «, and 5), namely 


= 4 Ax’. sec 8, 
corrected x=%+Aa®; 85=5,+ Ad’. 


(25) 
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Even if the corrections are only of the order of 0’-2 it may still be desirable to 
repeat (23) and (24) not only as a check but because the position lines are really 
ares of circles and, unless the errors plotted are small enough compared with 
the distances 6,, 0, and 6; to permit the position lines to be drawn as straight 
lines, the corrections will not be found exactly enough. Also the need to correct 
cos 8, in (23) must not be overlooked. 


- 14°10'- 


Fic. 3. 


Part 2. Factors affecting measured distances 


2.1. The position triangle must now be considered more closely. It arises 
not only from the errors of the catalogue star positions used in (23) but also from 
factors now to be detailed which affect the positions of the images and thus the 
angular distances deduced by (12). The effects of these factors are automatically 
taken into account in most rectangular coordinate methods of reduction, owing 
to the fact that the differences between the measures of the positions of the stars 
themselves are involved in the reductions and take account of all disturbing 
variations of a linear character. In the present method no measures between 
the stars are made and the disturbing factors have therefore to be considered 
separately, and the corrections for them applied through the quantities A@ in (24). 
They prove to be serious only in the case of refraction. 

2.2. Changes of focus, which are generally small in photographic refractors, 
are liable to be appreciable enough in reflecting telescopes to alter the scale, that 
is the value of f in (12). Since however such a change affects all measured 


2 - 
° 
14° 20'- 
s/s 
-14°30'- - 
| 
D 
D, 3 
Zenith 
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distances proportionately, its only effect is to alter the size of the position triangle 
and, if necessary, it can easily be adjusted by changing the scale value (that is /). 
The correction to « resulting from this change in f, and any other small error in 
the assumed value of that constant, can however be neglected, since they will 
affect all values of g in a similar manner and, as has been shown, these do not 
require to be known to any high accuracy. 

2.3. Errors from a small tilt of the photographic plate will be too small to 
affect results within the limited area of a plate used. The images on plates taken 
with fast wide-angle lenses, and Schmidt cameras, are very sensitive to focus 
an¢ thus to tilt in the off-axis zones. Consequently the tilt in a properly adjusted 
instrument is unlikely to exceed a small fraction of a degree. Even in slower, 
longer focus instruments, such as the standard astrographic, we may neglect it, 
as the error caused is proportional only to the difference between the cosine of 
the angle of tilt and unity. Even a tilt of }° would at most cause an error of 0”:2, 
at right angles to the axis of tilt, for a value of @ as great as go’, and only a fraction 
of this error would appear in the resultant position of the object from three or 
more comparison stars. 

2.4. The orientation of the plate, which is a factor in some methods of 
reduction, does not of course concern us in this one. 

2.5. There are two external factors, aberration and refraction, which affect 
the positions of images on the plate. The maximum differential effect of 
aberration occurs in the direction of the tangent to the Earth’s path, that is, on 
the ecliptic 90° from the Sun, and might there affect @, in the unlikely event of 
it being as large as go’, by 0”-5. As however in the case of refraction at the 
zenith discussed later, this can be shown to represent but a small change of 
scale which is uniform in all directions in those areas and therefore affects all 
distances there proportionately. Like uniform changes of scale from other 
causes it can easily be taken into consideration in interpreting the position triangle. 
Elsewhere the differential effect is smaller and the error caused by it in the position 
of an object resulting from a combination of three position lines will be but a 
fraction of 3”. For practical purposes therefore we shall ignore this factor too, 
although by taking account of the ecliptic position a correction could fairly easily 
be applied. 

2.6. Differential refraction.—The differential effect of refraction however 
cannot be neglected. With 6=1° it may diminish the value of @ by 9” at 20° 
altitude, and even at the zenith by 1”. 

2.61. For the purpose of finding the correction \@ to be applied to the 
apparent distance 0, the zenith distance, if it be greater than 15° or 20°, and the 
position angle there of the zenith must be known approximately. If not 
available from tables or diagrams for the purpose at the observatory, such as 
the diagram here shown (Fig. 4)*, they can be found as follows. From the 
spherical triangle (Fig. 5) formed by the object at O, zenith Z, and the pole P, 
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* The curves for z and 7, plotted in H and 8, were computed by 
sin H =sin 7 s1n 2 sec ¢, (a) 
and, by Napier’s Analogies, 
tan $(90°— 5)=tan $(90°— ¢+2) cos 4(n +H) sec 4(n—A), (6) 
or, for numerically large values of (7—H), 


tan $(90°—8)=tan }(g0°—¢—2) sin +H) cosec $(n—H). (c) 
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H (Hour angle west) 
Fic. 4.—Distance and position angle of zenith. 
Diagram for the latitude of Oxford, giving z and n, simultaneously, with arguments H and 6. 
For hour angles east, reverse the sign of 7. 


and given the latitude of the observatory ¢, declination 5), and the hour angle # 
measured westwards, we calculate the zenith distance 2, by the cosine formula 


COS = sin sin + cos cos dy Cos A, (26) 


or more conveniently, and usually with sufficient accuracy using a three-figure 
table of cosines and a slide rule, by 


COS = cos — 59) — (I — cos cos cos By. (27) 


The parallactic angle, the position angle » of the zenith (measured from the 
pole eastwards), is then calculated by 


cos n =[sin (6 — + (I —cos cos ¢ sin 59] cosec 2, (28) 
or by the sine formula 
sin = Sin COS COSseC Zp. (29) 


Generally only a low order of accuracy is required here, unless the altitude is 
very low, and the three-figure accuracy will be ample. For the same reason we 
may neglect the small differences between the apparent and true angles in the 
above equations. Indeed, within 30° or 40° of the zenith it may be sufficient, as 
will become apparent, to estimate the required angles, z) and », at the telescope ; 
and beyond that to 60° or 65° an accuracy of I° will suffice. 

2.62. Consider now the effect of refraction on 6. In Fig. 5 let S, be the 
observed position of the comparison star at distance 6, from the object at O, 
p its position angle from O, z, and 29 their respective observed zenith distances, 
Ay, the difference of their azimuths, and S,’, O’, the positions cleared of the 
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H 


effects of refraction, with corresponding zenith distances z,' and 2)’, and 0,’. 
We shall omit the subscript in 0. 
From the spherical triangle OZS,, 


cos 6 = cos 29 COS 2, + SiN 3g Sin 2, cos Ag, (30) 
and 
— + 48in sin z, sin? 


where for present purposes @ and therefore (z)—2,) have been assumed to be 
small angles and the half angles have been written for their sines. Since the 
azimuths are unaffected by normal refraction we have, similarly, for the angles 
freed from refraction 


— 3,')?+4 sin 2,’ sin z,' sin? Ap). (32) 
The effect of refraction is always to decrease the zenith distances; therefore it 
makes the spherical triangle O’ZS,’ smaller in all cases, and decreases 6’. The 
correction Aé to be applied to @ derived from the measured s, to give 0’, is therefore 
always positive. 
The mean refraction zenith distance is given, with sufficient accuracy to 75°, 
from the apparent zenith distance z, by the well-known formula 


R, =58":294 tan —0"-06682 tan’ z, (33) 


and we will write this 


) 
R,=rtanz, 


where (34) 
r=K-— K’ tan*z, 
and, in circular measure, K =0’-000 2826, K’ =0"-324. 10-8, 
Consequently we have 
Zp 2,'=2,+7, tan (35) 


Since refraction does not exceed a few minutes of arc, except on the horizon, 
we may, in what follows, neglect r* and (A@)? terms, and make the usual 


approximations for small angles. ‘Thus from (35) 
sin 3) (I+79) Sin Zo, 


and similarly sin z,’, (36) 
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and, after substituting this and (35) in (32), and using (31) to eliminate 
Sin Sin 3, sin? we obtain 
20. AO = (0? — +71) + tan — 7, tan 2,), 
where we have put (37) 
€= — 2}. 
To simplify this we use the mean apparent zenith distance 
Z=}(%+%), (38) 
and expand both tanz in terms of €, by Taylor’s theorem, namely 


tan =tan (z+ 
=tanZ+}€sec?Z+ (39) 


and similarly for tanz,=tan(Z—4£). Then to the second order in é for the 
K terms, and the first order for the K’ ones, (34) gives 


ro +r, = 2(K—K' tan*2Z), (40) 
and, by appropriate cross-multiplication of terms, also 
ry tan 3y—7, tan 2, > sec? Z— 3K’ tan? Zsec? 2). (41) 
Substituting these in (37), we get 


2 
2+ 3tants) | tants tants (42) 


As the last term does not reach half a unit of the fifth place at s=75°, and in 
any case represents only a minute uniform change of scale over the area measured, 
it can in practice be omitted. (42) may therefore be applied in the convenient 
form 
0K (1+ 50), (43) 
where V is tabulated in Appendix B with mean zenith distance as argument. 
2.63. It has been accepted that the formula for refraction given in (33) is 
sufficient to zenith distances of 75° and therefore (43) may be taken as 
adequate to that limit. The values of K and K’ are however for the mean refraction 
at a standard barometric pressure (P) of 30 inches, and temperature (7') of 
50 deg. F. At other values the K used in(43) needs to be multiplied by the fraction 
197 
400+ T° (44) 


A double-entry table for K for this purpose is given in Appendix B. It is 
assumed that K’ is affected in the same way as K, so that the ratio K’/K remains 
constant and V is unaffected by changes in temperature and pressure. 

2.64. The shift of all celestial objects towards the zenith by refraction is 
obviously equivalent, at the zenith, to a uniform shrinkage in the scale, and (43) 
in fact reduces there to 

= (Z=0). (45) 


The tan*z terms in (42) are not appreciable within 10° of the zenith, and for 
most practical purposes one will be able to ignore them as far as 20° and treat 
the scale as the same in all directions within that distance of the zenith. For 
the rest, if a line be drawn on the diagram through the position of the object 
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and at right angles to the position-angle line of the zenith, it will be parallel to 
the horizon, and é may be measured from it for each comparison star with sufficient 
accuracy (see Fig. 3). 


Alternatively, since, nearly enough, 


cos (p—n) = €/6, (46) 
we can in each case read the angle (p— 7), between the direction lines of the star 
and zenith, from the diagram and then use 
Aé = 6K[1 + V cos? (p—)}. (47) 
2.05. In this graphical procedure the change in the angle (p—7) due to 
refraction will be insignificant. That moreover the error in obtaining and using €, 
or (p—»), as if the area was a plane surface, may be neglected may be shown by 
solving the spherical triangle OZS, (Fig. 5) for = 2)—2,: 
COS 21 = Cos cos + sin sin cos (p — 7) (48) 
and 
= cos(p—n) — £2) cot 29. (49) 
By differentiating (43) with respect to €, it is found that the maximum error 
in A@, caused by neglecting the last term in (49), occurs when €=6/4/3. For 


zenith distances to 75°, with @ as great as go’, the maximum error is less than 


2.7. Deformation of curved photographic plates, and optical distortion 

2.71. We have assumed so far that the distances measured on a flat plate 
(or film), after its development, exactly equal what they were at the time of 
exposure; and that, in the case of a plate which was spherically curved at exposure, 
the subsequent measures of distance on the flat plate equal the lengths of the 
great circle arcs they were on the curved plate. This assumption clearly cannot 
be true, and simple geometrical considerations indicate that the deformation 
must be appreciable in terms of the accuracy we aim at. 

A knowledge of what happens to the positions of images, on plates which 
were deformed at the time of exposure and then allowed to recover their flat 
shape, requires physical investigation with the materials used. So far as is 
known to the writer no such investigation has yet been made. However, since 
this paper was originally prepared, Dr W. M. Shepherd of Bristol University 
has kindly informed the writer of a preliminary analysis he has made, applying 
the theory of elasticity to the problem; and an assumption based on his results 
will be made in the course of what follows.* 

2.72. Let it be assumed that when a plate (or film) is pressed against the 
spherically-shaped former it is everywhere substantially in contact with the 
former and that the surface assumes the spherical shape of it; and that, moreover, 
the plate deformation is symmetrical about the plate centre. Both square and 
circular photographic plates, and circular films, are at present used in Schmidt 
cameras; and, although these assumptions are justifiable for a circular plate or 
film, the case of the square plate will need separate investigation. 

In the derivation of the equation (12) for @ it was assumed that the measured 
distance s was the true measure of the great circle arc between the images. Let 
us then now assume As to be the correction to be applied to s, measured on the 


* His paper on the subject has since been completed and is published, immediately following 
this, on p. 450. 
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flat unstrained plate, to give the true arc between the images on the plate when 
spherically curved, which is the distance that should be used in (12). Similarly, 
denoting the radial displacement by u, let (g, + u,) and (gq + ug) be the corresponding 
great circle distances from the plate centre on the surface of the sphere. From 


the spherical triangle S,CS, (Fig. 1), and denoting the angle S,CS, by C, we 
derive by the cosine formula 
sin? = sin? $(A, — Ag) + sin A, sin A, sin? $C, (50) 
where instead of equations (1) and (8) we now have 
= (6+ As)/r, (51) 
A, =(g1+%)/7, similarly A, (52) 
and, for the ordinary Schmidt camera, r=f. 
The equation for the plane triangle on the unstrained plate, as measured, 
corresponding to (50) is 
= (81 — 82)" + $C. (53) 
Since radial symmetry in the deformation is assumed, the angle C is unaltered, 
and we may use this equation to eliminate sin? $C in (50), which then becomes 


sin? = sin? $(A, —A,) +s? ] 2. (54) 
This may be simplified, in order to obtain an expression for As, if we make the 
same assumption as will be made in the working formulae derived in Appendix A, 
namely that the angular distance between the measured images does not exceed 
and if the values of u are very small. Then the cube and higher powers of 
A, and (A, —A,), in the expansions of the sines of the half angles, and also (As)? 

and u?, may be neglected. Using on and (52), (54) then reduces to 


r? sind, sind, 
£182 

2 
=I- E — + (Art + Ast) | 

+2). (56) 

In this it has also been assumed that A will not exceed about 7°, so that the fourth 
and higher powers of A, and the uA? terms, are negligible. 

2.73. Shepherd’s theoretical analysis of the problem indicated, on certain 
reasonable assumptions, that in a circular plate there would be a radial stretching 
of the material, which would be uniform in the case of a celluloid plate, and very 
nearly so ina glass one. For present purposes we will assume that this conclusion 
is substantially correct and that, to the order of accuracy required, the radial 
displacement in both cases is given by 


(57) 
where y is constant for any given plate. According to Shepherd’s results, y equals 


the radius of the field covered by the instrument, that is, the maximum value 
of g/r; and the value of u/r at the edge of the field, 4°, is independent of the 


where 
2G=I1-— 
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material, and amounts to 2”-0 for a field of 34° radius. The interpretation of 
y is, however, not material, and the only assumption to be made here is that it 
is a constant. 


Assuming (57) then, we obtain from (55) and (56) 
S ~g.\2 

— +A?) E } (58) 
Since the 7* term represents a change of scale, and will be constant for a given 
instrument, it can be omitted from the As here given and included as a correction 
to the scale constant (see equation 60, Appendix A.2). The remaining term 
then represents the difference between the measured distance and the corre- 
sponding arc on the sphere. This term is included in an approximate formula 
for 6 in Appendix A.4. 

2.74. It will be noted that, if y is proportional to the radius of the field, 
(57) implies that in a square plate there may be about twice as much stretching 
of the plate along the diagonals as there is along the short radii vertical to the 
sides, and hence probably there is distortion on account of which the angle C 
in each case does not remain unaltered when the plate returns to its flat shape. 
Equations (55) and (58) cannot therefore be accepted as satisfactory, for cases 
where a square photographic plate is used, until justified by further investigation. 
Meanwhile the safest procedure, when such a plate has to be measured, would 
seem to be to select the comparison stars so as to keep the angle C as small as 
possible; and this will also make the last term of (58) small. As it is unlikely 
that the differences between the displacements (wu) involved will exceed about 
}”, this procedure should in practice be adequate. 

2.75. Errors due to optical distortion and other residual aberrations have 
not so far been mentioned. Here too we may assume that the displacements 
of the images are symmetrical with respect to the optical axis. Other methods 
of reduction assume these effects are taken up in the plate constants, and 
therefore in the present method we assume they are small enough, at any rate 
over the limited area of the plate measured, to be absorbed by a linear change of 
scale which will not affect the position deduced from the position triangle. 
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APPENDIX A 


Approximations and computational modifications 

A.1. The method of reduction given in the paper has been developed for 
the general case, and the leading equations* can be applied directly as they 
stand to find the photographic position of any object. ‘The computer may 
however prefer to use formulae which are arranged to facilitate the numerical 
work, and these may readily be derived using approximations which are valid 

* Equations (2), (9), (10), (12), (23), (24) and (42). 
31 
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for the practical cases which may be met. Thus the angular dimensions of the 
field covered by existing telescopes do not exceed 15° square (the Franklin—Adams 
camera) and are normally less than 8° square. Further, there are few places in 
the sky where suitable comparison stars cannot be found within one degree of 
an object for which the position is required, and one may therefore also assume 
that @ will not exceed go’ of arc. Moreover, owing to uncertainties in the 
catalogue positions of comparison stars, it will be sufficient to ensure that errors, 
arising from the modifications now to be introduced, shall not affect the final 
result by more than 0”-I or 0”:2. 
A.2. Under these circumstances equation (12) takes the form 


= Vcos $1 COS E — (sec —sec (59) 


for which the error in 6 will be at most 0:05 at 6 =60' and 0”-15 at go’. Further, 
provided r is not appreciably less than f we may, from (2) with the same error, 
substitute s for s,. Introduce now a scale constant o defined by 


cosecI’ —3437°75 
=> = , 60 
f ( ) 


so that measured distances may be converted to angles in minutes of arc, and 
equation (59) may then be written 


—_— f 2 12 
> oSV COS COS Pg (sec — sec | (61) 
In this expression 


COS COS }(COs + COS I — }( py + Pa) tan? 3( — 


or 


V cos COS 4(cos + COS Yo), (62) 
with an error in 6, at most, of only 0’-03 when y=10° and 6=go’. Further we 
can write 


2sin msin? 
SEC — SEC = sin sin }( — P2) 


cos (63) 
where P= 3(P1 + 
But, from (62), to the same order of accuracy, 
COS COS Py = Cos? Cos? — (64) 


so that 
SEC — SEC Po 2 sec Ptan Psec — Pq) tan — Po) 
= (v1 — Sec tan (65) 
and the maximum error due to this approximation, for our limiting values of 
¢ and 9, is less than 0”-02. Equation (61) thus becomes 


—@,\2 


Since (a- Y2)/os attains a maximum of 4/2 (in the case of a Schmidt), and 
since @ is less than 10°, we may expand the square bracket by the binomial 
theorem, and - 


— = (cos P1+COS Po) E sec? tan? p 


P2 
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the neglected terms contributing, for our standard limitation on 6 and gq, an 
amount less than 0-05. Multiplying out this expression, and substituting for 
COS P, COS Py in the second and third terms the value 
}(cos +Ccos =cos }(g,— Ya) = cos |, (68) 
we finally obtain 
0 ~os—as [ (2-4) +N’ ] 
os os 
where 1 (69) 
H, =hav ¢,=4(1—cos Hy=hav }(1 —cos 


N=h4sec ptan®? and N’= }sec* ptan*y, 
are tabulated, with argument qg, in Appendix B in Table I. The N’ term is 
negligible in most cases. 

A.3. Equation (69) can be simplified still further, when g does not exceed 6° 


and %— ¢ is numerically less than 85’, as will usually be the case, by putting 


N= H,+H,, (70) 

which it is found is sufficiently accurate (0”-1 in @) within these limits. Then * 

(73) 


A.4. In the case of the spherically curved photographic plate we may 
conveniently modify this equation so as to incorporate the correction As to s 
on account of the deformation of the plate, assuming the expression (58) is 
substantially correct. If the term representing the change of scale be omitted 
in (58), as there suggested, and the scale constant o in (60) be modified accordingly, 


we shall have 
2 


Now from (8) and (10) 


tan sin },, (73) 
and expanding this, and neglecting cubes and higher powers of 9, 
Am xr (74) 


Further, since (A,?+A,*) is a small quantity we may neglect u inside the square 
brackets of (72), and in place of (52) use 


(75) 
Substituting these in (72), 
— ¥2) being in angular measure. But from (69) 
(77) 


* By tabulating slightly adjusted values of H, this equation could be made serviceable as far 
as about = 8° when @< 50’. 
31* 
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so that, neglecting y* and higher powers, we can write * 


Adding this expression to os in (71), and neglecting the a of (H,+H.), 
we obtain 
\ 
6 os(H, + 4 + + (:- 3 (79) 


(y<6°); and where, from (9), f/r=«?—1I. 


For plates of the ordinary Schmidt camera, where r=f and x?=2, this reduces 
to the simple formula 


os—os y+ H,)| + "1, (80) 
(k= 1/2, 


A.5. It will not be necessary to compute ¢ in each case by (10). We can put 
=Kog— (81) 


and tabulate the small correction g’, once and for all, with argument xag. 
Table I in Appendix B includes, for this purpose, g’ for the case of the flat plate, 
and for the ordinary Schmidt. To compute g’, gy was obtained direct, in the 
case of the flat plate, by equation (13). In the case of the convex curved plate, 
using (9) to eliminate 7, (10) was re-written 


4 
tan p=25—— sin og, (82) 
which, for the case of the Schmidt plate, gave 


tan p=27/2sin jog, (k= 2). (83) 

Other cases may be treated similarly, and tables made for the particular type of 
telescope used. 

A.6, The accuracy to which g must be measured may be deduced from (20a). 

If we lay down the condition that the effect on @ of an error in the measurement 

of each g of the pair shall not exceed half that for the estimated precision of s, 


we shall have 
Ag. £i-8 As 


This may be simplified as before by putting (g,—g,.)/s=1, its maximum value 
which occurs when the two images lie on the line joining them to the plate centre. 
In that case too one may assume, in multiplying the main bracket by s, that 


which is true for a flat plate and a sufficient approximation in other cases. 


As may be given a value +0-003 mm, for in practice the accuracy with which 
star images can be measured is normally not better than this and indeed, for 


* Even when q>6° it will be accurate enough to use (78) to correct os in (69), provided the 
expression in square brackets in (72) is substituted for the ((~,—q@z) approximation of it in (78). 


9 
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off-axis images, trailed images and comets, it may be much worse. Thus we 
may safely write for the required accuracy in g 


Losec? 
Ag < 0°003 (86) 
2 
< 0-003 ¥ mm near enough. (86a) 


Substantially the same result is obtained more easily from the approximate 

expression for 6, equation (71). From it 
2 
Ag < 0-003 coms mm. (87) 

For the extreme case of a flat plate (ck?=1) with g=10°, these indicate that 
the required precision Ag is 0-05 mm; _ but most cases will fall within 6° of the 
plate centre and for that limit Ag is 0-14 mm. For a Schmidt plate (x? =2) with 
a field of 6° diameter, p=4°-2 at the edge of the field, and Ag at that limit 
is 0-14 mm. Thus it is clear that generally it will be sufficient to measure g with 
a good millimetre ruler. Moreover the measures need not be made from the 
exact centre of the plate (optical axis) so long as they are all made from the same 
point, for any small error in the assumed centre will affect all values of y more 
or less equally and will not appreciably influence the final result. 

A.7. Although in computing equation (23), for the distance p, it will not 
always be possible to substitute the angles for their sines, because at high 
declinations the right ascension differences may be large, the work can be 
simplified by calculating it in the form 

(py — Sq)? = COS COS 34 (%q — — S,)? + (89 — 8, — So)’, 
where (88) 
S,=x—2sin}xcosect’; x=p, (%)—%,), (89 etc., 


in minutes of arc, and S is given in Table III. S can be neglected provided p does 
not exceed go’, and the right ascension difference does not exceed 6 min of time 
whatever the declination, or if the declination does not exceed 40° or 45°. It 
will be noted that the errors from this approximation, since they occur on both 
sides of the equation, always partly cancel out. 


APPENDIX B 


Formulae and tables 
Data 
Comparison star positions (x, 5,, etc.); measured distances from object 
(s,, etc.); approximate distances from plate centre (g,, etc. and gy for object). 
Scale constant o =3437°75/f; instrumental constant «=(1+/f/r)¥*; where 
f =focal length and r = radius of curvature. Also approximately the temperature, 
pressure, and hour angle. 


Angular distances (;, etc.).—Table I 
Po , Similarly Pa Ys. 
(a) y<6°, and | <85’. For a flat plate case 
ong ll 
6, =08, —o5,(H, + E + (ee) ] similarly 63, 


oS, 
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and for a bent Schmidt plate (measured flat) * 


6, = 05, —05,(H) +H) E + (4) ] , similarly 45. 


os, 
(6) p>6°, for a flat plate case 


—q,\4 
0, =05,—o5, He +H,+N, +N,’ (24) ] , similarly 65. 


1 oS, 
TABLE I 
Quantities depending on g 
g’ with arg. «og; H with arg. p; N and N’ with arg. P (mean ¢) 


H 10° Hx108 | Nx10° | N’x10° 
q 
K=+/2 
° ° 300 o'8 190 384 I 
13 26 
10 310 203 410 I 
13 28 
20 I 320 216 438 I 
30 2 330 1°O 230 | 466 I 
4° 3 340 rl 495., I 
50 5 350 1:2 259 524 I 
15 31 
60 8 360 1°3 2 
70 10 37 1°4 289 2 
80 14 380 rs 305 620 2 
go 17 390 321. 653 2 
100 21 400 18 338 688 2 
4 17 35 
110 26 410 355,, 3 
120 o'l 30 420 759. 3 
130 orl orl 36 430 2:2 391 797. 3 
140 orl 41 440 2°4 409 835. 4 
150 48 450 428 874 4 
19 40 
160 orl o'l 54 460 a7 447, 914 4 
170 o'r 61 470 497 955,, 5 
180 o'2 69 480 487, 997, 5 
190 76 490 507, 1 040 5 
200 o'2 o'2 85 500 a5 528° 1 084 6 
21 45 
210 93 510 1129, 6 
220 102 520 3°9 571), 1175), 
230 112 530 1 222 8 
240 122 540 4°4 616 1 270 8 
250 O-4 O'5 132 550 4°6 639 I 319 9 
260 143 560 4°9 I 9 
270 154 570 686 1 420 10 
280 166 580 5°4 710. II 
290 178 590 5°7 735,, 1 525. 12 
300 08 o'8 190 600 6:0 760 1 579 13 


* Subject to certain assumptions. See Section 2.73 (p. 440). 
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Approximate position (x, 59) 

Plot comparison stars: 0, 5,; 52; 33; 
where cosé is an average value. With radii 0,, 0, 0, draw arcs to intersect. 
Read Aa’ and 8, of the intersection 


Aa®=4s8ec8. Aa’; a=a+Aa®; 8. 
Differential refraction correction (A@).—T ables and IIB 
(a) z<20°. A@,=0,K, similarly Ads. 
(6) z>20°. ¢=latitude, a=hour angle, z=zenith distance. 
cos z = cos (¢ — 59) — (I — cos cos cos 8p, 
cos 7 =[sin(¢d — 49) + (I — cos #) cos ¢ sin dy] cosec 2, 
or siny=cos¢sin cosec z. 
Draw in the parallel of altitude (7 + 90°), and read off € =2)—2 for each star. 
AO, = 6,K(1 + Vy €,?/0,?), similarly 
TABLE ITA 


Refraction factor K x 10° 
Arg . temperature and pressure 


Barom. Temperature Barom. 
pressure pressure 
mm_ |—40°C —30°C —20°C —10°C}] o° C |+10°C 20°C 30°C 40°C in. 
—40°F —22°F —4°F +14°F] 32°F 50°F 68° F 86°F 104° F 


500 22°5 216 207 20°0 18°5 17°9 16°8 19°7 


520 23°4 22°5 21°6 20°7 20°0 19°3 18-6 18-0 17°4 20°5 
540 24°3 234 20°8 20°0 19°3 18-7 18-1 
560 24-2 23°'2 ars 20°8 20°! 19°4 22°0 
580 25°0 24°1 ans 20°8 20°! 19°5 22°8 
600 270 249 239 21°5 20'1 23°6 


620 23'°8 23° 2273 20°8 24°4 
640 28:8 276 =. 265 25°5 24°6 239 229 222 25°2 
660 274 26°3 25°4 24°5 23% 22°90 26°0 
680 306 2904 $2832 26°1 244 236 228 26°8 
700 ors 30°2 29°0 27°9 26°9 25°% 24°2 23°5 27°6 


720 29°99 28-7 aa 26°7 25°38 24°9 24°1 28-3 
749 | 32° 30°7 | 284 | 274 265 25% 248 | 2971 
760 942 3375 307 29°2 23:2 272 29°9 
780 35°! 339 32% 30°0 28:9 27°99 26:2 30°7 
800 | 36:0 345 332 319 | 308 | 207 28:7 27°77 268 | 31°5 


Calculated distances (p) and intercepts (D).—Table III 
(py — S,)® =(% — — S,)® cos 5p cos 5, + (89 — 8; — Ss)’, 
D,=6,+ A0,—p,, similarly D,, Ds. 
Position-line correction 
Plot intercepts, draw position triangle. 
Read corrections Az)’ and Then Ao’ = 4 sec dy. Ax’. 
Corrected « =a» + Aap’, 5 = + Ady’. 
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TasLe IIs Taste III 
Refraction factor V 


z V z V x S, x 10% 
10° 003 56° 218 50’ I 
15 07 57 2°35 
20 2. 58 2°53. 60 I 
25 “22 59 2°74 7° I 
30 60 2-96 
35 49 25 
40 61 3°21 100 4 
62 3°49. 
41 63 3°79 110 4 
42 “81 64 4°13 120 6 
43 ‘87 ° 65 4°52" 130 8 
44 43 140 9 
45 0-99 4°95 150 12 
67 5°43 
46 1:07 68 5°98 160 14 
47 1'14 69 6-61 170 17 
48 1°23 : 70 7°34 ” 180 21 
po 7 8-17 28 
50 1°41 I 
10 72 
51 73 10°28 210 33 
52 1°63 74 220 38 
53 175 75 13°23 230 43 
54 1°88 240 49 
55 2°02 a 250 55 
260 62 
270 69 
280 77 
290 86 
300 95 
APPENDIX C 
Example 


The data for an image of the minor planet Bruchsalia (No. 455), about 54° 
off-axis on a plate taken on 1952 May 15-86875 with the Franklin-Adams 
camera at the Union Observatory, Johannesburg, were kindly communicated 
by the Director. The values of the distance s in millimetres were calculated 
from the rectangular coordinates given, and the values of g were deduced from the 
equatorial coordinates of the stars and the plate centre (« 15®50™, 5—19°). The 
data required for the present purpose were then 


BD star 1950-0 81 950-0 g 

h m 8 4 ca 
(1) —14° 4343 16.02 23°58 —14 57 10°71 8-062 99°4 
(2) —13 4337 16 03 43°49 —14 10 34°8 9°867 115°8 
(3) —14 4351 16 04 36°49 —14 37 27°4 4984 III‘l 


(Star places from Trans. Yale Obs., 12) Planet 107°4 
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The scale value was stated to be 1 mm=182":1, that is, o=3'-035. The 
computations and position diagram (Fig. 3, p. 434) gave 


(1) (2) (3) (o) 
kog 301°7 351°5 337°2 
—g’ —o'8 —1°2 
300°9 350°3 33671 325°0 
H 10° 191 259 239 223 

Approx. position 
os 24°468 29°946 15°126 a 16 368-4 
— A(as) —0'200 —0°'247 —o-108 do —14° 40'°3 
6 24°268 29 ‘699 15'018 cos 0°967 39 
+239 —16°2 —14°3 27°°4 

Aé +0007 +0:008 +0:005 n 241°°0 

K (assumed) 0-000 24 
Oj — —72°82 +7:-09 + 60-09 V 0°27 
cos 8 0'966 14 0'969 55 0:967 60 
5;—8, 16-868 +29-720 +2843 
24°378 29°770 14809 
D —0'103 —0'063 +0°214 


From the position triangle 
Aa = —0':228 = — 08-94, Ad’ = +0'-083, 

and the resultant position was 

03M 358-46, —14° 40’ 13”:0 (1950°0). 
A recomputation of p showed values of D= +0-003, +0-006, and +0-005 
respectively, and yielded a correction of —0’-002 to 6 and the final position 

a 16" 358-46, —14° 40’ 
For comparison, the position was computed by the standard coordinate method 
which gave 


16 03™ 35°47, 8 —14° 40’ 12-9, 
and by the dependences method which gave 
a 16" 03™ 35844, 8 —14° 40’ 13-4. 


The differences shown in this last case are due to an error of 0”-6, arising from 
the neglect of the square terms of the coordinates which affect the dependences 
increasingly with the off-axis distance. Most two- and three-star methods of 
reduction would be subject to error from this cause if used more than a few 
degrees off-axis, but the position-line method is free from it. 


THE DEFORMATION OF A PHOTOGRAPHIC PLATE 
IN A SCHMIDT CAMERA 


W. M. Shepherd 


(Communicated by E. H. Linfoot) 
(Received 1953 May 4) 


Summary 


When a flat photographic plate is bent on a spherical former, as in the 
Schmidt camera, strains are produced in the plate which must be allowed for 
in making measurements on the exposed plate. The necessary corrections 
are obtained. In the case of celluloid, which has a Poisson’s ratio of about }, 
these take a particularly simple form. 


The photographic plate of a Schmidt camera is, in the unstrained state, 
plane; but when in position to receive images it is bent by being forced by a 
ring at its edge on to a spherically curved former. The plate may be circular 


UNSTRAINED 
PLATE 


STRAINED 
PLATE 


Fic. 1. 


or it may be square, but in the following discussion will be assumed circular 
and, when bent, to form a cap of angular diameter 2y. The angle y is small and, 
in the example considered, is taken to be 3°, the radius of curvature being 
60 in. and the thickness of the plate about 0-04 in., ie. I mm. 


3 
| | | 
| 
| N 
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After the photograph has been obtained the plate is removed and measured 
when it has returned to the plane state. The problem is to determine the 
corrections which must be made in these measurements, in order to find the 
corresponding coordinates at the time of taking the photograph. 

A number of assumptions are necessary. ‘These assumptions are: 


(1) The plate remains elastic and isotropic (in the elastic sense) throughout. 
(2) The plate is in contact with the former over its entire surface. 

(3) There is no friction between the plate and the former. 

(4) There is no radial tension in the plate at the restraining ring. 


Of these it will appear that the second is not strictly true and this point will 
be reconsidered later. ‘The third is certainly very nearly true since, as will be 
shown later, there is little normal pressure between the surfaces. If there is 
any doubt about the effect of friction the interface could be lubricated. The 
last condition is perhaps the most dubious. Fig. 2 represents the configuration 
with the plate in the strained position. O is the centre of curvature of the former 
and ON is the axis of symmetry of the system. ABCD is an element of the 
strained plate and the coordinates (r, A, 4) of the point A are as shown in Fig. 2, 
r being the constant radius of curvature of the former. 


Fic. 2. 


The equation of equilibrium.—The tensile stress across AB is T and across AD 
at Ais 7’. On resolving parallel to the tangent to the meridian at the mid-point 
of the element, for the equilibrium of the element, and proceeding to the limit 
it is found that 


sind +(T—T"’)cosA=0. 


In this problem the greatest value of A considered is 3°, so that with A measured 
in radians we have with ample accuracy 


(1) 
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Displacements and strains.—In the unstrained state the plate is a plane disk 
with its centre in contact with the former at N. The displacement of A may 
be described by the change in the length of the arc NA from rA—u in the 
unstrained state to rA in the strained state. 

The mean principal strains in the plane of the element are s,, s, in the 
directions A, ¢ increasing and 
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u 
0A’ 


I 
(2) 


2 
In obtaining these formulae =, (5) , A* are neglected. It may later be 


verified that these quantities are of the same order of magnitude and are in fact 
negligible. 

The normal pressure g between the plate and the former may be obtained 
from a consideration of the forces acting on an element of the plate in the 
direction of the radius of curvature. It is found that g is given by 


h 
q= =(T+T’), 


(3) 


where h is the thickness of the plate. 
Since h/r is, in the cases considered here, about ;4 it follows that 

stresses in a direction normal to the plate can be neglected in comparison with 

stresses in the plane of the plate. It also follows that, unless the coefficient of 

friction is remarkably large, frictional forces may be neglected. 

The stress-strain relations give 


= ZIT’ 


E, « being Young’s modulus and Poisson’s ratio respectively. 
These may be written, on using equations (2), 


T= E - |. 
rA 


E du (5) 


On substituting these expressions for T and 7” in equation (1) the equilibrium 
condition takes the form 


d fA du u 7 _u Az o du 
Rivas tna’ 


du. du (30 — 1) 


The solution of this differential ~— is 


u= = 


pe — 30)r3, 


(4) 
) 
|__| 


No. 4, 1953 


Since u=o at the centre of the plate, B=o leaving A to be determined from 
the conditions at the edge of the disk. These conditions are difficult to determine 


but a reasonable assumption to make is that Z'=o at the edge of the disk, 
i.e. when A=y. 


On substituting from equation (6) into equation (5) 7 is given by 


so that A takes the value (1 —o)ry?/16. 
The displacement u is then given by 


= — — — 301%. (7) 
The greatest value of wu occurs when A=y and then u,,,,=47 y*, which is 


independent of Poisson’s ratio. 


It is now possible to reconsider the assumption that there is complete contact 
between the surfaces. Equations (5) give 


E fiduu <A” 


and using this equation with equations (3) and (7) it is found that 
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The pressure q is positive if A<y//2, i.e. A<o-7y approximately, but becomes 
negative for A>y/1/2. However, the conditions near the edge of the plate are 
uncertain in other respects and it does not appear likely that further investigation 
of this point would be of value. 


The strains in the plane of the plate are given by 


s,= = (1-0) (1-30), (8) 


r2 


1 


The necessary correction « to the off-axis angle A is given by «=u/r and may be 
evaluated in radians from the formula 


(10) 


This quantity is greatest when A=y and it then takes the value y* 24, which is 
independent of Poisson’s ratio. When y=3°, the greatest value of ¢ is 1”-23. 

Celluloid plate.—If the plate is made of celluloid, Poisson’s ratio may with 
sufficient accuracy be taken tobe 4. This simplifies equations (7), (8) and (9) and 
the displacements and strains are given by 


u=pryA, 


Thus the bending of the plate is seen to produce a constant radial stretch in the 
ratio I:1+5,, i.e. 1: 1+ y*/24 accompanied by a transverse stretch in the ratio 
I:1+sy ie. 1:1+5y?— If A<4y the material is stretched transversely, 
but if A>}y it is compressed transversely. The values of ¢ (in seconds of arc), 
s, and s, are given in Table I. 


a= 
_ —30 
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Glass plate.—In this case the appropriate value of Poisson’s ratio is about }. 


Ne O 


TABLE I 
Celluloid plate with y=3° 
€ $4 X 104 
0” -00 I “14 
o”-41 
o” 82 
1°-23 I 


Using this value equations (7), (8), (9) become 


The values of e, s, and s, are given in Table IT. 


The corresponding values of ¢, s;, s, for other values of y can easily be calculated 


TABLE II 
Glass plate with y=3° 
€ s,X 104 
0”-00 1°29 
-46 1°24 
88 1°09 
1°-23 0°86 


from equations (8), (9), (10). 


The University, 
Bristol: 
1953 May 1. 


$4 X 104 
0°63 

—o'89 
—3°43 


$4 X 10! 
1°29 
0°76 

—o'81 
—3°43 


°° 
2° 


THE ACCURATE CALCULATION OF APPARENT 
PLACES OF STARS 


F. G. Porter and D. H. Sadler 


(Communicated by the Astronomer Royal) 


(Received 1953 May 4) 


Summary 


The second-order terms in the standard formulae for the calculation 
of apparent places of stars are developed in detail, and various methods 
are suggested by which corrections can conveniently be applied; the 
magnitudes of the residual errors in various procedures are critically 
examined. 

It is concluded that the present standard of accuracy (0”-o1) may best be 
obtained by limiting the range of the formulae to six months, and working 
backwards and forwards from the beginning of the year, instead of forwards 
only as at present. Recommendations are made, incorporating various 
other proposals, to provide a consistent system in which errors can be 
reduced to acceptable limits without a prohibitive increase in the amount 
of calculation. 


1. Introduction—The standard formulae for the calculation of apparent 
places of stars were derived at a time when the limit of precision was 0”I. 
However, they have been used for the calculation of positions to anominal precision 
of o”-o1. Present requirements for a real precision of 0”-o1, particularly in 
connection with P.Z.T. (Photographic Zenith Tube) programmes, make it 
desirable to examine the limitations of these formulae and to indicate how they 
can be modified or corrected. 

The method of calculation and correction must be adaptable to systematic 
computation and, in this respect, the Besselian form using “star constants” 
and “day numbers” is generally satisfactory. The technique of the actual 
calculation will be dependent on the available computing equipment, but the 
errors and limitations of any particular procedure can be deduced from those 
of the standard formulae. Consideration is here restricted to minor modi- 
fications and corrections to this form. 

The individual expansions of the rigorous formulae for the reductions for 
precession, nutation and aberration are well known (1); but there does not appear 
to be a comprehensive expansion of the combination of the three corrections 
(2, 3,4). Accordingly the full second-order terms in the standard formulae 
have been derived, and are given here for several procedures and forms of 
corrections. The expansions are given ina general form, and are independent of 
the inclusion (or not) of short-period terms of nutation and of the precise method 
of calculating the aberrational day numbers. 

Second-order terms only are retained, so that higher-order terms will have to 
be considered for very close circumpolar stars; this is done in Section 12. 
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2. Notation.—Except where otherwise stated or implied, all angles are 
measured in radians. The standard notation is used for the Besselian and 
independent day numbers, and for the star constants with the omission of the 
factor ;; in a,b,c,d. ty is the epoch to which the mean place of the star is 
referred; t=t)+7 is the epoch of date. 

Suffixes are used to indicate particular frames of reference, or quantities 
referred to them, as follows: 


O mean equinox at epoch fp, 
I mean equinox of date f, 
2 true equinox of date ¢. 


No suffixes are used in the expressions for the second-order terms, or in 
contexts where the particular equinox used is not significant. 

3. General principles.—Since proper motion (including orbital motion), 
parallax and aberration cause changes in the actual direction in which the star 
is seen, the corresponding corrections must be calculated with reference to a 
particular equinox and applied to the position of the star for that same equinox. 
On the other hand, precession and nutation are changes in the system of axes of 
reference and do not affect the direction of the star; they must, however, be 
applied to the actual position of the star. The order of application of the two 
systems of corrections is immaterial provided each is calculated rigorously. 

Nutation is always calculated from the longitudes of the Sun and Moon 
referred to the mean equinox of date, so that precession must be applied before 
nutation. It does not appear practicable to refer nutation to a fixed equinox and 
to apply it before precession. If it were, it would be possible to use the true 
equinox of ft, as a basis for the calculation instead of the mean equinox; since 
nutation changes slowly the second-order terms would then be smaller. 

4. Fabritius terms.—If a star (~,5) is subjected to a small movement (real or 
apparent) on the celestial sphere of magnitude s in position angle S, then the 
changes in R.A. and Dec. are given (to second order) by: 


(1) 


cos dAx=s sin S+s*sin Scos Stand, 
=s cos S — }s* sin? S 


The second-order correcting terms, which increase towards the pole, may 
be written in the form: 


(2) 


inR.A. +tand Ax Ad, 
in Dec. —3sindcosd(Ax)?. 


In this form they have been discussed by Fabritius (§), the corrections being 
applied to the apparent position of a star calculated from the first-order terms of 
the standard formula. This procedure is inadequate in principle because: 

(a) the corrections (2) are only applicable to the approximate solution of a 
single spherical triangle, and so cannot be applied to the complicated geometry 
of precession and nutation ; 

(6) the Fabritius terms for the sum of two separate corrections are not identical 
with those deduced from the sums of the first-order corrections; cross-terms 
from the combination of the two corrections must be taken into account. 
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5. Proper motion.—If proper motion from the epoch of the catalogue to that 
of epoch f, has been rigorously applied to the mean place, there remains only 
the correction yz from the epoch ft, to date. From (1) it is seen that the second- 
order terms attain a maximum value of }s*tand, so that neglect of these terms 
will give rise to a maximum error of }$u?7?tand. For u=1" the error is 
0”-000 002 47 tan 6 and is negligible. 

A more important source of error arises from the cross-terms with the other 
corrections; these are of order g or h x ursecd, or 0”-0001secd for ~=1", and 
t=1. They are treated as negligible, but should be considered if, in special 
cases, large values of + are used. 

Orbital motion may clearly be equally treated as a first-order correction only, 
and is assumed to be incorporated with proper motion. 

6. Precession and nutation.—The complicated movements of the poles of the 
equator and of the ecliptic may be expressed by the rigorous formulae developed 
in the standard text-books (1). These formulae may be expanded to give terms 
of the second order as well as the first-order terms in standard form. The 
former may also be developed by more direct methods, but it must be noted 
that the Fabritius terms do not enter in any simple or direct way. The full 
corrections may be written: 


in R.A. U(a,5)+usec8, | G) 
in Dec. U’'(a,6)+w’, j 
U=Aa+ Bb+E=f+gsin(G+)tand, | 


U'=Aa’+Bb' =gcos(G+z), (4) 


and the terms in u and w’ are listed in ‘Table I (Nos. 2-5), in the line “No 
correction”. The most important terms are Nos. 2 and 3, which can give rise 
to a position error of fg (0”-003 for r=} and 0”-008 for 7=1), and No. 4, position 
error }g*tan6 (0”-o01 tand for t=} and 0”-oo2tand for =1). 

7. Aberration.—Since aberration consists of a simple movement of position, 
the second-order terms consist mainly of the Fabritius terms in (2); another 
term enters into the magnitude of the shift s, through the geometry of the triangle 
of velocities. ‘The full correction may be written: 


in R.A. V(x,5)+vsecd, 
in Dec. 


V=Cc+Dd =hsin(H +<«)secé, 
V’ =Cc' + Dd' =hcos(H + «) sind +icos6, 


and the terms in wv and wv’ are listed in Table I (Nos. 6-9), on the line ‘‘ No 
correction’. The most important terms are No. 6, which can give rise to a 
position error of $h?tané (0”-oo1 tan4, for all values of 7). 

The correction for stellar parallax may be incorporated with that for aberration, 
as described by Richards and Sadler (6), without affecting this discussion. 

8. Combination of the corrections.—There are two practical alternative methods 
of combining the corrections. 

32 


J 
where 
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TABLE I 
Second-order terms in apparent places of stars, before and after correction of standard formulae 
Term No. I 2 * 4 5 6 
Character | }f? fe hfe te? dg? 
Ma t=} 0”:0033 0”0016 0”-0004 0” -000 8 0” 001 0 
By 46 84 42 9 19 10 
Aa cos 8 
Function of « cos(G+«) cosGcosa sin 2(G+«) cosasin(2G+«) sin 2(H-+ a) 
No correction +sin § —sind +2 tandsind +cos 6 +secd 
(a) ° —sind ° +cos 8 +cos 8 
(d) +sind —sind ° +cos 6 +cos 8 
(c) +sin 8 —sind ° +cos 8 +cos 8 
Ad 
Function of « I sin (G+) cosGsina sin?(G+ a) sin? (H+ a) 
No correction] o —1 +1 —2tand —tand 
(a) +sin 28 —cos*d +1 —sin 285 ° 
(d) ° —1 +1 —sin 26 ° 
(c) ° —I +1 ° ° 
Term No. | 8 9 10 II 
Character hi hi? fh 
1.. | 00010 0”:0008 0” 0002 0” 8 0”-000 8 
Mag. 
8 2 60 12 
Aa cos 8 
Function of « sin (H+ a) cos (H+ a) I 
Method (i) Method (ii) 
No correction ° +1 ° ° 
(a) —sind +cos? 8 —sin 26 
() ° +1 ° ° 
(c) ° +1 ° ° 
Ad 
Function of « |cos?(H+«) cos(H+«) I sin (H+ a) 
Method (i) Method (ii) 
No correction] +sin25 +cos25  —sin 26 —sind ° 
(a) +sin25 +cos28 —sin 26 ° +sin 
(d) +sin25 +cos25 —sin 26 —sind ° 
(c) +sin26 +cos2d8 —sin 28 —sind ° 


1 
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TABLE I (cont.) 


Term No. 12 13 14 
Character gh gh gi 
©” 001 7 0” 001 7 0” 000 7 
Aa cos 8 
Function of « sin (G+ H+ 2a) sin (G+ «) cos (H+ sin (G+ a) 
Method (i) Method (ii) Method (i) Method (ii) 
No correction +tan 6 ° ° ° + 
(a) ° +4$sin25 +4sin 25 —sin?d +cos* 8 
(6) ° +4$sin25 +4sin 25 —sin®$ +cos?3 
(c) +tan 6(1 —sin 5) ° ° ° +1 
Ad 
Function of |sin (G+ «) sin (H+ a) cos (G—H) cos (G+ a) 
Method (i) Method (ii) Method (i) Method (ii) 
No correction —secd +cos 5 ° —sind ° 
(a) —cos 6 +cos 6 ° —sin6d ° 
(b) —cos 8 +cos 5 ° —sind ° 
(c) —sec 6(1 —sin 8) +cos 6 ° —sin $8 ° 


Table I lists the second-order terms in Aa cos 6 and A8 both before corrections 
are applied, and after the application of corrections (a), (6) and (c). For economy in 
presentation, the terms are split into three parts—the ‘“‘ character’’, the function of a, 
and finally the function of 8 by which the product of the first two is to be multiplied. 
Methods (i) and (ii) give rise to identical terms, with the exception of Nos. 10, 13, 14. 


(i) Precession and nutation first. Using suffixes to denote the reference 
system to which quantities are referred, the full correction is then: 


in R.A. tp, + U4(%, 59) + 5g) + (u+ v) secd 


oV 
= Tp, + Uj(%, 59) + Vo(%q, 59) + (u+v) secd + +U’ ap? (7) 
in Dec. + 59) + Vo'(%, 5g) +0" 
U4'(%, 59) + 59) + +U +U>-+U (8) 


In this case the day numbers A, B, E must be referred to the mean equinox 
of date and the day numbers C, D to the true equinox of date. Strictly, the 
values of the precessional constants m and m should be calculated for epoch 
ty +47, but it is only for large values of 7 that this refinement is necessary. The 
additional second-order terms are tabulated in Table I (Nos. 10-14). 

(ii) Aberration first. The full correction is: 


in R.A. ty, + Vo(%, 59) + + V, 59 + V’) + (u+v) secd 


32° 
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in Dec. + Vo'(%, 59) + + Vi 59+ V’)+u' +0" 
= + (a9, 39) + Vo'(%o, 5g) + +V’ (10) 

Here, the day numbers A, B, E are referred to the same equinox as in (i), 
but the day numbers C, D must be referred to the mean equinox of tf). The 
additional second-order terms, which differ from those in (i), are tabulated in 
Table I (Nos. 10-14, with modifications to Nos. 10, 13, 14). 

The most important term is No. 12, which is the same for both methods; it 
can give rise to a position error of ghtané (0”-002 tan 6 for 7 =4 and 0”-003 tand 
for 7 =1). 

9. Methods of correction.—It is clear that, to the accuracy required at present, 
correction for terms Nos. 4, 6 and 12, which are proportional to tan 4, is essential 
for close circumpolar stars. Three forms of correction are described, though it 
must be emphasized that the most convenient will be strongly dependent on the 
techniques used for the routine calculation of the standard formulae. With all 
three forms, the residual error does not increase towards the pole. 

(a) Fabritius correction. Although inadequate in principle, the application of 
the full Fabritius corrections does cater for those terms which depend on tan 6 
and sec 5; it introduces several other terms which, however, are independent 
of sec 5. If, therefore, the first-order terms are Ax and Aé, the full reductions are: 


in R.A. Ax +Aa Ad (11) 
in Dec. 


(b) Modified Fabritius correction. A slightly better representation of the 
second-order terms can be obtained by a similar device in which the term f is 
omitted from Ax. The reductions to apparent place may then be written approxi- 
mately in the form: 

inR.A. (12) 
inDec. 
where 
A’a 
+ Ap+Bb+Cc+ Dd; (13) 


N is the “ nutation in right ascension” and p is a—m. This form of correction 

"is particularly convenient in those special cases where the Atkinson proposal (7), 
to refer all apparent right ascensions to the hour-circle through the mean equinox, 
is used. 

(c) Tabulation of corrections. For large declinations, secd may be replaced 
approximately by tan 4, with an error that vanishes at the pole; when this is done 
the three main terms may be written: 

inR.A. }tan?d {g?sin2(G+«)+2ghsin(G+ H+ 
=Jtan?8, (14) 
in Dec. (H+«)} 
=J'tané. (15) 
The quantities J andJ’ may be tabulated in a simple table for each year as functions 
of « and date. An interval of 0-5 in « (24 entries only are required) suffices, 


if 
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without interpolation except in extreme cases; part of such a table for 1956, 
based on ty = 1956°5, is givenin Table II. The table can in principle be used with 
any value of r; the magnitude of the entries increases with 7, but not very rapidly 
for 7 less than, say, 3 or 4. 

For each star, values of J and J’ corresponding to « are taken from the table 
for each date, multiplied respectively by tan?d and tané and applied as corrections 
to the apparent place calculated from the standard formulae. 

In the case of stars for a particular P.Z.T. programme, 5 may be regarded as 
constant and the functions of 5 may be incorporated into the table ; all second-order 
terms may thus be included if deemed necessary. 

It is of course possible to extend the tabulation of corrections to include the 
terms which are independent of tané (or secd). Corrections could be given 
for afew ranges of 5. This complication has been regarded as involving too much 
work and is not further considered here. 


TABLE II 


Tabulation method: values of J and J’ for 1956 


Right ascension 
Date h h h h h h h 
J, unit 08-000 o1, to be multiplied by tan? 
Nov. 6] +15 +16 +16 ... —13. —15 =—16 9 +13 +15 
16 14 17 18 II 14 17 7 II 14 
26 13 16 18 9 13 16 +4 9 13 
Dec. 6 10 15 18 6 10 15 ° 6 10 
16 J 12 16 — 2 7 12 — 4 +2 7 
26]; + 4 9 14 +2 —4 9 —-8 —2 +4 
36 o +s + o —5 .. — 6 ° 
J’, unit 0”-000 1, always negative, to be multiplied by tan 6 
16 5 8 12 24 22 19 — 1 3 5 
26 4 % 10 27 25 23 ° 2 4 
Dec. 6 2 5 8 30 29 26 ° = 2 
16], — 1 3 6 32 31 29 ° 
26 ° I 4 32 32 31 — ° ° 


10. Residual terms after correction.—A complete list of the residual terms, 
arising from all three forms of correction as applied to both methods (i) and (ii), 
is incorporated into Table I. ‘These are all independent of tand. The “ character” 
has been chosen so that it is the maximum contribution of the residual term; 
maximum values of the magnitudes are given for 7 = } and r=1. 

In Table III are given estimates, deduced from the terms in Table I, of the 
maximum errors that can arise through the use of the various procedures described 


462 Jj. G. Porter and D. H. Sadler Vol. 113 


above. Values are given for 7 = } and 1, and for ’=0°, 45° and go°; for the case 
where no corrections are applied the error near 5 = 90° is expressed as A + B secd. 

The maximum errors, in either Ax cosd or Ad for the whole range of declin- 
ations, are collected together in Table 1V. No attempt has been made to combine 
the errors in right ascension and declination, so as to give maximum positional 
errors. 


Tas_e III 


Estimated maximum errors in various procedures, unit 0” -oo1 


Declination No corrections Form (a) 
and 


Aa cos 6 


a 
I 3 


8 9 
7 


6+3°4secd 13+56secd 
4--3°4sec8 9+5-6secd 


7+3°4s8ecd 15+5-6secd 
8+5-6secd 


TABLE IV 


Maximum errors in various procedures, unit 0” 001 


Form of correction, t= } Form of correction, T= 1 


Method (a) (b) (c) (a) (d) (c) 


(i) 7 16 
(1i) 5 II 


11. Comparison of methods.—There are two basic methods and three forms of 
correction, while the range of 7 is at choice. Examination of ‘Tables III and IV 
indicates that, unless more elaborate methods of correction are introduced, the 
use of 7=1 implies unavoidable errors of 0"-010. Moreover, corrections must be 


Form (6) Form (c) 
r=1 | r=} r=1 
(ii) I 3 
(ii) 4 8 4 8 
(ii) 5 9 4 8 4 9 
s= @) 7 14 4 10 4 10 14 
(ii) 5 II 4: 10 5 5 II 
. 45° (i) 10 21 5 10 7 14 8 16 
(ii) 7 14 6 13 5 11 + 9 
(ii) 5 9 3 8 3 8 
| 
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applied for all declinations if larger errors are to be avoided. This strongly 
suggests that 7 should normally be restricted to the range —4} to +}. 

When 7 is not allowed to exceed }, errors can be kept down to 0-005 by using 
either of methods (i) and (ii), with an appropriate form of correction. Method (ii) 
is, however, superior to method (i) in that, for a stated maximum error, it may be 
used without corrections for a larger range of declination; thus, if an error of 
0”-007 is acceptable, method (ii) may be used without correction up to 6= 45°, 
whereas corrections would have to be applied throughout with method (i). 

The choice of method determines the form in which the day numbers, to be 
used for all star reductions, are tabulated; it is therefore much more important 
than the choice of the form of correction. Of the three forms of correction, (c) 
involves the least work for both systematic and non-systematic calculation, 
provided a table of J, .J’ is available; this again favours method (ii), since the tabu- 
lation of corrections is more efficient than with method (i). Moreover, the adoption 
of method (ii) does not prevent the use of the forms of correction (a) and (6), 
whereas the combination of method (i) and (c) is poor. Method (ii) is also superior 
to method (i) when large values of 7 are used and all second-order terms are 
evaluated ; this may be seen from Table I in which more terms vanish with method 
(ii) thin with method (i). 

An apparent disadvantage of method (ii) is that the aberrational day numbers 
C and D (as well as h, H and 2) must be referred to the mean equinox of ty, so that 
there will be discontinuities at the changes of epoch. This is not a serious 
inconvenience, as the day numbers A, f, g and G have discontinuities at these 
points in any case. 

The consequences of restricting 7 to the range —} to +4, and the general 
adoption of method (ii), are discussed in Section 15. 

12. Higher-order terms.—The only higher-order terms that need be considered 
are the third-order terms of precession and nutation; calculation shows that 
these may be neglected, to the order of accuracy considered here, provided 7 tan 6 is 
less than 100. ‘This suggests the possibility of allowing 7 to reach 5, or even 10, 
with the considerable advantage of fixed “‘ star constants”’ for 10, or even 20, years. 
The tabulation method of correction, including all terms, is the only possible one, 
and this will be quite practicable in the case of P.Z.'T. stars for which 6 may be 
considered constant; in fact, the complete second-order corrections appropriate 
to method (ii), but not to method (i), are linearly interpolable at intervals of 50 days. 
Although the literal expansion is rather lengthy it is not difficult to compute 
the correction table, at an interval of 0-5 in «, for a period of 10 to 20 years. 

The use of large values of + is undesirable for non-systematic computation. 
Care must also be taken to see that other sources of error (see Sections 5 and 8 (i)) 
are negligible or allowed for. 

13. Aberrational day numbers.—In an earlier paper (8), a method has been 
developed for the accurate calculation of the aberrational day numbers C and D, 
referred to a fixed frame of reference. By oversight, the frame of reference used 
was not specifically stated ; it is that of 1950-0 and the quantities in Tables I and IT 
of that paper are (deliberately) referred to that equinox. Method (ii), in which 
aberration is applied first, can be used for such a fixed standard equinox provided 
corrections for precession and nutation are applied rigorously afterwards; as 
shown in Section 12, this method is both practicable and attractive, at least for the 
special case of P.Z.'T. stars. 
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The day numbers so calculated may be reduced to any other equinox by the 
use of the simple reductions : 


AC = —Dcose.AA—Ctane. Ae, 16) 
AD= + Csece. AA, 
where \) is the sum of precession and nutation in longitude from 1950-0. The 
term in Ac may be ignored, with negligible error. 

In method (i), the equinox required is that of date and Ad is variable; in method 
(ii) the mean equinox of fy is required so that AA is a constant. There is little 
difficulty in incorporating these corrections into the calculation of C and D. 

14. Reduction from apparent to mean place.—For meridian observations it 
may be necessary to reduce from the observed apparent place to the mean place 
referred to some standard equinox. Although this reduction can be done in one 
step (see Section 12), it will generally be easier to reduce first to the nearest epoch 
of the beginning of the year. To a reasonable accuracy, good enough for most 
practical purposes, it suffices to use the apparent place instead of the mean place 
as the basis of the main reduction, which is of course applied with reversed sign 
to the apparent place to give the mean place. The second-order correcting 
terms are approximately the same as with the direct calculation, and may be treated 
in exactly the same manner. 

Specifically, if R, R’ indicate the first-order reductions from mean to apparent 
place, then: 


( 
8, 3) + K’, 17) 


=%q + R(x, 59) + K, | 
J 


and, approximately, 


=%,—R(x4,54)+K, (18) 


59 — R’ (44,84) + K’, 


where ~,, 5, is the apparent place and K, K’ are the second-order terms. 
Corrections for the latter, if necessary, are thus to be applied with the same sign 
as for the reduction from mean to apparent place. 

15. Consequences of restricting the range of +.—-For systematic computation 
the epoch can be chosen at the most convenient point—for instance, at the middle 
of the year for 10-day stars and (in the northern hemisphere) the beginning of the 
year for P.Z.T. stars. The calculation will be done by means of Besselian day 
numbers, and there is no difficulty in deriving these for any particular epoch 
required. 

On the other hand, non-systematic computation implies the use of the indepen- 
dent day numbers published in the national ephemerides; there is no simple 
method of changing the epoch and (failing publication of day numbers for several 
different epochs) the user will have no choice of epoch or method. The epoch of 
the beginning of the year is too well established to be changed, but restriction of + 
to the range —} to +4 is too advantageous to be discarded. This means that 
the day numbers for the first half of the year will have to be referred to the epoch 
of the beginning of the year, and those for the second half of the year to the epoch 
of the beginning of the next year. As there will thus be a discontinuity in the 
middle of the tabulations for each year, there can be no objection to the adoption of 
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the preferred method (ii). Users should soon familiarize themselves with working 
from and to the nearest epoch, instead of the preceding one. 

The provision of a small overlap and of tables of the correcting factors J and J’ 
(as in Table II, but based on the epochs of the beginning of the year) should ensure 
both convenience and accuracy. Since the error in Ax cos 5 or Aé cannot greatly 
exceed 0”-o1 for declinations less than 60°, it will almost certainly suffice to ignore 
the correcting terms for all stars with declinations between —60° and +60°, a 
range which includes all projected P.Z.T.s. The maximum error will rise to 
about 0-023 at declinations of + 80° and it may well be that such errors could be 
tolerated, provided that a simple means of applying corrections is available (as it 
will be in the tables of J, J’); however, the work of including the corrections in 
the calculations is not heavy, and they will have to be included for circumpolar 
stars. 

16. Conclusions.—The main purpose of this paper is to put on record the errors 
that can arise in the various procedures for the calculation of apparent places of 
stars. It is, however, appropriate to make specific proposals, especially in the 
light of recent I.A.U. recommendations. ‘The relevant resolutions, formally 
adopted at the 1952 (Rome) meeting, are : 


(i) that the short-period terms of nutation be included in all apparent 
ephemerides other than those of 10-day stars (9) ; 


(ii) that the day number A be multiplied by and expressed in seconds of 
arc (2, 10). 

The Atkinson proposal (7) will not now be adopted. 

It is accordingly proposed that as from 1960 January 0: 

(a) the day numbers published in the national ephemerides be referred to the 
mean epoch and equinox of the nearest beginning of the year ; 

(6) short-period terms of nutation be included in all day numbers tabulated 
at daily intervals, but that f’, g’ and G’ (and possibly A’, B’ as well) be given 
separately in addition ; 

(c) short-period terms of nutation be omitted from the day numbers given at 
intervals of ten sidereal days; 


(d) the star constants a and a’ be redefined as: 
4 
a=m/n+sin«tand, a’ =cos4%, (19) 
where m/n = 87 + 0-002 367T, (20) 
and 7 is measured in tropical centuries from 1900-0. 


(e) apparent places of stars be calculated from these day numbers by the 
formulae (or their equivalents) : 


+ Aa+ Bb+Cc+ Dd+ tan? 4, (21) 

8 =89 + Aa’ + Bb’ + Cc’ + Dd’ +J’ tané, (22) 

where ./, J’ are taken as zero for —60° <5 < + 60°, and elsewhere obtained from 
tables, with argument «, published in the national ephemerides. 

The maximum errors of (21), (22) at 8=60°, with J, J’ taken as zero, are 


08-001 4 in R.A. and 0”-o10 in Dec. It is the present practice to tabulate the right 
ascensions of 10-day stars with declinations greater than 60° to o*-o1 only. If this 
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practice is continued, thenJ may be regarded as zero for all but circumpolar stars 
(5 greater than 80°); in fact, at 8=8o° the error will then have a maximum of 
08-008. The error in treating J’ as zero increases more slowly, being 0”-013 at 
6=70° and 0”-023 at 8=8o0°; for many purposes it may thus suffice to extend the 
range of declinations in which second-order terms are ignored (provided, of course, 
that 7 does not exceed }). 

17. Correction for short-period terms.—When daily values of A, B, C, D, E are 
used, «, 5 in (21), (22) will contain short-period terms of nutation. When 10-day 
values of the day numbers are used, short-period terms of nutation will not be 
included, though certain short-period terms will enter through the aberrational 
day numbers C, D; these have been discussed in a previous paper (8). 
Corrections for short-period terms of nutation to the apparent places of 10-day 
stars may be obtained, without multiplication or interpolation, with the aid of a 
permanent triple-entry table with arguments : 

g’ (interval 0"-02); G’+« (interval 0-5); (various intervals). 
g’ and G’ will of course be tabulated for every day. A typical extract from this 
table is given in Table V. 


TABLE V 


Corrections for short-period terms of nutation 
G’ + (o!-o for Ad) 


Aa—f’ (08-001) Aa—f’ (08-01) 
20° 30° 40° 50 65° 70° 75° 80° 


8 
10 
13 
15 


(The complete correction in R.A. is f’ +(Aa—f’)) 


For illustrative purposes the values of Aé are the maximum values corre- 
sponding to G’ + «=0"-o, ‘The full range of G’ + « need not be given, and negative 
declinations are catered for by changing the sign of Ax —f’. One of the advantages 
of the Atkinson proposal would have been that there would have been no term f’. 
A table of this form now offers little advantage over the Besselian form of correction, 
if the star constants are available; however, it provides a simple means of 
calculating short-period terms for the Sun, Moon and planets. 

18. Acknowledgments.—'Vhe derivation of the second-order terms in Table I 
is a straightforward piece of elementary mathematics; but the chances of error in 
such a complicated reduction are large. We have to thank Dr G. M. Clemence 
for pointing out an error in the first draft of this table. ‘The second-order terms 
have now been derived independently four times by Dr Clemence, Dr E. W. 
Woolard, and the two authors; it is hoped they are now correct. We are grateful 
to Dr Clemence and Dr Woolard for this verification. ‘Thanks are also due to the 
Astronomer Royal and Dr Clemence for helpful discussions on the choice of 
epoch and allowable errors in the apparent places. 
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EXPERIMENTAL AFTERGLOWS WITH BALMER DECREMENTS 
RELEVANT TO ASTRONOMICAL SPECTRA 


Martin Johnson 
(Received 1953 March 19) 


Summary 


The electron-capture form of Balmer decrement expected in some 
astronomical spectra cannot be isolated by ordinary laboratory use of fast 
electrons; the afterglow of a low voltage, high current, radio-frequency ring 
discharge is therefore examined sufficiently far away from glass walls. Electron 
multiplier measurement is used, with calibration to absolute intensities, and 
precautions are taken to obtain the afterglow free from survivals of applied 
voltage. Pure hydrogen and hydrogen-neon and hydrogen-helium mixtures 
are investigated. The portions of the decrement resembling a Boltzmann 
distribution of excitation are distinguished from those which suggest a 
recombination law, and divergences from either are examined for self- 
absorption, overpopulation of low energy states, radiationless de-excitation 
of higher states, etc., in view of likely distortions of astronomical Balmer 
decrement, notably in red stars, Be stars, and the ‘‘ Balmer progression’’ 
of peculiar stars. 


1. Variations in Balmer decrement, i.e. in the gradation of intensity along 
the visible series spectrum of hydrogen, are commonly observed in astronomy 
and experiment, but their correlation with physical conditions has remained 
uncertain. In a classic instance the intensity ratio of H« to Hf in the laboratory 
of Merton and Nicholson (1) altered from 1-7 to 4-4 when measured by an 
identical method, and there has been no possibility of deciding whether such 
a change implies that recombination into upper electronic levels is being super- 
seded by excitation from lower levels or whether upper levels are becoming 
impoverished by radiationless losses. The Balmer decrement discussed for 
nebulae by Plaskett (2), Page (3), Berman (4), Baker and Menzel (5), for the 
solar chromosphere by Menzel and Cillié (6), for prominences by Woolley (7), 
for red stars by Joy (8), Merrill, Scott (g), and in Merrill’s (10) ‘‘ Balmer pro- 
gression”’, could serve more powerfully as indicator of physical condition if such 
experimental data could be reduced to fewer controlling variables. For,instance, 
when true intensity changes are extracted from the blending of oxide bands in 
pulsating stars, the large apparent fluctuations of Balmer decrement with phase 
can locate the hydrogen layers only if we know how far any decrements do obey 
Boltzmann laws of excitation or the capture laws of Cillié or of Menzel, or whether 
no simple law without distortion can ever be expected. 

Laboratory approaches to astronomical Balmer decrement must face the 
obstacle that their spectra depend in general on wall-reactions and on a few high- 
speed unidirectional electrons, neither having astronomical counterpart, whereas 
many stellar atmospheres possess far greater concentration of slower electrons ; 
hence a pure capture spectrum with no “‘quenching’’, or an excitation obeying 
equilibrium conditions, or even the simplest astronomical fluorescence of Rosse- 
land’s type with Woolley’s uncompensated cycles, is unlikely to be isolated under 
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ordinary laboratory conditions. Deviations from these simplified mechanisms 
could, however, become an important contribution of experiment to astronomy 
if the experimental conditions were not too complex. 

In this paper we make an initial approach to such contribution by exploring 
the possibility that a simpler than usual deviation from astronomical conditions 
may be obtained for Balmer decrements in an afterglow: the luminosity is that 
which survives when voltage has fallen below maintenance potential in a high- 
current low-voltage electrodeless ring discharge at radio frequency. The fact 
that the afterglow is a time-transient may itself illustrate conditions in novae, 
non-steady envelopes, and the eclipsing atmospheres of binaries; but even in 
relation to steady stars its low-speed particles offer less of a departure from stellar 
physics than the high-speed electron streams which have made most discharge 
tubes irrelevant to astronomy. 

We shall retain the usual term “‘ decrement’’ for the instantaneous intensity 
gradation along the spectral series, and use “‘decay”’ for describing time-dependence. 

2. The very brief afterglow of hydrogen has been less accessible than those of 
atmospheric gases; it may occur for any of the following reasons: (a) survival 
of decaying voltage after cut-off, due to circuit inductance, (6) continued re- 
excitation during decay of kinetic energy of particles or during slow escape of 
imprisoned quanta, (c) recombination of free protons and electrons. Inhibition 
of afterglows may occur by the ‘‘ quenching”’ familiar in chemistry where an atom 
with high energy is precluded from radiating, the energy escaping in “collisions 
of the second kind’’ (superelastic) with other atoms or with free particles. 

It is notable that the principal experiments on hydrogen afterglows hitherto 
(Rayleigh, Craggs and Meek, Craggs and Hopwood (11), Ravenhill and Craggs 
(12), Lee and Fowler (13)) have variously ascribed the decaying luminosity to 
metastable states, to delayed re-excitation by hot gas, to propagation of super- 
sonic waves, as well as to a recombination process of the kind introduced into 
astronomical transients by Grotrian (14) and by the present author (15). 

The major discrepancies between astronomical and laboratory spectra arise 
firstly from the practice of using narrow tubes to intensify the light, which also 
intensifies wall-interference with gaseous processes, and secondly because the 
necessary isolation of Balmer lines free from molecular bands has commonly 
been obtained through high voltage. Other techniques, such as Merton’s 
dilution (16) with inert gas, Wood’s wall-poisoning, and Herzberg’s ring glows 
(17) driven by condenser impulses, are scarcely powerful enough to generate 
measurable afterglows, their intensities not competing with the recent condenser- 
flashing by Ware (18) or even the older high-energy discharges of Crew and 
Hulburt. A mechanism permitting high Balmer strength through high current 
but low voltage, and distributed over a large wall-free cross-section, will be needed 
to approach even remotely a stellar atmosphere; the requirements seem met 
by the discharge developed in Birmingham by Dr J. S. Symonds, following 
Thonemann and others, as a proton source for nuclear research. I am indebted 
to him for help in obtaining the gaseous and circuit conditions suitable for adapting 
to the present problem. 

3. A 350watt RCA oscillator is arranged for tuning from 3 x 108 to 14 x 10% 
c/sec, its output coupled by a loop to the inductance containing the hydrogen 
tube (Fig. 1) and its parallel variable capacity. At the resonance determined by 
pressure and frequency a Balmer luminosity about 15cm long by 6cm diameter 


470 Martin Johnson, Experimental afterglows with Vol. 113 


becomes very intense, and the axis of this electrodeless and wall-free plasma is 
focused through a revolving shutter on to the first slit of the spectrograph. 
A commutator of 8 segments, on the same motor shaft as the shutter disk 
running at I 500r.p.m., biases the output valves to cut-off for a millisecond 
at a time; during the central approximately 90 per cent of this interval the 
shutter slots admit the afterglow to measurement. 


350 watts at 3 


ohoto-rmultiplier 


galv. 
CRO beam | 
~ CRO beam 2 


FIG. 1. 


Double-beam cathode-ray pictures of the fall in energizing voltage and the 
decay of output from the spectrophotometer indicate that primary ionization 
ceases before the light enters the slit, having regard to the literature of low-voltage 
hydrogen currents. The narrowness of the gap between starting and extinction 
voltage is checked by experiments with the variable capacity. Care in phasing 
the shutter to open only at cut-off is also checked by frequent slow-motion re- 
examination of the light transmission. The hydrogen is admitted through Pd and 
liquid-air trap, and the neon and helium subsequently used are admitted from 
flasks sealed in through capillary lengths for measured expansion of each gas. 
The discharge in pure hydrogen extinguishes at about 0-08 mm pressure at the 
frequency most often used (3-2 Mc/s), but is steady at 0-1 mm after running long 
enough to exhaust initial adsorptive losses, and measurements are made during 
times of stabilization. Less brilliant discharges down to 10-3 mm were obtainable 
by retuning up to 14 Mc/s. 

The spectrophotometry was initially photographic, the Balmer decrement 
being obtained by microphotometer measurement of panchromatic plates which 
had received several hundred thousand such afterglow admissions. Calibration 
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was by means of a substandard GEC tungsten strip lamp with known temperature- 
current relation, and the characteristic curves of the plates for the different 
spectral regions. Accuracy suffered from the fact that a group of Balmer lines 
may reach in their intensity gradation to both extremes of the microphotometer 
and overspread the linear region of the characteristic curves. Exposures were 
also too lengthy for reliable maintenance of exact conditions in the discharge 
which generates the afterglow. Photography was therefore abandoned for an 
I1-stage EMI photomultiplier with linear galvanometer sensitive to 3 x 10-! amp. 
A camera case and second slit admitted light to the multiplier housing, arranged 
to allow a whole line at once to be integrated but free from adjacent lines. When 
calibrating with the tungsten strip and Planck curve to obtain absolute intensities 
in the Balmer lines, a dispersion curve for the spectrograph was applied to allow 
for the varying dv at each line. Selective absorption in the prisms, etc. is allowed 
for by identical transmission in the optical train for the hydrogen light and the 
light of the tungsten calibration. 

4. Fig. 2 exhibits a mean decrement curve for the afterglow in pure 
hydrogen, averaged over all sets in which constant conditions had been maintained 
throughout the traverse of the whole spectrum by a high-grade micrometer 
screw. Since errors in calibration to absolute intensity are greatest in the red 
with the red weakness of the SbCs photoelectric surface in the multiplier, while 
the absolute Balmer intensities are weakest towards ultra-violet, comparison with 
theories is most readily exhibited by the usual device of reducing all to a common 
value at Hf where the errors are smallest. 

We plot also the intensities in Cillié’s theory of electron capture (19) and 
in the one of Menzel’s theories which uses assumptions the most divergent from 
that of Cillié (5); these are at 5 000 deg., recollecting that temperature depen- 
dence is here slight. For thermal excitation both high and low temperatures 
are plotted, since the dependence is steeper and in the opposite direction; for 
this, Einstein coefficients A,,,,, statistical weights g,, and a Boltzmann factor 
exp (E,,,kT) are employed in the usual manner, where E£,, are excitation potentials 
of the Balmer levels. 

5. Although experience seems never to have shown either astronomical or 
laboratory decrement to obey a single law throughout the spectrum, it is of some 
significance that, even when so many of the complexities of an experimental gas 
are eliminated by this afterglow technique, nevertheless agreement with any 
simple theory only here appears in one portion of the spectrum at expense of 
sharp divergence in another portion. It may be useful to consider whether this 
implies unlikelihood of any undistorted law of decrement in stellar gases. 

Fig. 2 indicates that, if Ha/Hf alone were employed as data, approximation 
either to a high-temperature thermal distribution or to a capture law could be 
reasonably claimed by allowing for the fact that the red is the least accurate region. 
But the excess progressive fading towards violet invites physical discussion, since 
when it occurs commonly in ordinary discharges it is readily explainable either 
as (a) instrumental absorption, or (b) low temperature thermal distribution of 
electron impacts. Instrumental absorption is here eliminated by the black- 
body calibration being carried out with identical optical train. A very low 
temperature thermal excitation would roughly fit our violet experimental curve, 
but would demand an Hz so extremely high that the actual red intensity could 
only be possible through extreme self-absorption in Hz. A subsidiary experiment 
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on the contours of Hx and Hf was therefore carried out, in the full discharge 
light since the afterglow is too faint for contouring, and the amount of self- 
absorption in this geometry of plasma was found far too slight to lower sufficiently 
a concealed intensity of Hx. On the other hand, the high temperature thermal 
distribution which would provide the observed Ha/H ratio would fail for the 
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higher lines as radically as the capture law. It may be noticed that Boltzmann 
distribution of excitation by electron impact, with the artificially lowered E,, 
which would be one of the possibilities if Lyman self-absorption had overpopulated 
the second quantum level, would imitate a combination of very high temperature 
with normal E,, from ground level, reintroducing the violet discrepancy. 

We may classify Balmer decrements, astronomical or experimental, as (1) 
decided by pure electron capture at recombination, distribution of the states being 
as in Cillié or Menzel, (II) excitation from ground level distributed thermally, 
or (III) distortion of either law by overpopulation or underpopulation of certain 
Balmer levels. The above description of Fig. 2 therefore indicates a need to 
explore (III), which will also be required in astronomical instances. The follow- 
ing suggest themselves. 

(a) Reduction of higher Balmer levels by superelastic “ collisions of the second 
kind” which transfer energy, instead of by radiation, to either (i) molecules 
capable of dissociation or atoms accepting discrete quanta plus kinetic energy, 
or (ii) electrons or other free particles accepting kinetic energy. 

(6) Overpopulation of low levels by self-absorption of Lyman lines enabling 
either (i) slow electrons or (ii) low frequency radiation to excite adjacent levels to 
excess compared with capture or thermal excitation from ground level. We have 
mentioned a difficulty for (i) if it were to imitate a very hot gas. 

(c) Other effects of Lyman or Balmer self-absorption, and the incomplete 
cycles distorting the distribution, as applied in astronomy by Woolley and others. 
Experimental test of some of the possibilities was carried out as follows. 

6. It is probable that an afterglow is primarily a recombination spectrum, 
however distorted, and this is reinforced in our case by the cathode ray pictures 
which trace the luminosity up to a millisecond after voltage fall when little electron 
velocity for excitation can survive. This probability makes it worth discriminating 
between the above (III) (a) (i) and (ii), either of which would weaken progressively 
' the radiation from upper levels, as the observed data indicate, since collision 
cross-section of the higher states will grow rapidly with the quantum number. 

Such radiationless de-excitation is familiar in the chemistry of sensitized 
fluorescence and the quenching of resonance radiation, and all phenomena where, 
instead of radiating a line, an excited state loses its energy to a colliding atom 
with appropriate energy levels or to a free particle capable of acquiring velocity. 
There is in pure hydrogen itself the possibility of levels available for such transfer, 
as would steepen the Balmer decrement as observed; for instance, there is the 
known loss at about I1 volts yielding dissociation into fast atoms. Repetition of 
the experiments with hydrogen—-neon and hydrogen-helium mixtures admitted 
as in Fig. 1 enables a test to be made, since neon also has an approximately I1 volt 
level whereas helium possesses no level below 21 volts and so is incapable of 
undergoing quantum transitions by de-excitation of hydrogen from Balmer states. 
These two gases were selected because their own spectral lines leave Hf and Hy 
fairly isolated until higher ionization, whereas argon or nitrogen lines and bands 
would make impossible the spectrophotometry of Balmer lines. 

Since ionized inert gases add to the slow adsorptive and desorptive reaction 
of hydrogen with glass walls, intensities of lines must be compared with great 
rapidity; the method adopted was to await stabilization of a discharge in which 
the inert spectrum was at least as strong as the Balmer, and then to traverse 
rapidly from H8 to Hy and back repeatedly until the intensity ratio of the two 
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lines was determinate even if absolute values tended to wander. ‘The repeatable 
results are shown in Table I, expressed as absolute value of Hy referred to an 
H£ of 1000 units. It is noticeable at once that any radiationless depopulation of 
upper Balmer levels cannot be invoked as completely explainable in terms of 
those superelastic collisions which transfer the energy to discrete levels; because 
helium, rendered incapable by its 21 volt minimum, is at least as effective as neon 
in steepening the decrement when they replace hydrogen in the mixture. The 
intensity of Hy does, however, fall decisively with both inert diluents, and a feature 
common to both would be that electrons or other free particles accepting kinetic 
energy may either de-excite the hydrogen levels or assist hydrogen particles 
themselves to do so. It is possible that the inert gases merely provide additional 
electron pressure, a suggestion gaining a little support from the fact that in 
pure hydrogen itself at fixed gas pressure a lowered electron concentration raises 
Hy. Again, when by mismatching the circuit a high current in the oscillator 
fails to produce the usual predominance of Balmer over molecular spectrum, this 
increased Hy characterizing the low luminosity tends to drift back to its lower 
values, as if molecular ions and their electrons play a somewhat similar part to 
that of the electrons contributed by the inert gases. Loss of energy to free 
electrons instead of into quantum states was, of course, the original “superelastic 
collision’’ of the early theories, but has tended to be obscured by the chemical 
phenomena: that it may play a part in steepening the Balmer decrement of the 
denser among astronomical spectra is a first tentative suggestion from these 
experiments. 


TABLE I 
Intensities of Hy when HB=1 000 
Mean of all previous, high current in pure H 
Mean of new set at single pressure, high current pure H 
Mean at same pressure, low current pure H 


Mean of hydrogen in neon 
Mean of hydrogen in helium 


7. Since it is impossible to decide whether this qualitative account is adequate 
quantitatively, as the cross-section for Balmer states to become de-excited by 
collision is unknown, it will be of some importance to mention briefly a different 
explanation, linking astronomy and experiment, which must also to some degree 
contribute to the observed steepening of the decrement. We turn to III (6) 
and (c), processes depending upon radiation absorption and verging upon the 
definitions of fluorescence in astronomical spectra. Here it is worth recollecting 
that Woolley’s theory (7) of uncompensated Lyman absorptions, based on 
Rosseland’s cycles, alters a Balmer decrement from H«:H8:Hy:: 4:10: 1-00: 
0-44 (Cillié capture) to 4-4: 1-00: 0-34 for exposure to weak Lyman radiation and 
9-3: 1°00: 0°30 for exposure to greater Lyman strength. These are in Schwarz- 
schild’s units of intensity, which somewhat exaggerate H« but leave the other 
lines almost as on our absolute scale. 

Woolley’s theory neglects collisions; it rendered more comprehensible the 
Balmer decrements in collisionless regions, i.e. nebulae and prominences, and it 
offers an intermediate type between pure capture ignoring absorption and the 
excitations of thermal equilibrium; it has not been considered for denser gases 
subject to collisions. Our discharge is not likely to obey even the polychromatic 
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radiative equilibrium of Woolley, but its Lyman absorptions will distort the 
recombination law, as we listed in III (c), in the direction of Woolley’s figures, 
as well as permitting unbalanced radiative excitations as in III (6). Lacking the 
simplified conditions assumed for Woolley’s gases, no quantitative estimate can 
be made; but if we add to the collisional de-excitation of upper levels suggested 
by our experiment (Section 6) any consequences of Lyman absorptions in the 
direction indicated by Woolley, the two effects would seem to satisfy the 
peculiarities of both ends of our experimental decrement graph, better than 
other known processes. For instance, besides the ultra-violet fading, the excess 
in Hx becomes reasonable without the forced explanation of a thermal distribution 
which would carry the intensity of that line far higher than the experimental data. 

8. In tentatively introducing Lyman absorption cycles into experiment when 
they were first discussed in astronomy, the other experimental possibility, of 
de-excitation by collision, which we suggested, may in turn have some astronomical 
relevance: a balance between the two effects may be drawn as follows. 

(i) Among astronomical emission spectra, the solar chromosphere and 
prominences, nebulae, and those Be stars whose lines denote genuine expanded 
chromospheres, will tend to be controlled by radiative absorption more than by 
collision, and decrements will be distorted more according to the principles of 
Section 7 than of Section 6. 

(ii) With emission in Me stars, where Joy discovered Balmer bright lines 
below metallic bands, and in those Be where absorption also appears to cut into 
low-lying emission, the Balmer decrement may well denote gas pressures greatly 
exceeding the nebular or chromospheric, and collision as in our experiment may 
no longer be negligible. In particular, where McKellar and Odgers (20) re- 
cently locate the Balmer of Me stars in the rising of a photospheric convection 
zone, the conditions approach our type of experiment. At all these astronomical 
situations of higher pressure, where emission dominates at Ha fading to absorption 
components in the violet, such steep decrement indicates a physical counterpart 
of the ionic assemblage upon which we experiment; it will be worth considering 
whether there the consequences of fluorescence of Woolley’s and Rosseland’s 
type are likely to be reinforced or even superseded by de-excitation at collision, 

(iii) In erupting stars and the blue-red “‘ binaries’, the Balmer “ progression” 
of Merrill contains the important features of decrement varying with depth of 
material, a difference in radial velocity for different members of the Balmer 
series; when the true decrement is extracted from the effects of blending, a 
distinction between our two-fold experimental suggestions of fluorescence and 
of de-excitation could assist the essential but difficult separation between the 
chromospheric and the photospheric locations of the lines. The disappointing 
insensitivity of capture laws to temperature may be replaced by a useful sensitivity 
of decrement to foreign electron pressure, just as the opacity of cooler stars in 
the main sequence has been valuably attributed to negative H~ formed according 
to the excess electrons from metallic ionization. 

I am grateful to Professor P. B. Moon, F.R.S., for facilities and encouragement, 
to Professor James Sayers of the Electron Physics Department for discussion 
and criticism, and to Mr John Griffiths, technical assistant in astrophysics, for 
construction of apparatus and for sharing in all operations and measurement. 

Biv. gnam University : 

1 3 March 18. 
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ON THE PROBLEM OF THE H« EMISSION IN 
THE SHELL STARS* 


Anne B. Underhill 
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Summary 


Some observations of the emission features at Ha and Hf in the spectra 
of the Be stars 7 Tauri, 8 Canis Minoris, the peculiar Be star « Draconis, 
and the shell stars ¢’'Tauri and 48 Librae are presented. It is shown that 
the widths of the emission lines of the Be stars are consistent with the 
hypothesis, postulated some years ago by Struve, that the emission lines of 
Be stars originate in a rotating equatorial bulge surrounding the rapidly 
rotating star. However, this hypothesis is not adequate to explain all the 
emission features of the spectra of shell stars where Ha is shown to have 
very extensive wings. It is proposed that the excessive broadening of the 
emission at Ha in shell spectra comes from material in tenuous streamers 
at comparatively large distances from the rapidly rotating star. These 
streamers are thought to be moving at great speeds in random directions and 
are only visible in Ha radiation. The equatorial bulge and inner shell of 
the shell star are thought to give rise to the ordinary shell spectrum of 
absorption and emission lines. 


It has been known for many years that generally some of the Balmer lines of 
hydrogen are in emission in the spectra of the B-type stars that show shell spectra. 
Usually there is a strong emission line at Hx and weaker emission at H8 and 
sometimes at Hy. Since it would appear that with regard to the emission lines 
of hydrogen the shell stars are related to the Be stars, some quantitative measures 
of the emission features at H« and Hf have been made in order to investigate this 
relationship. Significant differences between the shell stars and Be stars are 
found. 

Investigations by Struve and others (1), based chiefly on observations of the 
spectra of the Be stars in the blue-violett, have shown that the hypothesis, that 
the emission lines originate in a rotating equatorial bulge or ring surrounding the 
rapidly rotating main-sequence B-type star, is valid for these stars. Important 
points in favour of this hypothesis are the symmetrical, double nature of the emission 
lines and the observation that the total width of the emission features expressed 
in velocity units is the same for all emission lines in the spectrum of any one star. 
The velocity of rotation of the region giving rise to the emission may be estimated 
from the total width of the emission features. The value derived is consistent 
with the hypothesis that the velocity of rotation in the equatorial extension 
decreases outwards from the star in such a way that angular momentum is 
conserved. Such behaviour is to be expected for material at the estimated 

* Contributions from the Dominion Astrophysical Observatory, No. 32. 
+ In 1906 R. H. Curtiss began investigating the Be stars, and his observations, particularly those 


in Publ. Obs. Univ. Mich., 3, 1, 1920, form the basis of the theories of the Be stars, proposed by 
Struve et al. 
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temperatures and densities of these emitting rings, for both radiative viscosity 
(2, 3) and eddy viscosity (4) would be very low, and therefore the gravitational 
field will dominate the motion. As regards shell stars, the rotation hypothesis is 
consistent with the observed shapes and widths of the shell absorption and 
emission lines that occur in the blue-violet spectral region, and this hypothesis 
has been used for example by Struve (§) and by Underhill (6) to obtain an estimate 
of the dimensions of the shells surrounding certain B-type stars. 

According to the simple rotational hypothesis the width of the emission at 
H« expressed in velocity units should be the same as that at Hf, both for the 
Be stars and for the shell stars. ‘To test this hypothesis, high-dispersion spectra 
of some Be stars and some shell stars were obtained in the spectral regions of 
Hf and Ha. For most of the plates the new universal spectrograph in its three- 
prism form with the M camera (focal length 71 cm) was used at the Cassegrainian 
focus of the 73-inch telescope. The linear dispersion of the spectra is 13-8 A/mm 
at Hf and 45-6A/mm at Hx. Grating spectra of the Be star » Tauri obtained 
by K. O. Wright were also available for study. ‘The dispersion on these plates 
is 7-2 A/mm at Hf and 6-8 A/mm at Hx. Each plate was calibrated for intensity 
work by means of the rapidly rotating sector and auxiliary spectrograph regularly 
used at this Observatory for this purpose. Eastman IID plates were used for 
the Hf region and IIF plates for the Hx region. The plates were developed 
in Dig. The most intense parts of the emission features were not over-exposed 
on any of the plates. 


"TABLE I 


Spectral Mean Date Mean Date 
of of 


Type of Observation eee of Observation Plates 


23630 | » Tau Bse 1949 Jan. 4 1948 Nov. 
37202 | ¢ Tau Bap 1953 Mar. 6 1953 Mar. 
58715 | BC Mi B8e 1953 Mar. 26 1953 Mar. 
109387 | x Dra Bse 1953 Feb. 28 1953 Feb. 
142983 | 48 Lib B3p 1950 June 19 1951 June 24 


Microphotometer tracings magnified 100 times in the direction of the 
dispersion were made of the regions near Hf and Hz, and these were reduced 
to direct intensity tracings with the semi-automatic intensitometer. ‘The 
direct intensity tracings for each star were all reduced to the same scale so they 
could be superimposed and a mean profile obtained. ‘The stars studied, the 
number of plates used, and the mean date of observation are given in Table I. 
The mean profiles of Hx and the central portion of Hf for the Be stars » Tauri 
and £ Canis Minoris are given in Fig. 1. ‘Those for the peculiar Be star « Draconis 
are given in Fig. 2, while those for the shell stars ¢'Tauri and 48 Librae are given 
in Figs. 3 and 4. The ordinate is intensity in terms of the continuous spectrum 
at the line, and the abscissa gives the displacement from the centre of the line in 
km/s. Each plotted point gives the mean intensity at that wave-length. ‘These 
profiles are true profiles since the instrumental broadening is negligible for lines 
of this width. 
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It is clear that in the spectra of the Be stars 7 Tauri and 8 Canis Minoris the 
emission lines are pretty well symmetrical and that the width of the emission at 
Hz in velocity units is the same as that at Hf, in agreement with the hypothesis 
that the emission lines originate in a rotating equatorial bulge or ring. Since the 
stellar hydrogen absorption lines are much stronger in a B8 star than in a B5 
star, the dip of the continuum due to these lines is greater in 8 Canis Minoris 
than in 7 Tauri. Although only small portions of the profiles of the stellar lines 
are shown, it is evident from the shape of these lines that the underlying stars 
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Fic. 1.—Profiles of the emission features at Hx and H® in the spectra of the Be stars, 
B Canis Minoris and » Tauri. The ordinate gives intensity in terms of the continuous spectrum. 
The continuous spectrum is indicated by a straight line at intensity 100. The abscissa gives 
displacement from the centre of the line expressed in km/s. At Hx one angstrom=45°68 km/s, 
while at HB one angstrom= 61°67 km/s. 


are rotating rapidly. It should be recalled that all the observed profiles are 
composite, consisting of a shallow dish-shaped absorption line from the rapidly 
rotating underlying star and an emission feature. It is clear that the true widths 
of the emission features may be greater than the apparent widths, since it is nearly 
impossible to delimit accurately the filling in of the shallow absorption feature 
at the edges of the emission. In most cases, at H« all trace of the absorption 
feature is lost, and the emission appears to rise from the continuum, which is 
drawn as a straight line at intensity 100. 
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In the case of the peculiar Be star x Draconis, the emission feature at Hx is 
somewhat wider than that at Hf, and it does not show the typical double peak of 
a Be star, although the emission at Hf does. The portion of the stellar line 
profile shown at Hf indicates that the underlying star is rotating rapidly. Both 
the radial velocity and the intensity of the emission features of « Draconis are 
known to vary in a somewhat irregular manner (7). Slettebak (8) has estimated 
that the projected velocity of rotation of « Draconis is 280km/s. Since the wings 
of the emission feature at Ha extend to about + 400km/s, it would seem that this 
emission feature is broadened by some cause in addition to rotation. 

The profile of the emission feature at Ha in the spectrum of the shell star 
¢'Tauri has very extensive wings, which, expressed in velocity units, are much 
wider than the wings of the emission at Hf. It has been pointed out previously 
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Fic. 2.—Profiles of the emission features at Ha and HB for the peculiar Be star, « Draconis. 
The coordinates are the same as for Fig. 1. 


by the author (6) that the width of the emission at Hf and the detailed appearance 
of the shell spectrum in the blue-violet are consistent with the simple hypothesis 
that the absorption and emission features of the shell spectrum originate in an 
equatorial bulge of the rapidly rotating star, and that the material in this bulge 
rotates so that angular momentum is conserved. However, this explanation is 
quite inadequate to account for the extensive emission wings of Hx. The under- 
lying star is rotating with a projected velocity of about 350 km/s (6), but even this 
velocity is too small to account for the emission wings of H« which extend to 
about + 800km/s. 

The emission feature at Ha in the spectrum of the shell star 48 Librae is also 
much wider in velocity units than is the emission at Hf, although the wings do 
not appear to be quite so extensive as those for Tauri. (The observations for 
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45 Librae are not quite so good as those for ¢ Tauri, and it is possible that the 
extent of the emission wings at H« and at Hf have been underestimated in the 
spectrum of 48 Librae.) Slettebak (8) has estimated that the underlying star of 
45 Librae rotates with a projected velocity of about 400km/s. It is evident that 
for this shell star also the simple rotation hypothesis that was adequate to explain 
all the emission features of the Be stars breaks down. ‘The width of the emission 
at Hf is of the order of +250km/s, and this width could be explained on the 
rotation hypothesis in the same way as for ¢ Tauri. 


-1000 


-800 


-600 


-400 -200 ° #200 +400 +800 #1000 «M/SEC 
qT T 


¢ TAURI 


“1000-800 -600 -400 -20C © 4400 +600 +800) +1000 KM/SEC 


Fic. 3.—Profiles of the emission features at Hx and HB for the shell star € Tauri. 
The coordinates are the same as for Fig. 1. 


The profiles of the emission features at Hx and Hf have been given here for 
only two shell stars. However, it may be seen from the pictures and intensity 
tracings of the spectra of shell stars published by P. W. Merrill (g) that in general 
for the shell stars the emission at Ha is very much wider in velocity units than 
that at Hf. It appears for the B-type shell stars that the emission feature at H« 
extends to at least +650km/s, while that at Hf extends to about +250km/s. 
The underlying stars rotate with projected velocities of the order of 350 km/s, hence 
it is possible to explain the width of the Hf emission on the rotation hypothesis, 
which is successful for the ordinary Be stars, and which accounts for the appearance 
of the shell spectrum in the blue-violet spectral region, but it is not possible in 
this way to account for the great width of the emission feature at Hz. Some 
other cause of broadening must be operative for Hx. ‘The observed width of 
the emission is not caused by photographic errors or by instrumental broadening 
by the spectrograph. 

Since, on the whole, the emission features are symmetrical and the absorption 
component, if it exists at all, is more or less centrally placed, there can be little 
likelihood of a uniform outward (or inward) motion of the emitting material of 
the order of several hundred km/s, such as is postulated for the PCygni or 
Wolf-Rayet stars. It is possible that the extra broadening of the emission at 
Hz« is caused by large chaotic motions of the emitting material. ‘These motions 
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would have a root-mean-square velocity of the order of 100 km/s. It is significant 
that the excessive broadening of the emission feature is only shown by the leading 
member of the Balmer series and not by the second member Hf. It is suggested 
that the extra width of the emission at H« for the shell stars comes from tenuous 
streamers of gas moving with high speeds and at relatively large distances from 
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Fic. 4.—Profiles of the emission features at Hx and HB for the shell star 48 Librae. The 
coordinates are the same as for Fig. 1. The fluctuations of the photographic grain have been 
smoothed for the profile of HB, but not for any of the other profiles. 


the rapidly rotating star, which has a fairly dense equatorial bulge giving rise to 
the regular absorption and emission features of the shell spectrum. The extra 
emission is not recorded at Hf because the density of the material in the streamers 
is so low that, when account is taken of the Balmer decrement, there is not enough 
light emitted at Hf to record on a photographic plate. 
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Observations of the K-type giant star 31 Cygni (10) have indicated that large 
irregular motions do occur in the tenuous outer atmosphere of a late-type giant 
star, hence it is not difficult to imagine similar motions occurring in the tenuous 
outer regions of a shell star. ‘The difference between shell stars and Be stars is 
probably one of degree, the Be star not having sufficient material in its equatorial 
bulge and surrounding tenuous atmosphere to produce the typical shell absorption 
and emission features. The peculiar star x Draconis is probably intermediate 
between a Be star and a fully developed shell star. 


Dominion Astrophysical Observatory, 
Victoria, B.C. : 
1953 April 21. 
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SOME CALCULATIONS INVOLVING 
THE RESTRICTED THREE BODY PROBLEM 


K. N. Dodd 


(Communicated by W. H. McCrea) 


(Received 1953 May 7)* 


Summary 


The results are given of some calculations made on the Manchester 
University Electronic Computer concerning the angular momentum exchange 
between the bodies in the classical ‘‘ restricted three body problem”. From 
these calculations, a numerical estimate is made of the torque experienced 
in a particular case by a binary star in moving through a cloud of interstellar 
material. 


1. Introduction.—If a binary star is in motion through an interstellar medium, 
there is in general an interchange of angular momentum between the binary and 
the medium, with consequent effects upon the orbital elements of the binary. 
A full mathematical treatment of the hydrodynamical problem is quite beyond 
the range of existing techniques. In order to estimate the effects, use has been 
made of known results for simpler problems whose adaptation has to be justified 
by certain physical arguments (1). In this simplified treatment, the interstellar 
material is supposed to consist of particles which move under the gravitational 
influence of the components of the binary but do not influence one another. 
This supposition is equivalent to assuming that the interstellar material is a gas 
at zero absolute temperature. On this supposition, if the orbits of the binary 
components are circular, the determination of the motion of a given particle 
of interstellar material is reduced to the classical “ restricted three body problem ”’. 
A considerable amount of computation on this problem was performed by 
E. Strémgren (2) and his collaborators between 1919 and 1931 in order to deter- 
mine the forms of periodic orbits of the bodies. These calculations were almost 
entirely concerned with the motion of the bodies in a plane. 

The orbits (for various initial conditions) of a particle of negligible mass in 
the field of a binary star have recently been determined by numerical integration 
on the Manchester University Electronic Computer. ‘The binary was taken to 
consist of a pair of particles of equal mass, pursuing a common circular orbit. 
The present paper is an account of this work and its application. In Section 2, 
the numerical method which was used to determine an orbit is described. 
Section 3 is a statement of the specific problem to be solved and of the orbits 
that were determined. In Section 4, the results are used to estimate the torque 
experienced by a binary in a particular case of its passage through a cloud of 
interstellar material. In Section 5, this torque is used to determine the resulting 
rate of reduction of the separation of the binary components. 

It should be stated explicitly here that the particle-orbits are not restricted 
to plane motion but are in three dimensions. 


* Received in original form 1952 October 30. 


No. 4, 1953 Calculations involving the restricted three body problem 485 


2. Numerical method.—We wish to determine the motion of a particle C of 
negligible mass in the gravitational field of two other particles A and B which 
are of equal mass and are revolving under their mutual gravitational attraction 
in acommon circular orbit of radius a, the centre of which is at rest in a Newtonian 
frame of reference. The two particles A and B are always at the ends of a 
diameter of the common orbit and each moves with the velocity 


V=(u/4a)™, 


where ». = mass of A (or B) x constant of gravitation. 

If we place a frame of reference Oxyz with its origin O at the centre of the 
circular orbit and such that this orbit is in the x, y plane, the equations of motion 
for the particle C are 


x=u, y=v, z=w, 0=Vt/a, 


where (x, y,2) are the coordinates and (u,v, w) the velocity components of C at 
time ¢ and 


R,? =(x—acos 6)? + (y—asin 6)? + 2’, 
R,? =(x+acos 6)? + (y+asin 6)? + 2%. 


In the calculations by Strémgren the bodies A and B were reduced to rest by 
a suitable change of variables but this would be of no help for our purpose. 
So the equations of motion were integrated as they stand, except for a transforma- 
tion to make the variables dimensionless. This was obtained by measuring 
distances in terms of a, velocities in terms of (./a)" and time in terms of (a3/)"?. 
The variables were also scaled down so as to be accommodated in the number 
range of the machine in the manner customary in this type of computation. 

The method of integration used was a modified form of the Runga-Kutta 
process which has been developed by S. Gill (3) especially for use in electronic 
machines. We used this process in the following form: Let the independent 
variable be yo and the dependent variables be yy), ({=1,2,...,”) and the 
differential equations be 


dy; = 


Let Y be the initial value of y) so that y, 9=y9,(Y) are given. Let A be the inte- 
gration step in yo, then we calculate the following quantities with « =o, 


Ria =hfyo, Vn, 

— (Ag+ 2Ba)Ry as 


| 
u= acos 6) 
v= Ri” asin @) 
w= RS” 
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where q,9=0 and A,, B, are given in the following table: 


A, B, 
3 
I 
2 I+} 
3 


From the values for « =0, the values for « =1,2 and 3 are calculated. ‘Then 


+h). 


The process is repeated, starting with y(Y+h) to obtain y(Y+2h) but the 
values taken for g; will be the same as g;, of the previous step. The error 
involved in this process is of order A®. To obtain the values of the variables at 
a given point, the method requires a knowledge of the variables only at the 
previous point and not at several previous points. ‘This removes the need for 
storing previous values in the machine and also the need for shifting these values 
after each step. Another advantage is that no special starting procedure is 
necessary. The quantities g,, are used to give an indication of the previous 
behaviour of the variables but they are usually small and there is practically no 
loss of accuracy if they are taken to be zero at any point. ‘This enables changes in 
the magnitude of h to be made without the need for any special calculations. 
Although the calculations described here were performed with a fixed integration 
step h, it is possible to make the machine automatically adjust the size of the step 
so as to keep a specified accuracy. It is also obvious that this integration method 
on the machine can be used for the general problem of three bodies or in fact 
of any number of bodies subject to the storage limitations in the machine, but the 
time required for each integration step would be increased. 

The cosines and sines in the equations of motion were calculated directly 
from the power series using terms up to that involving ©!*, where © is such that 


6=nn+© (n being any integer) 
and 0<O<z/2. 


The machine took 11 seconds to advance the integration by one step. This, 
however, involved two actual calculations and a comparison of the results to 
eliminate the risk of ‘“‘random”’ machine errors. At intervals of two or four 
steps, the values of x, y, 2, u, v, w, t and (x? + y?+ 2") were printed. Each printing 
required 10 seconds and involved six decimal digits and a sign for each variable. 

If the particle C were to approach too closely to either A or B, the rapid 
variation in the gravitational field would considerably affect the accuracy of the 
solutions and this might not be noticed by examining the printed results. To 
remove this risk, it was arranged for the machine to stop if R, or R, became less 
than 4a. This did not occur very often, so that few orbits were lost on this 
account. 

3. The problem.—The main problem consisted of finding the angular 
momentum gained by the particle C when it is projected parallel to the z-axis 
towards the revolving bodies A and B. More precisely, the initial conditions 
were 


x=Rcos¢, y=Rsin¢d, 
u=0, v=0, w=—W, 
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where various values of R and ¢ were taken; ¢ being measured from a direction 
parallel to the x-axis. Orbits were obtained for one value of W only. This was 
1-6(u/a)"?. This value of the velocity is of the right order to give a result for 
comparison with our other methods of estimation. The particular factor 1-6 
was selected from the relevant range because, after the variables have been 
scaled down, it gives a convenient number in the scale of two, for use in the 
machine. The orbits were computed in most cases until the z coordinate of 
C was less than —4a. The initial and final values of z (i.e. + 4a) were arbitrary 
but were taken to be sufficiently large, so that beyond these limits the angular 
momentum of C would not have been appreciably changed. 

For each of the orbits obtained, the final value of the angular momentum per 
unit mass of C about O (i.e. xv—yu) was then calculated on a desk machine. 
The initial conditions were, of course, such that the initial angular momentum 
about O was zero in all orbits. 


TABLE 
(xv 
0-754 a 1'25a 2a 3a 4a 
+0:3000 +0°5170 +0°6308 +0°6265 +0°5518 +0°3594 +0°1260 +0°0481 


+0-7710 +1:0000 +1°0070 +0°8742 +0°5504 
+0°4577. +1°0187. +1°4981 +1°5183 +1°2370 +0°6902 +0°2232 +0°0876 


+1:1260 +1°8744 +2°4128 +1°5078 +0°6712 

+0°2458 +0°7457 [+1°90] [+290] +0°8266 +0°3272 +0°0978 +0°0410 
—o0'0502 [—o'25] —o'8070 —o-2821 

—o'1382 —o-4189 —1°'3243 —1°4007 —0°6559 —0'1536 —o-0502 
—0°4642  —0°7987  —1°0435 —1°0385 —0°6574 

—0'1674 —0°3426 —0'5140 —0'6193 —0'6270 —0°0845 
—0'2278 —0'2776 —0°2653 —0°2673 

+0'0404 +0°0549 +0°0461 +0°0203 —0'0099 —0'0532 —0'0579 —0°0338 
+0°2784 +0°3241 +0°3118 +0°2647 +0°1556 


+0°1230 +0°2567 +0°38 +0°40 +0°'0928 +0°0297 +0'0058 +0-'0014 


In the accompanying table the angular momentum per unit mass gained 
by C is given in the dimensionless form. To obtain these orbits, a step of 0°25 
in dimensionless time was used. About 20 steps were required for each orbit. 
Most of the values given are correct to within one per cent. This limit of 
accuracy was sufficient for our requirements but the machine could be made to 
work to a considerably higher accuracy. Orbits were not obtained for the four 
cases: R=a, 1°25a and ¢=60°, 75°, owing to the close approach of C to either 
A or B resulting in the machine being stopped as mentioned above. The values 
given in the table (in square brackets) are obtained by graphical interpolation 
between the numbers in the columns. The values in the range ¢ = 180° to 360° 
are a repetition of those given in the table. Inspection of the numbers in any 
column shows that they lie approximately on a sine curve and the oscillation is 
about a mean value that is small compared with the amplitude of the oscillation. 
This situation makes it difficult to obtain an accurate value for the mean. Since 
the function is necessarily periodic (period 7) the best method (5) is simply to 


R 
$ 
30° 
45° 
60° 
75° 
go” 
105° 
120° 
150° 
165° 
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take the arithmetic mean of the values in each column. The mean values, f(R), 
are given in the bottom row of the table. ‘These values are accurate to about 
5 per cent except for the cases R=a and 1-:25a where only 10 per cent accuracy 
is claimed. 

4. Calculation of the torque experienced by a binary.—From the results of the 
work described in the last section, it is possible to calculate the torque experienced 
by a binary moving in a certain way through a cloud of interstellar material. 
For consider the binary consisting of A and B to be moving in the direction Oz 
with velocity W through a cloud of particles which are at rest at a great distance 
from the binary. By giving the entire system a velocity — W, we have the situa- 
tion described in Section 3, neglecting the effect of the gravitational field upon 
the angular momentum outside z= + 4a. The mass of material which crosses 
the plane z = 4a between the circles 

x?+y?=R*? and x7+y?=(R+dR)* 
in unit time is 
27RdR. Wp, 
where p is the mass of particles per unit volume in the undisturbed cloud. ‘The 
angular momentum gained by this mass in its encounter with the revolving 
system is 
27RdR. Wof(R), 


so the total angular momentum gained by the material in unit time is 


2nWp | RA(R) AR. 


This is the rate of loss of angular momentum by the rotating system and so is 
equal to the torque acting on it. ‘This expression was evaluated for W = 1-6(p/a)!* 
by plotting Rf(R) and obtaining the integral by counting squares. The result 
was 

2mpa*u x 0-6. (1) 


The function Rf(R) rises and falls steeply near its maximum value. Conse- 
quently, the value of the integral is not appreciably affected by an inaccuracy in 
the maximum value of the function. ‘The maximum value of Rf(R) was estimated 
from the graph to be 05 at R=1-25a. A 20 per cent error in this maximum value 
would only have caused a 5 per cent error in the value of the integral, but as we 
have stated in Section 3, we do not expect this maximum value to be in error by 
more than 10 per cent. [As an alternative to the graphical method, an attempt 
was made to fit a curve of the type Ax*exp(— Bx). This was fitted to two of 
the accurately known values; namely, at R=0-75a and R=1-5a. The values 
of A and B were 1°89 and 2:27 respectively. ‘The result of the integration was 
then 13 per cent less than that obtained by graphical integration, but as the curve 
did not fit very well, the result from the graphical method is considered to be the 
more accurate.] After this discussion of various errors we consider that 10 per 
cent accuracy can be claimed for the final result (1). 

A question arises as to whether the effect of particles passing very near the 
revolving bodies is large compared with the effect of the remaining particles. 
In the immediate neighbourhood of one of the bodies, a particle will behave 
almost as if it were influenced by this body alone. The effect of the particles 
on a single star has been examined in (4) and equation (3.6) of that paper gives 
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the force on the star due to the particles passing within a given distance of the star. 
It is seen that this force tends to zero as the given distance tends to zero. It 
follows that particles, passing very close to either component of the revolving 
system considered here, do not contribute appreciably to the total effect. 

It will be noticed that we have neglected any change in the motion of A and 
B. This is justified because, as we shall see in Section 5, the orbital elements 
of the binary are only appreciably atfected by the torque after a period which is 
several times the period of revolution of the binary. 

4.1. Comparison of the result with a theoretical formula.—It is possible to 
obtain a theoretical formula for the torque in cases where W/V is large compared 
with unity. From equation (3.6) of (4), the force exerted on a single star due 
to its passage through a cloud of interstellar material is 

In (1 + (2) 
where p is the density of the undisturbed cloud, U is the velocity of the star 
relative to the cloud and » = mass of star x constant of gravitation. & is a “cut- 
off distance ”’ which is of the erder of half the distance to the nearest neighbouring 
star. ‘The force acts in the opposite direction to the velocity. If we assume that 
this formula holds for each component of a binary star which is moving perpendi- 
cular to the plane of its orbit, we put 

X=a and U?=W?+V? 
in (2) to obtain the force on each component. ‘The force on each component can 
be resolved into a force parallel to the direction of motion of the binary relative 
to the cloud and into a force tangential to the common circular orbit. The forces 


of the latter type, acting on the two components, together form a torque of value 
T = In (1 +a? W4y-*) (3) 
if W/V is large compared with unity. This condition is imposed by the 
mechanism which produces the force (2). Each star “tunnels” through the 
interstellar material and this condition ensures that one star does not tunnel 
through a region which has already been disturbed by the other star. 
Putting W =1-6(u/a)!? we obtain 


T = 27pa*u x 0°49, 


so that the formula is in satisfactory agreement with the numerical result (1). 
In this case, W=3-2V, so this agreement is in the region of W where we might 
hardly expect the formula to hold accurately. It is also to be noted that the 
formula does not overestimate the torque. 

It should be noted that this agreement between the formula (3) and the 
numerical result is only offered as additional support for the formula. Formula 
(3) rests on the reasoning given above and so its application is not limited to the 
value of W used in the numerical work. 

5. Rate of reduction of separation of binary components.—If we neglect any 
variation in the masses of the components of the binary, the effect of the torque 
is to reduce the separation of these components. To determine the rate of 
reduction of this separation, let E be the total energy of the binary at time ¢, 
ie. E=}mV?+}mV?—}hym?/r=—ym?/4r, where m=mass of A (or B), 
y=constant of gravitation and r=the radius of the orbit of the binary at 
time ¢t. In a time interval dt the loss of energy is 

— dE = ym dr/4r*. 
34 
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But this also equals the work done by the torque T in dt, i.e. 
TV dt/r. 

2T 


Hence 


so that a 10 per cent reduction of the separation will occur in the time 
m(yma)"?/20T, 
approximately. 

We may illustrate this formula by a specific example, using the value of 7 
given by (1). Consider a binary consisting of two equal stars, each of 5 solar 
masses and separated by a distance of o-I parsec. If this binary moves through 
a cloud of interstellar material of density 10-*? g/cm? in a direction perpendicular 
to the plane of its orbit with velocity 1-05 km/s, the separation will be reduced 
by 10 per cent in about 2°68 x 10’ yr. The period of revolution of the binary 
is 9°36 x 10° yr. 

Y/a 


Pic.. x. 


6. Another set of orbits.—A set of less accurate integrations were performed 
to find the orbits of the particle C after being released from rest at various points 
on the sphere 

x? + y? + 2? = (4a). 
For some positions of release, the particle C moved once through the circular 
orbit of the bodies A and B and then receded to a distance greater than 4a from 
the origin. For other positions of release, the particle C receded to a distance 


| 
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of less than 4a from the origin and then returned for one or more encounters with 
the revolving system before escaping from the sphere. One of the more com- 
plicated orbits was computed accurately and projections of it on the planes 
y=o and z=o are shown in the accompanying figures. In one of these, the 
positions of C are marked when AB has turned through various multiples of $7. 
The accuracy of the curves is not guaranteed beyond the point D where the 
particle has a close encounter with one of the revolving bodies. In the orbits 
which were computed the particle C eventually escaped from the sphere and it 


Z/a 


Fic. 2. 


seems certain that this would always occur. For, it is impossible for C to settle 
into a periodic orbit around the revolving system; otherwise, by reversing the 
time coordinate, it would be possible for periodic orbits of C to become non- 
periodic. When C escapes from the sphere of radius 4a, it has gained energy 
from the revolving system. It will also have gained a proportional amount of 
angular momentum. ‘This follows from the fact that the revolving system as a 
whole has no translatory motion and so any forces (as opposed to torques) acting 
on it during the encounter with the particle can do no work. The angular 
momentum per unit mass gained by C (i.e. xv —yu) while inside the sphere was 
evaluated for the various orbits and was considered as a function of the point 
of release of C. Then by a crude integration over the surface of the sphere, a 
mean value was obtained. This was 


This result can be used to obtain the torque exerted on the rotating system when 


material falls steadily from the surface of the sphere and there are no particle 
collisions. We may compare this result with that obtained on the assumption 


34* 
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that the particles collide and are eventually all captured by one or other of the 
revolving bodies. In this case the angular momentum per unit mass gained by 
a particle C is aV (V being the velocity of A or of B = }(y/a)"*) which is 


We thus see that these results, and hence the torques, are of the same order 
when all the material is captured as when no material is captured. It would then 
appear that the torque will be of the same order whatever fraction of the material 
is captured by the revolving bodies. 

7. Concluding remarks.—Although, in the cases we have considered, the inter- 
action of non-revolving material with a binary causes a reduction in the angular 
momentum of the latter, we do not claim to have shown that this is always the 
case. In fact, there appears to be no dynamical principle which would be 
violated if such were not the case. 

The work described in the present paper was undertaken with the object of 
checking the mathematical basis of the work in (1) but the problem of angular 
momentum exchange in the gravitational interaction of a particle with a revolving 
system is in itself a problem of purely dynamical interest and may be worthy of 
further investigation when high-speed computing facilities become more generally 
available. 


I would like to express my thanks to Professor W. H. McCrea of the Royal 
Holloway College who is supervising my work and to Mr R. Brooker of the 
Manchester University Computing Machine Laboratory for permission to use 
the machine. Also, I am grateful to the referee who suggested some improve- 
ments in the paper. 


Royal Holloway College, 
Englefield Green, Surrey : 
1953 May 6. 


References 


(1) Mainly work not yet published. 

_ (2) E. Strémgren, Kopenhagen Meddelauk, No. 39, 1922 (also later numbers). 

(3) S. Gill, Proc. Camb. Phil. Soc., 47, 96, 1951. 

(4) K. N. Dodd and W. H. McCrea, M.N., 112, 205, 1952. 

(5) E. T. Whittaker and G. Robinson, The Calculus of Observations, p. 264, Blackie, 1940. 


ON A PROBLEM IN THE THEORY OF FORMATION OF 
ABSORPTION LINES 
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Summary 


The paper gives an exact solution to the problem of determining 
the intensity of radiation emerging from a semi-infinite plane-parallel 
atmosphere when the Planck function B(v, T) varies linearly with the optical 
depth, the ratio of line to continuous absorption is constant with depth, and 
when the line scattering is partly coherent and partly non-coherent. The 
non-coherent part of the scattering is taken to be proportional to the coefficient 
of line absorption and to the mean intensity of radiation at the centre of the 
line, a process which is a simplified form of completely non-coherent 
scattering. The combination of both scattering processes in a suitable 
ratio describes the formation of subordinate lines and, to a certain degree 
of approximation, resonance lines broadened by collisions. The solution 
in the general case involves auxiliary functions which satisfy two 
simultaneous linear integral equations, but when the non-coherent part of 
the scattering alone is present, the solution only depends on the H-functions, 
their moments, and the quadrature of a simple integral form. 


1. Introduction.—\n the theory of formation of absorption lines exact solutions 
have been obtained for the intensity of radiation emerging from a semi-infinite 
plane-parallel atmosphere in which »,, the ratio of the coefficients of line absorption 
/, and continuous absorption k, is constant with depth, for two cases, one in which 
the line scattering is coherent, and the other in which it is completely non-coherent. 
Chandrasekhar’s (1) well-known exact solution for coherent scattering was origin- 
ally worked out for the case in which the Planck function B(v, T) could be written 
down as a linear function of the optical depth, and it has recently been extended by 
Su-Shu Huang (2) to include higher-order terms in the expansion for B(v, T). 
The problem becomes more complicated when the emission is completely non- 
coherent, since we have to find the solution of an integral equation over all 
frequencies. ‘The exact solution for this problem has been given by Sobolev (3) 
when B(v, 7) is a constant, and by Busbridge (4) when B(v, 7) is linear, but the 
complexity of the solution suggests that it might be worth while to look for a 
simplification of the non-coherent process which would preserve the main 
features of completely non-coherent scattering and yet make the solution more 
tractable. ‘This simplification becomes even more necessary if one wishes to 
consider the more general case in which the scattering is partly coherent and 
partly non-coherent. 

It is convenient to describe non-coherent emission, in the general case, by 
means of the function p(v,v’). This function gives the probability that the average 
atom in an assembly of line absorbing atoms will emit radiation in the frequency v 
after the absorption has taken place in the frequency v’. It follows from 
equilibrium considerations that p(v, v’) satisfies the relation 
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The emission j, in the frequency v per unit solid angle per unit mass of atmosphere 
is given by 


Jy=(1-€) | p(y dv’, (2) 
“0 
where ¢ is Eddington’s collision constant and J,, the mean intensity of radiation, 
I 
J,=5| (3) 


Let 7 be the optical depth in the continuous spectrum; then Eddington’s torm of 
the equation of transfer becomes 


is detined by 


“=z 


We take Biv, T) =a + br, (5) 
where a and b are slowly varying functions of v which are, in fact, quite insensitive 
to v over the line-width. 

When the scattering is the sum of coherent and non-coherent processes we have 


=2y(v, v’) + v’), (6) 
where « and f are constants and 6(v—v’) is Dirac’s delta function. ‘The first 
term of equation (6) gives the non-coherent component of the emission probability 
and the second term the coherent emission. 

If we let /,°(u) be the solution of equation (4) when «=o and f=1, and 
the solution when x =1 and =0, it follows that 
where «+ 8=1, is not a solution of equation (4). ‘Thus, in any combination of 
coherent and non-coherent scattering in the form (6), the solution of equation (4) 
is not simply a linear combination of the two solutions mentioned. We therefore 
have to solve the equations as they stand, but in so doing we may treat certain 
functions (for example, the H-functions) which occur in the coherent problem as 
known functions. 

2. The scope of the paper.—In this paper we give the exact solution of equation 
(4), with the condition (5), for the emergent intensity /,(0, ~) when the non- 
coherent part of the emission probability y(v, v’) is simply (7,/9)4(v% — v), that 


M being the value of 7, when v=vy. ‘The equation of transfer is then 

+ — (1 )(a + (5) 
where « + 8 = I according to equations (1) and (7). 

In Section 3 a new form for the law of diffuse reflection is given in equation (15). 
This is used to obtain the emergent intensity /,(0,) in terms of the scattering 
function for an atmosphere which scatters radiation according to equation (7). 
Section 4 deals with the scattering function. A general form for the scattering 
function for any p(v, v’) is found in Section 4.1 from the invariance of the law of 
diffuse reflection. The general form is given by equation (28), and the scattering 
function when f(v, v’) has the form (7) is deduced from it in Section 4.2. ‘The 
solution for /,(0, 1) is completed in Section 5, and is given by equation (50) for the 
general case x and f not zero. ‘The case 8 =0 is considered in Section (5.2) where 
the residual intensity 7(u) for radiation emerging in the y-direction and the 
residual intensity R in the emergent flux are found. 


: 
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The physical significance of the form of the equation of transfer discussed in the 
paper is as follows :— 
(i) Subordinate lines. When a line is formed by transitions between two 


broadened states, it may be shown (5) that the line emission in the frequency v by 
unit mass of atmosphere is 


3, By 
—e)l, (5 (9.1) 


where 6; and 6, are the half-widths corresponding to the lower and upper states 
respectively. ‘lhe emission process (9.1) is described by equation (7) when 
8; 

Oj; +O, 


p= 


(1i) Resonance lines. When the upper state of a resonance line is broadened by 
collisions, the line emission, according to Zanstra (6), is given by * 


5, 
where 6, and 6,, are the half-widths due to collisions and radiation damping, and 
J =§l,J,dv/fl,dv. Now 1, is extremely large near the centre of the line and falls 
off rapidly as vy — vp increases. Consequently, most of the weight in averaging over 
the line absorption profile comes from the region of the line centre. ‘This means 
that J =./, is likely to be a good approximation for the non-coherent part of the 


emission (9.2), and the total emission in the frequency v is described by equation (7) 
when 


6,40,” 6, 


(11) Lines formed by completely non-coherent scattering. ‘The line emission in 
this case is given by j, =/,J so that with J=-./,, the scattering is described by (7) 
when x=1 and 8=0. The equation of transfer with the line emission given by 
1,J,, has been discussed qualitatively by Spitzer (8). ‘The approximation only 
applies to strong principal lines with large collision damping. 

3. The method of solution.—Consider the equation of transfer in its general 
form 


dl, te 
(1 + ~ (1 €) | PY, v y dv + en, (a + br), (10) 
where /,(0,—’)=0 for o<p’ <1. ‘The angular distribution of the emergent 
intensity /,(0, «) can be found without actually solving equation (10) by appealing 
to “ an invariance arising from the asymptotic solution atinfinity”’. ‘This principle, 
which was enunciated by Chandrasekhar (9), is the idea behind the method used 


recently by Busbridge (4) to whom its introduction into the theory of line 
formation is due. 


Let =a+b (x1) 


represent the solution of equation (10). ‘This consists of the solution for an 
infinite unbounded atmosphere and a term /,*() which represents the departure 


* According to Miinch (7) it is not reasonable to include (1—e) in this case. 
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of the asymptotic solution from the value of /,(u) as we approach the boundary. 
Substituting (11) in equation (10) and using the relation (1), we find that /,*() 
satisfies the equation 
dl,* 


with the boundary condition 
T,*(0, —p’) (0<p' <1, O<v< 00). (13) 
From equation (II) we have 


1,(0, 4) =a + bu/(t +,) + 1,*(0, 1), (14) 
where the intensity /,* (0, «) arises from the diffuse reflection of the inward 


radiation given by equation (13). Let the intensity of diffuse reflection from the 
atmosphere be given by 


1,*(0, = =|. v’) dv 


where S(v,v’; u, p’) is the scattering function for the atmosphere. ‘The law of 
diffuse reflection as given in equation (15) is a useful general form which enables 
one to deal with any redistribution of frequencies by the diffusely reflecting 
atmosphere. This is made possible by the introduction of p(v, v’) under the sign 
of integration. If p(v,v’) is a single delta function we recover the usual 
formulation. 

Equations (14) and (15) complete the solution for the emergent intensity 
7,(0, w) in terms of the scattering function. When p(v, v’)is given by equation (7), 
we have 

Ny 


, 
T,(0, =a+ + | —a)S(v, v9 du 
No-o 


= f,—a)S(y, v; du’, 


where I I 


bo 


At the centre of the line where v = vy, equation (16) becomes 


I rl 
T,,(0, 4) =a + + | —@)S(v9, ¥93 dp’, (17) 


since x+$8=1. Equations (16) and (17) give the formal solution of equation (8) 
in terms of the scattering function for the atmosphere. 

4. The scattering function.—In this section we find the form of the scattering 
function S(v, v’; u,u') using the principle that the intensity of diffuse reflection 
from a semi-infinite plane-parallel atmosphere is invariant to the addition of a 
layer having the same optical properties as the atmosphere itself. ‘The method was 
clearly outlined by Ambarzumian (10) in his original paper, and has been discussed 
at length by Kourganoff (11) in the conservative case. An elegant method of 
deriving the equation satisfied by the scattering function has been given by 
Chandrasekhar (12), but in this paper we have used Ambarzumian’s method. 
We first establish the equation for the scattering function for any p(y, v’). 


& 
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4.1. A general form.—Consider the diffuse reflection of radiation by a semi- 
infinite plane-parallel atmosphere illuminated by a parallel beam of mono- 
chromatic radiation making an angle cos"! with the inward normal to the 
surtace, and let F,, be the flux in the beam. The intensity distribution of this 
field of illumination can be represented by 


where 6 denotes Dirac’s delta function. ‘The intensity of diffuse reflection is 
given by equation (15), namely 
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\/ 
Vv 


Fic. 1.—A_ layer of thickness dz overlying a_ semi-infinite plane-parallel atmosphere 
illuminated by a parallel beam of radiation of which one ray is drawn. The diagram shows the 
intensities emerging from the layer due to scattering in the layer and by diffuse reflection from 
the atmosphere. Intensities of order dt*® are not shown. The separation of the layer from the 
atmosphere is artificial. 


According to equation (12) the atmosphere whose diffusely reflecting properties 
we are considering has an extinction coefficient k in terms of which the optical 
depth 7 is measured, and a line-absorption coefficient /,(=ky,). The only 
emission which takes place in the atmosphere comes from line-forming atoms 
themselves. 

If we now add a layer of thickness dr containing absorbing and scattering 
material with the same “‘ optical properties”’ as the atmosphere itself *, the parallel 
beam of incident radiation described by equation (18), and of which one such ray 
is shown in the figure, will be diffusely reflected again according to equation (12). 
But the original ray traverses the path AB in the layer, is diffusely reflected at 
depth dr and then transmitted along CD finally to emerge from the layer with an 


* This restricts the discussion to the case in which 7,(=/,/k) is constant with depth. 
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intensity distribution .¥ ,“() given by 


The second term on the right-hand side gives the depletion, to the first order in dr, 
due to absorption in the paths AB and CD in the added layer. But the intensity 
J ,“(u) is restored to /,4() by the appearance of additional radiation in the 
p-direction due to scattering in the layer. ‘The components /,(j) and /,'?() 
shown in the figure arise from radiation scattered in the path AB, whilst /,(,) 
and /,() come from the path CD. ‘These components combine to restore 
JS ,“(u) to the intensity /,“(u). ‘Thus we have 


=F + + + + (22) 
From equation (2) the emission within the layer is 


dr(I+n, 1+, 
( 


j=) 


where J/, is the intensity of the absorbed radiation averaged over the directions of 
incidence. Now the mass of atmosphere within an elementary cylindrical volume 
with its axis in the x-direction is dr uk. Hence the intensity emerging from the 
layer, due to scattering within it, is 


dr 


Of the four intensities required in equation (22), 7,%) and J,®(yu) arise from 
upward scattering from the paths AB and CD respectively, and /,®(u) and 
T,®() from downward scattered radiation which is diffusely reflected at depth dr 
and then transmitted through the laver. ‘There are additional scattered intensities 
present, but they are of the second or higher orders in dr and are not needed. We 
now calculate the four components mentioned. 


dr 
(i) Tw) =(I | dv’, 
where Jy 1,(0, 


dr) (* 


(24) 
dr ;* 
I 


“1 
where | T,"(u') du’ + O(dr), 
0 


and /,“(y’) is given by equation (20). Hence 


=(1 —e)F,, =| Pt, v )p(v dv S(v » M4) > + O(d7*). (25) 


Ae 

‘ 
3 
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(iv) The calculation of /,() is similar to that for J,°(). We have 
Pl», v’) dv S(v, v's P(r’, v")p(v”, vy)n,- dv 
” ” du” 2 
x | a) + O(dr?). (26) 
0 
Now from equation (21) we have 
&fi+n, 1+”, 
+ + + LH), (27) 
by equation (22). Summing the four components given by (23), (24), (25) and 


(26), equation (27) gives, in the limit as dr—o, the following general equation 
satisfied by the scattering function : 


P(r, v;) (=* + + (v Vise Hy) 


by 


, , , , d 4 
0 “0 
+| pl», v’) dv | S(v, v | P(r’, pv", vy)n,- dv 
0 0 2h Jo 
1 ” 
x | S(v", v5", (28) 
0 


4.2. A special case.—In applying equation (28) to the special case in which 
P(v, v’) is given by equation (7), namely 


P(v, — + v’), 
we make use of the following properties of the delta function given by Dirac (13): 


— v')d(v, — v’) dv’ = — v’ )d(v, — v’) dv’ 
= — ¥1) = — 


My (1 —e) 


Thus we obtain 
/9)8(v9 — v4) + »4)] (= 
= v,) + — 


“1 du’ 
0 
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If we write 
(- =) S(v, Ho) 
+ 
Ho No(I —€) 


= A(v, v9) + {1 +B S(v, v5 Hy + 5 H's Ho) 
+ {1 + S(v, S(v, 
+ Bie S(v, v SE} sto) (30) 


and 


2) 


equation (29) for the scattering function may be written 


a A(v, v9)5(v% — vy) + BB(v, v)8(v — =0. (32) 
Take the following representation for the delta function : 
— v,) =0, vy >v+ 
(33) 


=O, vy <v— 
where Av is as small as we please, and let us first suppose that 
— 
We then have 8&(v—v,)--0 provided vy+Av<v<vy—Av. Under these 
conditions, equation (32) gives A(v, vy) =0. Now suppose that 
— vy> v9 + $Av. 
We then have 5(v) — v,) = 0, and from equation (32) B(v, v)=0. Since Av may be 
made as small as we please, equation (30) with A(v, v)) =o gives the scattering 
function S(v, vp; 4, 49) for all yAvy and equation (31) with B(v,v) =o gives the 
scattering function S(v,v; 4,49) again for all vAvy. To find the scattering 
function S(vp, v9; 4, 49) we only have to integrate equation (32) with respect to 
v,. This integration gives 
2A(v, v9) + =0, (34) 
and the required function follows from equations (30) and (31) with v = v9. 
4.2.1. The scattering function S(v,v;,') for all v.—When v+vy we have 
B(v, v) =o, and equation (31) gives 


=) S(v, v5 Mo) 
Ny Fo 

and when v =v, we obtain a equation (34) 


I+ (- + =) S(¥%, Yo; Mo) 
To Fo 


1 d 1 d 
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Equations (35) and (36) define coherent scattering functions. These functions 
are symmetrical in pu and py. 
The scattering function S(v,v ; 4, 49) for all vy may be written in the form 


S(v, v5 Mo) =(I—€)(I— 2» Ho) (37) 
where 


where 6, ,, is zero except at v = v) when it is equal to unity, and 
1 du’ 


4.2.2. The scattering function S(v, vo; 4,9) when v~vo.—If we impose the 
condition vv», the non-coherent scattering function S(v, v9; 4,49) follows at 
once from equation (30) with A(v,vy)=0. Now from equation (38), and the 
symmetrical property of the coherent scattering functions, we obtain 


(39) 


I+ ¥05 H's a) = Ho), S(v, v; 


where 


Hence 


= H( Ho) + He Ho) | 


(— 2) S(v, 3 Bs Ho) 


+ 0's 0) (40) 
No Jo 


Hu)= St, (41.1) 


S(v, Ho Ho) = 1) H (2, Ho) K Ho)s (42) 


where 


* 


Introducing the scattering function from (42) into equations (41.1) and (41.2), we 


obtain 
cl , 


ul) du! 
0 +p’l, 


(43) 


Wu) =H (ay, — 1] + 1) | (45) 


* When we have — =(1—e)(1 — 
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where 
(wo, du’ 


TH(w,, dp’ 
o 


I 


Writing 
fw) = — 1], (47) 
equations (44) and (45) become 


“1 
Hu) + | We’) dy’ (49-1) 


1 
= + | dy’ (49.2) 


Equations (42), (43), (49.1) and (49.2) complete the solution for S(v, v9; 1, Ho) 
when vr. 


5. The emergent intensity 


5.1. The solution in the general case.—The solution of equation (8) for the 
emergent intensity /,(0, 1) is given by equations (16) and (17) in terms of certain 
integrals of the scattering functions S(v, v; and S(v, v9; The formal 
solution may now be completed. Evaluating the integrals involving S(v, v; u, 4’) 
we find 


rl 


where £,! is the first moment of H(w,,). Thus 
| — a) S(v, v5 pe’) dp’ 


= bbw f,H(a,, (a+ bpl,)[ H(a,, — 1]. 
Furthermore 


— | — a) S(v, vy 5 de 


2H No-o 


Hence the result 


{ 


where w, = —€)(1 =(1 —€)(I — %), B, is the first moment of H(w,, 1), 
and K(, py’) is given by equation (43). The intensity at the centre of the line is 


given by equation (17) and is simply the emergent intensity for coherent scattering 
in the frequency vy, namely 


T,,,(0, = [2b a + (a + — (51) 
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where «,' is the first moment of H(w»,). This completes the solution in the 
general case. 

5.2. The solution in the special case 8 =0.—When f =0, y(u) =0 by equation 
(41.2), and K(y,») is a function of ~ only. From equation (44) we have 


#(u) =A(w)—1 since H(w,,4)=1. Hence K(u, =h(u) and equation (42) for 
the scattering function is then 
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5 Hs Ha) = H Ho): 
We now have 


i b 0 


where satisfies 
(1) =1 + | (53) 
We now evaluate 
Writing 
pl, + ply pl, + ply’ 
wt? 


7 7 C.+ "bes 
we obtain 


pl, + 


+ 
‘This gives 
1 j 1 — 
— + — — al 1 — wo)! 
Hence 


where ¢, =1/(1 +m), ¢,=1/(1+7,), and a! is the first moment of the H-function 
for an albedo =(1—e)(1—%). ‘The function is given by 
and satisfies the relation 

<h(u) <(1- am) 
the equality applying when ¢, =o, (I — aw)!” being the asymptotic form of h(u) 
for small (,/Z,. 


5.2.1. The residual intensity.—In the continuous spectrum outside the line we 
have ¢,=1, and the angular distribution of the emergent intensity is simply 
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a+ bp, as it should be for a linear form for the Planck function. The residual 
intensity in the y-direction is given by 


=(a+ bu) * E 
1—¢, 
The emergent flux in the frequency v is given by 
“1 
F(o) =27 | T,(0, dy’ 


From equation (54) we obtain 


0 0 


a 
+ (1 —wo)"(ax,! + %,2) + (- 
( o) ( 27 3 
where «,! and «,” are the first and second moments of A(j). Since the value of the 
integral in the continuous spectrum outside the line is $a + 44, the residual intensity 


R in the emergent flux is given by 


+ (1)! + =e (; = a 
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Summary 


During 1951 and 1952 radial velocities were observed of the components 
of « Centauri near the time of periastron nodal passage which occurred 
1952°15. Special attention was given to obtaining accurate results for 
the difference in velocity between the components. From this difference 
and the known visual orbit a parallax is derived of 0”-776+0”-008 (m.e.). 
Radial velocities with respect to the Sun were determined for each component. 
By combination with early results obtained at Lick and the Cape a rough 
value for the mass. ratio is derived of m:/m,=0°4. The velocity of the centre 
of gravity of the system is determined at —- 22-7 km/s. 


Introduction.—The first radial velocities of « Centauri were determined by 
W. H. Wright (1) at the Observatory of the D. O. Mills expedition to the Southern 
Hemisphere at Santiago, Chile, in 1904. From the difference in radial velocity 


between the components, in 1904 about 5 km/s, combined with the known visual 
orbit, Wright obtained an absolute parallax of 0”-76 + 0-04 (m.e.), thus confirming 
beautifully the then only existing trigonometric result by Gill and Elkin (2), 
0"-76 +0"-015 (m.e.), determined with the heliometer, which result is again in 
good agreement with the more modern photographic determinations (3), 
0"-751 +0"-016 (m.e.). 

The next study of the radial velocity of « Centauri was made by J. Lunt at 
the Royal Observatory, Capetown, in the years 1904-15 (4). ‘The radial velocity 
did not change appreciably during this interval and the difference between the 
components was of the order of 4km/s. Lunt’s prediction of a difference in 
radial velocity at the periastron node near 1953 of the order of 13km/s is of 
special interest for the present article. With the arrival, early in 1951, of the 
Cassegrain spectrograph at the Radclitfe Observatory it was decided to make new 
observations of radial velocity, the system being so near the periastron node. 
As the percentage accuracy of the parallax derived from the radial velocity 
difference is about the same as that of this difference, it is evident that the present 
phase is the most favourable for a parallax determination by means of Wright’s 
method. It was hoped to obtain a result at least comparable in accuracy with that 
of the modern trigonometric determinations. 

The observations; relative radial velocity.—The spectrograph camera giving 
the highest dispersion (22 A/mm at Hy) has been used throughout. Ilford 
Process plates were chosen because of their high contrast and convenient low 
speed. ‘The projected slit width was 14 microns. Exposure times were of the 
order of 4 and of 20 minutes respectively for the brighter and fainter companion. 
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These were varied proportionally to an extent depending on the transparency 
and the seeing. The exposures on the two components were taken in quick 
succession in order to ensure to the utmost identical conditions of the spectro- 
graph. The observations were always made near the meridian; the extreme 
hour-angle occurring in the material is 1"40™. Even so the minimum zenith 
distance of « Centauri at Pretoria is still 34°, and a red fringe at the south side of 
the image due to atmospheric dispersion could always be discerned and was 
found troublesome for accurate guiding during the 1951 observations. ‘The 1952 
observations were therefore all made with a dark blue filter across the guiding 
eyepiece, which device made the accurate centring of the image on the slit very 
much easier. ‘The iron arc comparison spectrum was as usual exposed before 
and after the exposure on the star. It was purposely exposed much denser than 
customary in order to provide a great number of reference lines. Because it 
was specially intended to obtain an accurate result for the velocity difference 
between the components, the observations, the measurement and the reduction 
were arranged to achieve this by making the determination entirely differential. 

Measurement and reduction.—Only those pairs of spectra were measured in 
which the spectra of the components of stars and arc were sensibly of the same 
density. Only the sharpest lines were selected for measurement and the same 
lines were measured in each spectrum of the same pair and this practice was 
followed for the star lines as well as for the arc lines. | Contrary to current practice 
the spectra were measured in only one direction, which was the direction of 
increasing wave-length. As for both spectra the same part of the measuring 
screw and the same part of the wire were used, any systematic effect existing 
between the stellar and the arc lines should cancel in the velocity difference. ‘The 
differences in reading between the measures of the same arc line as measured in 
the two spectra were plotted against the reading. ‘The resulting reduction curve 
was then used to reduce the differences found for the stellar lines to differences 
in microns due to velocity. ‘These differences in microns were then converted 
into km/s by multiplication with an appropriate factor that was calculated as a 
function of the reading. In the whole process no accurate knowledge of wave- 
lengths was required either of the stellar lines or of the arc lines. Knowledge of 
the wave-length was needed only approximately for the determination of the 
conversion factors and this could be done with sufficient accuracy by taking these 
factors from a table prepared once and for all as a function of the reading. ‘The 
results for the velocity difference from nine pairs of spectra taken in 1951 and of 
five pairs taken in 1952 are shown in Table I. 

‘lhe first column gives the epoch. ‘The second column shows the observed 
velocity difference in the sense: (fainter star) minus (brighter star), whereas the 
number of stellar lines measured is shown in the last column. ‘The small effects 
due to diurnal motion did in no instance amount to more than a few hundredths 
of a km/s and have been taken into account. ‘The effects due to orbital motion 
of the Earth around the Sun were negligible in all cases. 

Professor Hertzsprung kindly sent me in advance of publication his latest 
orbit of « Centauri, for the determination of which the accurate series of photo- 
graphic positions derived from plates taken at the Union Observatory in 
Johannesburg and the Bosscha Observatory in Lembang had been included. 
By means of this orbit the radial velocity ditlerence was calculated for a parallax 
of 1” of arc for the epochs given in Table I. The results are shown in the third 
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column. ‘The parallax was then found as the factor which changes the values of 
column 2 into those of column 3. The result is: 0’*776 + 0"-008 (m.e.). 
Column 4 of ‘Table I gives the velocity difference as derived from this parallax 
and Hertzsprung’s orbit and column 5 shows the residuals. It is somewhat 
surprising that the 1952 residuals are so much smaller than the 1951 values. 
Part of the improvement may be due toa readjustment of the spectrograph done 
early in 1952 or to the use in 1952 of the blue filter over the guiding eyepiece 
or to both. However, it seems more likely that the effect is due to chance, and 
the value for the parallax has been derived giving equal weight to all the obser- 
vations. It can be shown that the mean error of the parallax is almost entirely 


due to the uncertainty of the radial velocity determinations, the uncertainty due to 
the orbit being very small indeed. 


Radial velocity and parallax of « Centauri 507 


TABLE I 


The radial velocity difference (km/s) 


Obs. 

diff. diff. diff. 

f-—6b f—b 
Epoch 


1951-411 — 12°37 
“411 — 13°27 —12°64 —0'63 31 
“419 —12°22 —9g'81 —12°64 +0°42 68 
“419 —9g‘81 — 31264 39 
“419 —12°86 —9g'81 —12°64 —0'22 67 
“419 —12°12 —9g'81 —12°64 +0°52 47 
"452 —12°49 —g°82 —12°66 +0°17 58 
452 — 12°36 —9g°82 —12°66 + 0°30 62 
616 —13°51 —9g°86 —12°71 —o'80 35 


1952°192 —12°77 
433 —12:66 moe ce) 69 
“435 —12°74 —12°76 7O 
“485 —12°72 —12°76 +0°04 88 
“507 — 12:66 —12°76 +o-10 42 


Clearly any difference in wave-length of the same lines in the spectra of the 
two stars due to causes other than velocity would affect the present determination 
of the parallax; for instance, a differential relativity redshift may easily attain 
o-1km/s. Unfortunately the uncertainty in the masses and radii of the stars 
would make a prediction of this effect illusory. In any case the possible 
differential effects due to causes other than velocity could be eliminated by means 
of observations extending over a complete period. In this connection obser- 
vations of the radial velocity difference at those phases where the relative velocity 
is known to be zero are of interest. 

The radial velocity with respect to the Sun.—In order to determine the radial 
velocity of « Centauri with respect to the Sun, the 1952 spectra were measured 
in the conventional way, i.e. by measuring twice, once in the direction increasing 
and once in the direction decreasing with wave-length. It is known that by this 
practice some errors are eliminated. ‘The lines and wave-lengths as recommended 
by R. M. Petrie (5) for late-type stars and appropriate for the dispersion used 
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have been employed. As mentioned before, the iron arc comparison spectrum 
had been purposely exposed rather more heavily than usual. Accordingly only 
the fainter I.A.U. standards have been measured for the derivation of the 
reduction curve. 

In Table II the results of these measures are given. ‘The Cape series of 
observations was reduced on the Lick system by adding + 0-9 km/s to all velocities 
of the brighter component, whereas — 0-2 km/s was applied to the Cape velocities 
of the fainter star in accordance with the differences found between Lick and 
the Cape as given in Lick Observatory Publications, Vol. 18, XII. 


TABLE II 


Radial velocities of « Centaurt 


Brighter Fainter 
Epoch component component 


1952°192 18-6 33°1 
433 18-6 32°0 
435 
“485 20°6 32°7 
"507 18-4 32°0 


The relation between the radial velocity system of the Radcliffe Observatory 
and that of Lick cannot be determined with great accuracy yet. From a dozen 
radial velocities determined at Pretoria of stars for which results have been 
published by Lick, a mean difference of 0-2 km/s has been found. Consequently 
to the values of Table II +0-2km/s was added to reduce them provisionally 
on the Lick system. 


—I5 


4 
4 
—30 —20 km/s 
Fic. 1.—Radial velocity of brighter component (ordinates) against that of fainter (abscissae). 
O Lick (1904). 
@ Upper left: Radcliffe (1952). 
@ Lower right: Cape (1909-1915). 
Full-drawn line: relation between radial velocities of the components. Its slope gives m,/m,=0-4. 


Dotted line is the locus of equal velocity. It intersects the full-drawn line at the point representing 
the velocity of the centre of gravity. 


In the accompanying figure the radial velocities of the brighter component 
are plotted against those of the fainter star for all the results available. An 
attempt to derive the mass ratio from the slope of the full-drawn line in this 
diagram gave m,/m,=0-4, the fainter star being the lighter. For the radial 


—20 
km/s e 
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velocity of the centre of gravity the value —22-7km/s is obtained. The present 
result for the mass ratio is rather smaller than that found from meridian obser- 
vations (6) which is m,/m, =0-82. 

Clearly the mass ratio as determined above from radial velocities cannot be 
very accurate as a consequence of the uncertainty in the systematic corrections 
to be applied to the results obtained at different observatories in conjunction 
with the relatively small range of radial velocity due to orbital motion. Therefore 
there is no doubt that a more accurate value for the mass ratio will eventually be 
determined from parallax plates, though at present these results do not seem to 
be very consistent (7). 


I am indebted to Dr A. D. Thackeray for his interest and for discussion. 
My thanks are due to Dr Hertzsprung for allowing me to use his accurate orbit 
of « Centauri. I am grateful to Dr W. H. van den Bos and Dr W. S. Finsen for 
valuable information on the orbit and to Professor Stratton for discussion on 
spectroscopic questions. 


Radcliffe Observatory, 
Pretoria : 
1953 June. 
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THE ABSOLUTE MAGNITUDE AND SPECTRUM OF THE 
CLASS S STAR 7! GRUIS 
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(Communicated by the Radcliffe Observer) 


(Received 1953 August 31) 


Summary 


Spectrograms have been obtained of the typical class S star 7‘ Gruis, 
and its close, faint companion. The two stars have a common radial 
velocity. This and other evidence strongly suggests that the stars form a 
physical system. The companion is a Go V star and from this fact it 1s 
deduced that 7! Gruis has an absolute magnitude of between —1 and o 
(a normal giant luminosity). An investigation of the spectrum of the 
S-type star shows that the supergiant characteristics indicated by lines of 
SrII, Cal and Ball are probably illusory. The star is classified as of 
equivalent spectral type M7 III. This binary system apparently consists 
of a late-type giant of abnormal chemical composition and a solar-type 
dwarf of normal composition. 


Introduction.—The stars of spectral class S are characterized by the 


appearance of bands due to the zirconium oxide molecule in their spectra. 
As regards temperatures the S, N and M stars appear to form parallel 
sequences. Since the S stars are very few in number it has not been possible 
to use statistical methods (on proper motions, radial velocities, etc.) to obtain 
an estimate of their absolute magnitude. However, it is generally believed 
that these stars are supergiants. ‘The main reason for this belief appears to 
be the great strength of the Sr II 4077, 4216 and Ba II 4554 lines in their 
spectra compared with the strength of these lines in other late-type stars. 
The strength of these lines is generally a good indication of luminosity. In a 
theoretical study Wurm (1) deduced that the appearance of ZrO and the 
suppression of TiO in S stars as compared with M stars could be explained 
in terms of the supergiant characteristics of the S stars (i.e. low atmospheric 
pressure) and that it was not necessary therefore to assume an abnormal 
abundance of Zr in the atmospheres of these stars. ‘The work of Keenan (2) 
which showed an absence of supergiant irregular variables in the later 
sub-divisions of class M apparently supported this argument. More recently 
however Merrill (3) and Merrill and Buscombe (4) have shown that not only 
the ZrO bands but also the ZrI and Zr II lines are enhanced in S-type stars 
so that a high abundance of Zr seems to be definitely indicated. Nevertheless 
the assumption that class S stars are supergiants has still been maintained in 
order to explain the Sr II and BaIlI line intensities. ‘That the class S stars 
are very peculiar, and possibly relatively transient, objects is shown by 
Merrill’s discovery (§) of lines due to the element technetium in their spectra. 
No stable isotopes of technetium are known. Furthermore Bidelman (6) has 
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deduced that the ionized lines of the rare-earth elements (except europium) 
are much enhanced in the S star R Andromedae compared with the M star 
R Leonis. 

An opportunity for determining the absolute magnitude of the S-type star 
m Gruis occurs, since this star has a faint companion whose spectroscopic 
parallax (if the star forms a binary with 7! Gruis) may be used to determine the 
luminosity of the S star. This paper reports the results of such an investigation. 

a! Gruis and its companion.—The star 7! Gruis, « 22" 16™ 378, 8 —46°27'+1 
(1900), is an irregular variable of from 5-8 to 6-4 visual magnitude (7). It is the 
type star of class S as defined by the I.A.U. (8) and was used as such in the 
Henry Draper Catalogue (9). The star has a 10™-g companion 2”-71 distant, 
position angle 202°-1 (1943°49) (10). Observations since 1900 show no 
indication of any relative motion between the two stars. It was found possible 
in conditions of good seeing to obtain spectra of the companion free from that 
of the main star. Spectra of the S star have been obtained with dispersions 
of 30 A/mm, 50 A/mm and 96 A/mm at Hy and of the companion with 
dispersions of either 50 A/mm or 96 A/mm at Hy, the former for radial velocity 
determination, the latter (with the spectrum considerably widened) for spectral 
classification. ‘These observations were all carried out at the Cassegrain focus 
of the 74-inch Radcliffe reflector. 

‘The companion was measured for velocity using the solar-type wave-lengths 
recommended by Miss McDonald (11). The velocity of 7!Gruis was 
determined for both emission and absorption lines. In the latter case the 
wave-lengths used were selected from the list given by Merrill (12). The 
velocities so determined are entered in Table I and it will be seen that they 
all agree well with each other. 


TABLE I 
Velocity 
Emission Lines —18-8 km/s 
Absorption Lines —19°6 km/s 
Companion Absorption Lines —17:1 km/s 


7 Gruis 


In the case of two stars so close together and of great difference in magnitude 
it is possible for the scattered light of the bright star to cause trouble in 
obtaining a pure spectrum of the fainter one. It may therefore be as well to 
point out firstly that the observations were only made in the best seeing 
conditions and great care was taken both with the focusing of the stars on the 
slit and with the guiding, and secondly that tests showed the ZrO bands of the 
S star could be recorded with very short exposures indeed. The absence of 
these bands therefore forms a good check on the purity of the spectrum of the 
companion. 

The agreement of the radial velocities, the small separation of the two 
stars and their lack of relative motion all strongly suggest that they are 
physically connected. That the two stars form a binary system is therefore 
assumed and the results deduced in this paper are based upon this assumption. 
The probability that it is incorrect must be very small indeed. 

The spectrum of the companion having been provisionally classified as 
dwarf F8—G2, spectra of a number of standard stars in the system of the 
Yerkes Atlas (13) near these types and of various luminosity classes were 
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obtained and compared with the spectrum of the companion using a Hartmann 
spectrocomparator. ‘The companion was thus classified as Go V. Plate 4 
shows a spectrogram of the companion together with some stars of standard 
spectral type. According to the Yerkes Atlas the absolute magnitude of sucha 
star is between +4 and +5 (visual). Since the apparent magnitude of the 
companion is 10-9 and the apparent magnitude of the main star at maximum is 
5:8, the absolute magnitude of 7! Gruis at maximum must lie between — 1 and o. 
These figures may be slightly in error due mainly to the difficulties in estimating 
the apparent magnitude of a companion lying so close to a bright star. However, 
such an error will usually result in the brightness of the companion being 
underestimated, which would result in the luminosity of z' Gruis being 
overestimated. Furthermore 5-8 is probably an upper limit for the apparent 
magnitude of 7! Gruis. Mr R. P. de Kock considered the star much brighter 
than usual when he estimated it as near 6™-o in 1952 May. Thus it is unlikely 
that any revision of apparent magnitudes will make 7! Gruis more luminous 
than absolute magnitude —1. It is clear therefore that so far as luminosity is 
concerned z! Gruis must be regarded as a normal giant. It is certainly not 
a supergiant. 

It may easily be shown that if the binary orbit lies in the plane of the sky 
and its total mass is about 10 solar masses, orbital motion with a period of about 
5 600 years should be observed. No significant change in the position angle 
or distance between the two stars has been observed since 1900. If future 
observations did show orbital motion they would be of great interest, since a 
direct determination of the mass of an S-type star would be of considerable 
importance. 

The absorption spectrum of ~ Gruis.—Whilst a detailed investigation of the 
spectrum of 7! Gruis must be carried out with a high-dispersion Coudé 
spectrograph, an examination of spectrograms at 30 and 50 A,mm at Hy has 
been undertaken with the prime purpose of investigating those features which 
have previously led to the conviction that the S stars were supergiants. 

The M-type stars are generally classified by using the strength of the 
TiO bands in their spectra. However, Keenan (2) has given criteria whereby 
irregular and semi-regular variables of small range and type M may be classified 
both for type and luminosity using atomic lines in the ultra-violet region of 
their spectra. An attempt has been made to classify 7! Gruis with these 
criteria using the slight revision of Keenan’s types as recommended in the 
Yerkes Atlas (13). In general the ultra-violet regions of M and § stars are 
quite similar so that Keenan’s criteria should be applicable at least as a first 
approximation. The classification was made by preparing prints of 7! Gruis 
and comparing them with the spectrograms reproduced in Keenan’s paper. 
Whilst such a procedure is in general less accurate than the direct comparison 
of plates with a spectrocomparator it has the advantage in this case that since 
Keenan’s stars are all variable in light they may also vary in spectral type, and 
it is thus preferable to rely on the spectra actually used in the original 
classification of the stars. Furthermore Keenan’s prints show that some of 
the lines exhibit very marked changes with type and luminosity, so that no 
difficulty was found in using this method of classifying 7! Gruis. Spectra of 
several M-type stars classified by Keenan were also obtained and some of 
them are reproduced together with spectra of 7! Gruis in the plate. 
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Keenan gives the following ratios as indicators of type : 

(1) Ti 3982, 3989, 3999, 4009/Fe 4005. 

(2) Ca 4226/Fe 4143. 

(3) Cr 4254/Fe 4250, 4260. 

Criteria (1) and (3) were found very useful and indicate an equivalent 
spectral type of M7. ‘The line Ca 4226 is mentioned later; it is not desirable 
to use it to assign a spectral type to the star. Nevertheless it is apparent from 
the spectra which Keenan reproduces that Fe 4143 may be used to assign a 
spectral type simply by comparing it with the large number of other lines in 
its vicinity. In this way an equivalent spectral type of M7 was also obtained. 

Keenan lists the following line ratios as being sensitive to luminosity : 

(1) Sr II 4077/Fe 4045 or 4063. 

(2) Sr II 4216/Ca I 42206. 

(3) The blend 4172-5/Fe 4143. 

(4) Fe 4376, 4389/Fe 4383. 

Criteria (1) and (2) involve the Sr II lines and as noted previously it is the 
great strength of these lines and of Ba II 4554 which has led to the placing 
of S stars in the supergiant class. An examination of the spectrum of 7! Gruis 
shows clearly that if the ratios (1) and (2) are taken as true luminosity indicators 
the star would be of very high luminosity, certainly brighter than » Ceph 
(luminosity class Ia, absolute magnitude -—6-0). That such a _ high 
luminosity is probably illusory is seen when the ratios (3) and (4) are considered. 
They indicate that the star should be placed in luminosity class III. This 
would make its absolute magnitude about —0-2, which agrees well with the 
absolute magnitude for the star derived from the companion. 

Whilst it would be of interest also to apply these criteria to other class 
S stars it appears safe to assume that 7! Gruis is typical of the class and that the 
high luminosity indicated by the ratios (1) and (2) above are spurious for the 
class in general. It is of interest to note that considering 7! Gruis as a star of 
equivalent spectral type M7 III, not only are the SrIl and Ball lines 
abnormally strong but the Ca I 4226 line is abnormally weak, which of course 
enhances the spurious luminosity effect. One of the few possible causes for 
these abnormal intensities is an abnormal relative abundance of the elements 
though not necessarily of Ca, Ba and Sr. An abnormal abundance of some 
other element may radically affect molecular dissociation in a cool star and hence 
affect the intensities of the atomic lines of other elements. Bidelman (14) 
for instance tentatively suggests that an abnormal carbon abundance may 
explain the presence of ZrO and the absence of TiO in the spectrum of the 
peculiar star GP Orionis. Whatever the solution to the problem may be the 
present work appears to lend support to the suggestion of Bidelman and 
Keenan (15) that the class S stars are related to the Ba II stars (peculiar stars 
of class G and K showing abnormally strong Ba II, Sr II lines and carbon 
bands). 

Emission lines in the spectrum of 7 Gruis.—Neither the Henry Draper 
Catalogue (16) nor Merrill (17) record any emission lines in the spectrum of 
7! Gruis. | However, several emission lines have been observed during the 
course of the present investigation. The results are recorded in Table II. 

Mr R. P. de Kock reports that the star was 0™-4-0™-5 brighter than normal 
(that is to say about 6™-0) about the middle of 1952 May (about JD 2434150). 
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He has observed no brightening since then. From the magnitudes which he 
has kindly supplied it appears that the emission lines were observed during a 
period of decreasing brightness. The behaviour of the lines as judged from the 
available data is the same in 7! Gruis as in other Se variables (17). ‘The radial 
velocity derived from the bright lines is recorded in Table I. It does not differ 
significantly from that given by the absorption lines. 


TaBLe II 
Plate No. JD Emulsion Exp oe Region Remarks 
2434000°0 + (minutes) 
Cb 799 154°6389 K1o3a F 5 6750 A- Noemission lines seen. 
4500 A 
Cc 823 160-6083 K1o3a F 30 6750 A- HB emission possibly 
4oco A present. 
Ce 905 1936146 O 120 3000 HB, y, 8, €; Si3g05, Fe 3852, 
3850 A 4202, 4308 all in emission. 
Cc 998 228-5278 Kro3za O 10 5000 A~_——- Fe 4202, 4308 in emission; 
4000 A HB possibly weakly in 
emission. 
Cd 1612 581-4778 K1o3a O 30 5000 A~ 4202, 4308 weakly in 
3900 A emission. 
Cb 1619 -§85°6132 Kro3za O 105 5000 (Poor observing conditions.) 
4400 A HB possibly very weakly 


in emission. 


Conclusion.—F rom the results presented above it appears that class 5 stars 
must be considered as giants with visual absolute magnitude about 0 or —I 
and not as supergiants. ‘This reduces the possibility that the distinctive features 
of these stars can be explained as due to abnormal physical conditions in their 
atmospheres, and an assumption of abnormal relative abundances of the 
elements seems the only possible hypothesis. In view of this it is perhaps of 
interest, from the point of view of stellar evolution in general and the evolution 
of binary stars in particular, to note that there exists a binary consisting of an 
S star showing evidence of abnormal relative abundances of the elements and 
a solar-type dwarf whose spectrum appears quite normal and which thus gives 
no indication of any deviation from the normal solar-type relative abundances. 

Acknowledgments.—\ should like to thank Mr R. P. de Kock who kindly 
supplied visual estimates of the magnitude of 7! Gruis and Dr W. H. van den 
Bos who kindly supplied unpublished data on 7! Gruis as a double star. I am 
deeply indebted to Dr A. D. Thackeray for some very helpful discussions. 
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THE PHOTOELECTRIC LIGHT CURVE OF V525 SAGITTARII 


G. G. Cillié, Jnr. and E. M. Lindsay 
(Received 1953 April 1) 


Summary 


Photoelectric light curves have been obtained in blue and yellow light 
for the star V525 Sgr of B-Lyrae type. From observations of primary 
minima it is found that Min=JD 2 429 6624-459 3+-0%-705 12200 E. 
The difference between the maxima is larger in blue than in yellow light 
and the star is very slightly redder at secondary maximum than at primary 
minimum. 


The f-Lyrae star V525 Sgr (CPD —30°5841 =HD 177768) has been the 
subject of two papers by D. O’Connell. In the first (1) O’Connell gives a photo- 
graphic light curve based on 273 estimates, from which he derives the following 
light elements : 

Min=JD 2 427 1049-987 + 09-705 125 E. 


In his second paper (2), in which he uses 875 photographic estimates, he gives the 
following light elements : 


Min =JD 2 429 6621-459 3 + 09-705 121 OE, 


and confirms the essential features of his previously published light curve. In the 
second paper he uses his mean light curve to derive the orbital elements of the 
system. 

C. Payne-Gaposchkin (3) gives A2 as the spectral type of this star and from I 293 
estimates derives a period of 04-705 121 40 for it and gives the following extreme 
magnitudes: M,8-13, M,8-20, m,8-96 and m,8-37. 

V 525 Sgr was observed on 27 nights during 1951 and 1952 with the Linnell-— 
King photoelectric photometer with 1 P21 tube attached to the Rockefeller reflector 
of the Boyden Station of the Harvard Observatory. Altogether 1 445 comparisons 
between the variable and a nearby standard star were made, of which 1 008 were in 
blue and 437 in yellow light. The blue filter was a Corning No. 5551; the yellow 
a Corning No. 3385. ‘The comparisons are expressed as magnitude differences — 
Am, and Am, for blue and for yellow light respectively. 

Primary minima of the variable were observed on 7 nights, and these runs were 
used to compute the following light elements : 


Min = JD 2 429 6624-459 3 + 09-705 122 00 E. 


Table I gives the differences between the observed and the predicted minima in 
heliocentrically corrected time. The weights in the third column, assigned to 
each determination of the primary minimum according to the number and the 
distribution of the observations of the variable in the vicinity of the minimum, run 
from I (very poor) to 5 (excellent). We regard the residuals to be within the errors 
of measurement. 


a 
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The above light elements were used to compute the phase of every observation 
of the variable using primary minimum as zero phase. ‘The observations, collected 
into normal points in the usual way, are given in Tables II and III. In both 
tables the mean errors of each normal point are given in the third columns and the 
number 2 of observations used to form each normal point in the fourth columns. 


TABLE I 
Observed and predicted minima of V525 Sgr 
Observed Minima E Wt. Calculated Minima 


d d 
2 433 854°408 2 2 433 854°409 6 
3 864°279 6 3 864-281 3 
3 866-396 4 3 866-396 7 
3 885-4379 3 885-4350 
3 890°371 1 3 890°3708 
4157°613 6 4157°6121 
4 203°4449 4 203°445 0 


TABLE II 
Normal points of the blue light curve of V525 Sgr 


Phase .e. Phase A 


©°000 4 0°559 2 0°677 

‘O10 3 593 4 643 

020 8 6317 “586 

4 ‘ ‘680 0 539 

"040 0 1 "505 

"050 I "749 5 493 

"069 3 ‘777 © 486 

9 ‘ ‘801 1 513 

"1259 ‘ ‘824 5 ‘531 

“154 6 ‘ 844 2 562 

“188 8 ‘ ‘868 9 ‘612 

"231 3 ‘ 5 681 

276 8 ‘922 1 “787 

"336 6 "940 8 ‘857 

"451 3 ‘ 965 1 

"486 3 ‘974 5 1°073 

‘ 0-986 7 +0°004 
0°539 0 + 0-003 


Our mean light curves are fairly smooth throughout, small irregularities being 
probably due to accidental errors. Our curves confirm at least qualitatively the 
essential features of those given by O’Connell, the most important difference 
being that we find primary minimum deeper in comparison with primary maximum, 
but then not quite so deep as was found by Mrs Payne-Gaposchkin. We find 
the difference between primary maximum—which occurs after primary minimum 
—and secondary maximum somewhat larger than O’Connell, but not quite so 
large as was found by Mrs Payne-Gaposchkin. 

On account of the flatness of both maxima and of secondary minimum and 
because of small irregularities in those neighbourhoods, we used a graphical 


| 

O-C 
d 
4 
7 
—0-000 3 
+0002 9 
+0000 3 
-O°OO! § 
—0'000 I 

M.e. of 

16 
35 
29 
38 
24 
33 
35 
35 
28 
36 
39 
37 
38 
16 
18 
17 
19 
23 
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method (4) with our mean light curves to find the exact phases of the maxima and 
the secondary mimimum. This was done for both colours separately, the results 
being given in Table IV. In the bottom row of the table the means of the phases 
are given, determined by giving double weight to the determinations from the blue 
curve. If, following O’Connell, we take zero phase at secondary maximum, the 


III 
Normal points of the yellow light curve of V525 Sgr 

Phase M.e. of Phase a, M.e. of 

(Period) Norm. pt. (Period) ¥y Norm. pt. 
m m m 

4 +0004 0°564 5 0°530 +0-002 

026 9 005 "592 9 -496 "002 

"052 4 005 626 3 442 003 

092 I -006 680 5 -389 “002 

"133 5 ‘005 ‘722 2 "364 

"171 6 "004 ‘7597 "353 

6 003 ‘792 7 +362 

‘2522 823 5 “402 

‘3147 ‘003 ‘851 8 ‘439 

*396 0 "002 881 4 

4437 ‘ "003 ‘9169 622 

‘475 8 "003 946 8 ‘753 

"504 6 "003 °°977 3 
0°535 5 + 0-003 


TABLE IV 
Phases and magnitudes of maxima and minima of V525 Sgr 


Wt. Primary Primary Secondary Secondary 
for Minimum Maximum Minimum Maximum 
Phase 


Phase Am Phase Am Phase Phase 


Blue 0'760 


Yellow 0-261 0-764 
Mean 
Phase 0°257 0-761 


‘TABLE V 
Magnitude ranges of V525 Sgr 


Prim. Max.| Sec. Min. | Sec. Max. |Prim. Min. 
—Sec. Min.|—Sec. Max.|-Prim. Min.| —Sec. Min. 


m m m m 


Blue 0'259 0200 ©°459 


Yellow 0°243 0640 0-448 
Payne- 
Gaposchkin 0°24 0:76 0°59 


n 
19 
19 
17 
17 
16 
22 
22 
19 
21 
15 
16 
12 
15 

Am 

0°353 

Z Prim. Min. Prim. Max. 
—Prim. Max. —Sec. Max. 

0:07 
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following phases are obtained for primary minimum, primary maximum and 
secondary minimum, respectively: 0-239 (0-244), 0-496 and 0°753(0-760). The 
values in brackets are those given by O’Connell. According to our results, 
maxima and minima follow each other very nearly at quarter intervals, the only 
significant exception being the primary minimum. 


AN. 


0-4 


\ 
\ 


1-0 
0-0 0-2 0-4 0-6 0-8 0-0 
Fic. 1.—Upper part: Blue light curve of V 525 Sagittarii, Lower part: Yellow light curve of 
V 525 Sagittarit. 
Ordinates: Magnitude differences. 
Upper part: Amp. 
Lower part: Amy. 
Abscissae: Phase. 


With two different maxima and two different minima, it is possible to compute 
six different amplitudes for V525 Sgr. ‘These are given for blue and for yellow 
light in Table V and compared with the corresponding amplitudes given by 
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Mrs Payne-Gaposchkin. From Table V it appears that primary minimum is 
slightly deeper in blue light thanin yellow light. Our results in the last column of the 
table, giving the differences between the two maxima, in the two colours, bear out 
a conclusion reached by O’Connell in a study of a number of other B-Lyrae stars 
(5), namely that the difference between the maxima is bigger in blue than in yellow 
light. According to our results the star V 525 Sgr is slightly redder at secondary 
maximum than at primary maximum, but only very slightly. Our results for this 
star do not appear to fit O’Connell’s empirically established relation for this effect 
very well. 


In conclusion, it is a pleasure to record our thanks to the staff of the Boyden 
Station and especially to the Chief Assistant, Mr M. J. Bester, who, amongst other 
things, observed the variable for us on a number of nights. 
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MONOCHROMATIC MAGNITUDES OF MARS IN 1952 
R. v. d. R. Woolley 


(Received 1953 September 16) 


Summary 


Monochromatic magnitudes of Mars were determined from spectra taken 
on eleven nights near the opposition of 1952. ‘lhe plates were measured at 
seven wave-lengths from A 4050 A to A 6360 A. Monochromatic coefficients 
of phase and values of the Russell-Bond albedo are calculated. 


At the suggestion of Dr G. de Vaucouleurs, the planet Mars was put on the 
programme of observations of monochromatic magnitudes of bright stars carried 
out on Mt. Stromlo in 1952 and to be described by Woolley, Gascoigne and A. de 
Vaucouleurs. Spectra of the planet were photographed on 15 plates taken between 
1952 March 8 and 1952 July 19. The spectra were measured at the six wave- 
lengths used in the stellar programme (AA 4050 A, 4250 A, 4550 A, 4945 A, 5980 A, 
6360 A) and also at A5430A. ‘This last wave-length was discarded in the stellar 
programme as most of the stars were too faint to record a suitable photographic 
density at this wave-length, at which the plates used were slow: but for Mars 
and for the bright stars used as standards for Mars the densities were quite 
sufficient. 

The stars with which Mars was compared are listed in ‘Table I, which shows 
the magnitudes adopted at the seven wave-lengths chosen. 


"TABLE I 


Magnitudes of comparison stars 


Star 4050 AN 4250A 4550A 4945A 5430A 5980A 6360A 


B Car 1°59 1°58 1°60 1°58 1°58 1°58 


1°56 
A Cen 1°97 201 2°02 2°05 2°02 2°02 2:01 
« Vir 0°68 0'73 0°76 0°95 
B Cen 0°32 0°37 O'55 0°63 
a Oph 2°16 2°53 2°09 1°97 1°93 1°87 
e Sgr 1°70 1-72 1°72 1°70 
Aql 0°94 0-90 082 0°77 067 0°59 0°55 
a Pav 1°65 1°69 1°82 1°88 1°88 
a Gru 1°54 1°59 1°63 1:66 1°73 
a PsA 1°16 1°09 1-06 1‘07 


‘The magnitudes of the comparison stars are all taken from the main stellar 
programme except those for 45430 A which were found by interpolation among 
the magnitudes for the other wave-lengths. ‘The observed magnitudes of 
Mars, in the seven wave-lengths measured, are shown in ‘lable II. ‘The weights 
were computed as follows. In any particular wave-length, the weight of a plate 


36 
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is (1 m+1/n) 1, where m spectra of Mars and n spectra of comparison stars are 
measurable at this wave-length on this plate. If there are two or more plates 
on a night the weight for the wave-length is the sum of the weights of the plates: 
and the weight shown in the table is the average weight of the seven wave-lengths, 
rounded off to the nearest integer, except on June 3 where it is the average of 
tive wave-lengths. 


Taste II 
Observed magnitudes of Mars 
Date 
1952 4050 A 4250 A 4550A 4945 A 5430A 5980A 6360A Wt. Plates 


March 8-74 1°55 1:36 —O'17 —0'°93 —1'16 2 56 
12°71 1°65 1:28 0°73 —0°30 —O'99 —1°32 
25°74 0:98 085 0:28 —0'23 —o-69 —1°28 —1°58 2 76 


N 
wm 


May 3:64 —o-'12 —0:30 —o-91 —1°37 —1°84 —2:24 —2°66 5 89,90 


5°58 0°05 —o-62 —1°:28 —1:66 —2:07 —2°56 1 95 
21°66 0°32 —0O°52 —1'05 —2°08 —2°30 5§ 100,101 
June 3°47 0°65 0°40 —0'25 —0°67 1°29 on . 4 103, 104, 105 
8-46 o'88 0°63 —0O°39 1°58 —1°87 2 106 
20°46 1°09 o'84 0°16 —0°65 1°40 
July 19°46 1°66 5°34 +0°72 0°26 —I°I4 —1°59 2 128 


Plates 103, 104, 105 overexposed in the red. 


‘The magnitudes were reduced to unit distance of Mars from the Sun and 
also to unit distance of Mars trom the Earth by adding —5log RA, where R is 
the heliocentric and A the geocentric distance of Mars, both in astronomical 
units. ‘lable 11] shows the reduced magnitudes. 


III 


Magnitudes of Mars reduced to unit distance from both Earth and Sun 


Longitude 
Date of central 
1952 slog RA 4050A 4250A 4550A 4945A 5430A 5980A 6360A Phase meridian 
w 
March 8°74 +0°759 + 0°03 —0-93 1°69 1°92 31-9 204 
+0670 0.98 0°06 0°33 —0-97 —1°66 1°99 30°5 157 
25°74 +0:362 o62 40°29 —0°08 —1°36 —104 25°71 47 
May 364 —0'280 016 —o'02 —0°63 1°09 —1°56 —1°96 2°38 2°5 27 
5°58 0°34 —0°33 —0'99 —1°37 —1°78 —2°27 346 
21°66 —0-262 0°58 0°40 —0°26 —I1'19 —1°82 —2°04 17'0 237 
June 347 0°78 —1°16 25°60 51 
20°46 O99 O74 O06 —O29 —0'74 —1°50 1°57 25 
30°55 0'224 0°79 0°05 —0'29 —0'77 —1°45 1°72 364 19 
July 19°46 1:16 084 0°22 —0'24 —0°63 —1°64 —2:'09 39°6 341 
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‘The magnitudes at unit distance, shown in Table III, can be analysed accord- 
ing to the phase of Mars, ~, equal to | 4g — Aol. 

If in each wave-length m=g +x, then yp is the coefficient of phase, and g 
the magnitude reduced to unit distance and zero phase. Least squares solutions 
of the data given in Table III give the coefficients of phase and reduced magnitudes 
shown in Table IV. ‘The residuals from m=g+j% were plotted against the 
longitude of the central meridian, which is listed in Table III, but the plots 
showed no well-defined relation, even in the two red wave-lengths, where a 
variation of the magnitude of the planet depending on the longitude of the central 
meridian with an amplitude of 0™-2 has been reported (cf. Guthnick and Prager, 
Veréff. Berlin-Babelsberg, 2, 1918). 


TABLE IV 
Reduced magnitudes and coefficients of phase 


Wave-length g pe 


m 


4050 0°03 -0'0246 + 0013 
4250 —0'02+ 0°02 0°0215+ 0011 
4550 —0'63+0'02 00212 + 0007 
4945 —1'14+0°02 0:0237 + 000g 
5430 —1:60+ 0°03 0°0213 + 0013 
5980 —1'98+ 0:03 0-0115 +0012 


6360 ~2°31+0°05 0°0134 + 0020 


The albedo of Mars can now be calculated. ‘The work given here follows 
Russell’s classical work on the subject (Astrophysical Journal, 43, 177, 1910). 
If A denotes the albedo of a planet, as defined by Bond, we have 

2 2 nm 
A= 2| d(x) sin « dz, 


0 


where 7 is the mean radius of the planet’s disk, R the distance of the planet from 
the Sun, A the distance of the planet from the Earth, M, the ratio of the apparent 
brightness of the planet at full phase to that of the Sun, 4(«) the ratio of the 
brightness at phase « to that at full phase. Then if o is the angular semi-diameter 
of the planet seen from the Earth, 


R?A?/r? = R® cosec?c, and if we write p= R*M,A?/r’, 
p= M,R?* cosec? oc. 
Now if o, is the semi-diameter of the planet at full phase, in seconds of arc, and 


g is stellar magnitude at full phase, reduced in both cases to unit distance from 
the Sun and from the Earth, and if Gis the Sun’s stellar magnitude at unit distance, 


log p = 0-4(G—g) — 2 logsino,. 
Again, if we write q=2 [ 4) sinzdz we have 4=pq. Russell shows that for 
a number of likely forms of 4(x), g= 2°20 x 4(50°). Hence 
A =p x 2:20 x 6(50°) 


may be called the Russell-Bond albedo. Now 4(z) is the ratio of the brightness 
at phase « to the brightness at zero phase, so that if 


pes, 
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where yp is in stellar magnitudes per degree and « in degrees, we have 

log = 

log 4(50°) = — 20, 

log q = 0°342— 
To evaluate G, at each wave-length we use values given by Woolley and Gascoigne* 
for Sun—«CMa at four wave-lengths from 6209 A to 4510A and interpolate 
(or extrapolate) these, according to the interval in 1/A, to the wave-lengths now 
used. Magnitudes of xCMa relative to the standard stars in ‘Table I were 
determined in the main stellar programme (of Woolley, Gascoigne and A. de 
Vaucouleurs). By addition we have G, relative to the stars in Table I, and hence 
we can form G, —g, without raising any question about the zero of the Mt. Stromlo 
magnitudes. 

To calculate the Russell-Bond albedo we need a value of o, and we take 
o,=4"-70 following G. de Vaucouleurs. Then log sino, = 5*358, and 
log A, =0-4(G,—g,) — 204; + 9°627. 

Details of the calculation in each wave-length are shown in Table V. 


TABLE V 
Calculation of the Russell-Bond alhedo of Mars 


A 4050 4250 4550 4945 5430 5980 6360 

O—aCMa 24°84 —24:98 25°16 25°36 ~25°57 ~25°76 —25°88 
aCMa 144 -~ - 
0,G, 26:28 — 26°42 26°59 — 26°79 26°99 —27°18 — 27°25 
3; O'14 — 0°02 0°63 — 1°14 — 1-60 — 1:98 — 2°31 

—26-42 —26-40 —25:96 25°65 25°39 —25'20 —24°94 
(G,—g;) 10°568 -10°384 -—10:260 10'156 —10:080 — 
20f3 — — — 0°474 — 0°430 6°230 0-268 
log A, 2°567 2-637 2-819 2°893 1°317 1°383 
A; 0'037 0'043 0:078 O'110 O'21 0'24 


Since Woolley and Gascoigne give +0™-04 as the standard error of their 
magnitude differences Sun—«xCMa, and since the s.e. of determinations 
Mars—a€ Ma is probably about +0"™-03 or perhaps +0™-04, the standard error 
in G;—g, is probably about +0™-05 or +0™-06. Accordingly the s.e. of 
0-4(G,—g,;) is about +0-02. Since the average s.e. of 201, is also about + 0-02, 
the standard error of log A, is about +0-03. The figure in the third place of 
decimals in log 4; shown in Table V should therefore be rounded off, and 4, 
is shown to approximately the right accuracy. 

Mars, at all events at zero phase, has quite a good gradient: this is to say that 
when the magnitudes are plotted against 1/A the points lie close to a straight line. 
Calculations are shown for «=0 and «=30° in Table VI. At zero phase the 
residual at A4250A is outstanding, and at 30° phase the fit is generally poor. 
There is no reason at all why the magnitudes of Mars should fit a gradient, and 
the success of the fit at zero phase is a matter of accident. The gradient-magnitude 


relation is of course 
mj, 


and is on the Greenwich system, since the Stromlo standard magnitudes have 
been adjusted to this system. 
* R. v. d. R. Woolley and S. C. B. Gascoigne, M.N., 108, 491, 1948. 
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G = 2-42 corresponds to a K-type star (e.g. 9Oct, Ki, G=2-39; «Aps, K5, 
G=2-47) but G=2-84 is very red, the only object redder than this yet observed 
at Mt. Stromlo being o Pav, Mo, G=2-88. The gradients here quoted are from 
A5980 A to A 4480 A, observed with the coarse grating by Woolley, Gottlieb and 
Przybylski.* 

TABLE VI 
Gradient of Mars 


Wave-length Zero phase 30° phase 


O Cc O Cc 
Magnitude,g,; G=2-42 O-C Magnitude G=2°84 O-C 


4050 A +0°14 +0:06 + 0°83 +0'05 
4250 A —0o'22 +0°20 ‘ +048 
4550 A —0'63 —0'63 0°00 0°00 0°00 
4945 A —1°14 —0°'05 —0'54 
5430 A —1°60 —1°§7 —0'03 —O'l4 
5980 A —1°'98 —2:'00 +o-o2 —1-62 —o'ol 
6360 A —2°31 —2:28 —0'03 —1°93 +0°'02 


This account of the monochromatic observations has been restricted to a 
plain recital of what was in the measures themselves, without reference to previous 
photographic and visual measures, or to the theory of the Martian surface. Good 
accounts of both of these will be found in Dr de Vaucouleurs’ book Physique de 
la Planéte Mars. A forthcoming English edition of this work will refer to the 
present measures in an appendix to the general discussion. 


It is hoped to repeat measures of monochromatic magnitudes of Mars with 
the same apparatus on Mount Stromlo at the next opposition. 


In conclusion I wish to thank Dr G. de Vaucouleurs for calling my attention 
to the desirability of making these measures: I have found the reduction very 
interesting. I also wish to thank Mrs A. de Vaucouleurs for her help in measuring 
the plates and checking the arithmetic. 


Commonwealth Observatory, 
Mount Stromlo, 
Canberra, Australia: 
1953 August 4. 


*R. v. d. R. Woolley, K. Gottlieb and A. Przybylski, M.N., 112, 665, 1952. 
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ERRATA 
M.N., 113, No. 2: 
R. Hanbury Brown and C. Hazard, A radio survey of the Milky Way in 
Cygnus, Cassiopeia and Perseus. 


P. 113, line 32, for power emitted by the feeder, read power transmitted by the feeder. 


A. D. Thackeray, /dentifications in the spectra of Eta Carinae and RR Telescopii. 
P. 212, 5 lines from foot of page: For +25-0 km/s read —25:0 km/s. 
The wave-lengths in Table II are correctly reduced for direct comparison of the main 
spectrum with laboratory wave-lengths. 
P. 216, lines 8 and 35-37: The wave-lengths in the first column of the table should 
be corrected as follows :— 
For 4214°14 vead 4214°45 
4314°57 
4319°6 » 4319°61 
4326-0 4326-07 
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